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Povzetek 



V tern delu raziskujemo kako lahko razlicna teoreticna odstopanja od okvira Standardnega modela 
(SM) vplivajo na lastnosti razpadov kvarka top. Manifestacijo nove Hzike (NF) onkraj SM, katere 
energijska skala znatno presega skalo elektro-sibkih pojavov, parametriziramo v obliki vrste visje 
dimenzionalnih operatorjev in z uporabo metod efektivne teorije polja. 

V prvem delu obravnavamo NF, ki bi lahko vplivala na nevtralne tokove kvarkov top, ki sprem- 
injajo okus. Ti tokovi so v okviru SM zelo redki in v primeru, da bi jih v bodoce eksperimentalno 
odkrili, bi to pomenilo odkritje NF. Bolj konkretno se zanimamo za dvo-delcne razpade t — > qZ, 7 
in tro-delcne razpade t — > q£'^£~. V analizo dvo-delcnih razpadov vkljucimo popravke prvega 
reda kvantne kromodinamike. Preucujemo tako mcsanje operatorjev pod renormalizacijo, kot tudi 
koncne popravke matricnih elementov, vkljucno s tako imenovanimi "bremsstrahlung" procesi. 
Izkaze se, da so popravki kvantne kromodinamke, zlasti v razpadih t qj, lahko pomembni in 
vplivajo na interpretacijo eksperimentalnih meritev. Glavni motiv za obravnavo tro-delcnih raz- 
padov je povecanje faznega prostora koncnega stanja, kar pomeni, da imamo na voljo vccje stevilo 
kinematicnih opazljivk, ki nam lahko pomagajo razlikovati med razlicnimi strukturami NF v nev- 
tralnih tokovih kvarka top, ki spreminjajo okus. 

V drugem delu pa se zanimamo za NF, ki bi lahko vplivala na glavni razpadni kanal kvarka top 
t —7- Wb. Ta je tesno povezan z nabitimi tokovi kvarkov. Z vnosom novih operatorjev, ki spremenijo 
lastnosti tWb interakcije, v teorijo, ne vplivamo le na glavni razpadni kanal kvarka top. Indirektne 
posledice, ki se pojavijo v teoreticnih napovedih opazljivk v fiziki mezonov so lahko znatne, saj 
virtualni kvarki top v redkih procesih igrajo vodilno vlogo. V tem delu zato podrobno preucimo 
indirektne omejitve na NF v nabitih tokovih kvarka top, ki izvirajo iz fizikc mezonov B. Sele, ko 
odkrijemo intervale za parametre NF, ki so skladni z meritvami mezonske fizike, se osredotocimo na 
posledice v visje energijskih procesih, razpadih kvarkov top na masni lupini. Osrednje opazljivke 
nase analize so sucnostni delezi (ang. helicity fractions) bozona W, ki nastane pri razpadu kvarka 
top. Ker se kvantno kromodinamski popravki pri obravnavi teh opazljivk izkazcjo za pomcmbnc 

V analizi SM, take popravke obravnavamo tudi ob vkljucitvi NF. Iz eksperimentalnih meritev, ki 
za enkrat kazejo ujemanje s SM, lahko postavimo direktne omejitve na NF, ki jih primerjamo z 
indirektnimi. V tem se nam razkrije zanimiva povezava med fiziko kvaka top in fiziko mezonov B. 

KlJuCne besede: Razpadi kvarka top, nevtralni tokovi, kvark top in nova fizika, efektivne teorije in nova fizika, 
popravki kvantne kromodinamike, sucnostni dclezi bozona W , indirektne omejitve nove fizike 
PACS: 12.15.Mm, 12.38.Bx, 12.60.Cn, 12.60.Fr, 13.88.+e, 14.65.Ha 



Abstract 



We study possible theoretical deviations from Standard Model (SM) in top quark physics which 
alter the decay properties of the top quark. Using effective filed theory techniques we parametrize 
the effects of potential new physic (NP) of scales well above the electroweak scale in terms of 
effective operators. 

On one side we investigate NP manifestation in the form of flavor changing neutral current 
(FCNC) decays of the top quark, which are highly suppressed in SM and potential observation 
of which would undoubtedly signal the presence of NP. We examine the two-body t — > qZ, 7 and 
three-body t qt^t' decays. Our analysis of the two-body FCNC decays is performed at next- 
to-leading order (NLO) in quantum chromodynamics (QCD). We examine the effects of operator 
mixing under the QCD renormalization as well as the finite matrix element corrections along with 
appropriate bremsstrahlung processes. We find that the effects of FCNC operators mixing under 
renormalization can be substantial, especially in the case of t ^ qj decays and QCD corrections 
affect the way the experimental measurements are to be interpreted. In the three-body decays we 
aim to exploit the increased phase space of the final state by defining different types of observables 
which could help to discriminate between structures of the vertices governing the FCNC transition 
of the top quark. 

On the other side we examine possible deviations from SM predictions in top quark's main 
decay channel, which is governed by the charged quark current interactions. Introduction of higher 
dimensional operators that modify tWb interactions, however has additional consequences that 
have yet not been thoroughly analyzed in the literature. In particular, low energy observables of 
rare processes in B physics, where virtual top quarks and their charged current interactions play 
a dominant role, are expected to be affected. We therefore perform a detailed study of indirect 
constraints on the NP operator basis and only then turn to the study of the effects in the decays of 
on-shell produced top quarks. The observables of interest are the helicity fractions of the W boson 
produced in the main decay channel. Since for SM predictions of helicity fractions higher order 
QCD corrections prove to be crucial, we conduct the analysis of NP effects at NLO in QCD. We 
confront our predictions with the experimental measurements to obtain the direct constraints on 
NP further comparing them with the indirect constraints from low energy processes revealing an 
interesting interplay of top and bottom physics. 

Key Words: Top quark decays, flavor changing neutral currents, new physics in top quark physics, effective 
theory approach to new physics, QCD corrections, helicity fractions of W boson, indirect constraints of new 
physics 

PACS: 12.15.Mni, 12.38.Bx, 12.60.Cn, 12.60.Pr, 13.88.-|-e, 14.65.Ha 
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Chapter 1 

Introduction 



1.1 The Standard Model 

The goal of theoretical high energy physics is to mathematically describe phenomena occurring 
at lowest experimentally accessible length scales using as the main tool quantum field theory. 
Theoretical predictions are put to the test by some of the most sophisticated experiments in the 
world ranging from different kind of particle colliders, to satellites orbiting our planet. Continuous 
advances in theoretical insight and experimental techniques in the past century have led to the for- 
mulation of the Standard Model (SM), theoretical framework describing the content of elementary 
particles and their interactions with its formulation dating back to the 1960's [1]. 

SM is remarkable both in its simplicity and great predictive power which has been put under 
enemas scrutiny over the last century, most recently by the Large Hadron Collider (LHC). 

In this introductory section we aim to give a very brief overview of the main theoretical features 
of the SM. Our discussion is kept on a purely informative level providing however appropriate 
references to be consulted for many underlaying details. A bit more time is spent on the discussion 
of flavor aspect of SM, since the concepts encountered there prove to be crucial for our top quark 
studies. At the end of the section we try to motivate the need for theoretical explorations beyond 
the SM, since this is the frontier that we shall be crossing in our work. 

1.1.1 Particle content 

The corner stone of SM as a quantum field theory is the gauge group under which the theory is to 
be invariant 

SU{3)cX SU{2)lxU{1)y. (1.1) 

Here SU{3)c is the gauge group of quantum chromodynamics (QCD), SU{2)i is the weak isospin 
group and finally, the Abelian U{1)y is called the weak hypercharge. The corresponding coupling 
constants of the three groups are denoted hy Qs, g and g' respectively. Specifying the gauge group 
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entirely fixes the content of the gauge boson sector. On the other hand, there is more freedom 
in specifying the scalar and fermionic sector of the theory. In particular we have to assign what 
representations of the gauge group the fields are to be in. 

The only scalar field in SM is the Higgs boson. Its representations under the gauge groups can 
be written as </>(!, 2)_|_i/2) meaning that it is a singlet under SU{3)c, a doublet under SU{2)l and 
it carries a hypercharge of +1/2. At present it remains the only quantum of the SM that has not 
been experimentally confirmed, however the recently discovered new particle at LHC [2, 3] poses 
to be a very strong candidate. Higgs boson plays a very specific role in the SM, since by acquiring 
a vacuum expectation value (VEV) it spontaneously breaks the SM gauge group 

5[/(3)c X 5C/(2)l X C/(l)y ^ SC/(3)e X t/(l)Q , Q = Y + T3, (1.2) 

where U{1)q is the gauge group of quantum electrodynamics (QED), with Q the electromagnetic 
charge and T3 the eigenvalue of the diagonal SU{2)l generator. QED gauge coupling is e = 
gsinOw = g' cos 6w, where 9w is the Weinberg mixing angle. The described pattern of symmetry 
breaking gives rise to mass terms for the three weak gauge bosons through the covariant derivative 
(Dfj^c/))^ (D^cf)) part of the Higgs Lagrangian and to mass terms for the fermions through the so-called 
Yukawa interactions, which we shall return to shortly. 

Turning to the fermionic sector of the SM, we first note the subscript L of the weak gauge 
group, which describes an important postulate of the SM, that only left-handed fermions carry 
weak- hypercharge . 

Further we classify fermions depending on what representation of SU{3)c they are in. Singlets 
are called leptons and quarks are postulated to be in its fundamental representation. 





Quarks 


Leptons 


left handed 
right handed 


ur{3, 1)+2/3, dR{3, l)_i/3 


^l(1,2)_i/2=(;J) 
^r(3,1)-i 



Table 1.1: SM fermions and their representations under the SM gauge group (1.1). Subscripts L 
and R denote the chirality of the fields. 



The fermionic sector of the SM is further enriched by making three repetitions (generations) 
of gauge representations described above and given in Tab. 1.1. Each of the repetitions is assigned 
a flavor allowing us to distinguish among them. We define {u, c, t} up-type quarks, {d, s, b} down- 
type quarks, {e,/i, r} charged leptons and the accompanying neutrinos {z/g, i^^, p't}- Since it is not 
crucial for our studies, we shall not deal with the leptonic part of flavor physics and concentrate 
only on the quarks. 
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1.1.2 Flavor 

The only part of the SM Lagrangian that is not flavor universal, meaning that it distinguishes 
between difl^erent flavors, is the Yukawa interaction term 

Cy = -Qi[Yd]ijct><ffj^ - QWuWWr + h.c. . (1.3) 

Indices i and j denote the flavor and we have introduced the "up" and "down" 3x3 complex 
Yukawa matrices and = \o^(\). While there are two Yukawa matrices for the quarks, there is only 
one for the leptons, because as evident from Tab. 1.1, there are no right-handed neutrinos in the 
SM. 

In the SM the Yukawa sector is the only source of flavor physics. This statement can be put 
in a more formal group theoretical form by saying that Yukawa interactions break the big global 
symmetry of flavor 



[/(3)q X C/(3)„ X C/(3)rf , (1.4) 



obtained when three generations of fermions are introduced to the theory. V(^)q^u4 are groups of 
rotations in flavor space Vq^u^d that can be applied to Q, n and d quark fields respectively 

Ql > Ql = vqQl , UR > uji = VuUR , dR > d^ = VddR , (1.5) 

where we have suppressed the flavor indices. Omitting the scalar fleld, a general rotation 
effects the Yukawa terms in the following way 

QLYddR = Q'l VgYdV} 4 , QlYuUr = Q'^ VgYaV^ u'r , (1.6) 

which can be seen as a change of basis in the flavor space. We can use the rotations of the broken 
symmetries to rotate away all the unphysical parameters of the Yukawa sector, since we know that 
out of all the parameters, there are as many unphysical as there are generators of broken symmetries 

-^phys. -^all -^broken gen. • 

We note that there is a remaining U{\)b global flavor symmetry even after inclusion of Yukawa 
terms, associated with the baryon number conservation. This means that we start of with 36 (18 
real and 18 imaginary) free parameters of Yukawa matrices and break 26 out of 27 generators of the 
global symmetry. 17 of which are rotations containing phases, and 9 are rotations with no phases^ 
this leaves us with 10 physical flavor parameters, 9 of which are real and 1 is a complex phase. 
Using Vq, 

,u,d we can rotate quark flelds to the basis where one of the Yukawa matrices is diagonal 
Q'l VgYdVj 4 + Q'lVqV^\W^u^ . (1.7) 



^U{3) can be written as U{1) x SU{3). The (7(1) transformation is rephasing e'^. SU{3) transformation can be 
written as e'""^ , where T'^'*''^ are imaginary and contribute real rotations, while T^'^'*'* are real and contribute 
rotations with rephasing. 
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In the second term we had to insert a unitary matrix Vx to obtain an expression with a diagonal 
Yu, which is however multiphed with a unitary matrix V, called the Cabbibo-Kobayashi-Maskawa 
(CKM) matrix, which was formulated in [4, 5] and will be subject of more discussion later. This 
basis is usually referred to as the "down basis", while we could have performed an analog diago- 
nalization of the Yu matrix, ending up in "up basis" . 

When the SU{2)l x U{1)y gauge symmetry is broken by the VEV of the Higgs field ((/>) = 
(0,f/\/2), Cy can be written as 

Cy = ~d'^Yt^-d'^ - ^u'^VYf^^-u'j, + .-.. (1.8) 

Note that we have obtained a mass term for every down-type quark. Since SU{2)l has been broken, 
we can rotate the ul fields with a separate rotation. By choosing the rotational matrix to be 
the CKM matrix = u'^^V, we obtain the mass terms for the up quarks as well. This rotation has 
important consequences in the charged current sector, where it generates flavor changing currents. 
Removing all the primes from the quark fields and reintroducing the flavor indices the charged 
current Lagrangian is 

^cc = [u^Ll''d,L] V,,W^ - -| [d,Ll''UiL\ V^.W- . (1.9) 

On the other hand, due to the unitarity of CKM matrix V ^ there are no tree- level Eavor changing 
neutral currents (FCNCs) in SM^. The elements of the CKM matrix are usually denoted as 



V = I Vcd Vcs Vcb , (1.10) 




where not all of the parameters are independent. We have accounted for 6 of the flavor parameters 
to be quark masses, meaning that the CKM matrix must contain the remaining 4 of which 3 are 
real and 1 is a complex phase. 

The CKM mechanism has proven to be very successful in describing flavor physics which has 
been tested with various high precision experiments. The so called CKM unitary triangle has been 
over-constrained by numerous measurements and is showing remarkable consistency. Furthermore 
a highly diagonal dominated form of the matrix has been experimentally established and perhaps 
most importantly, the complex phase has been measured to be non-zero, proving that the discreet 
symmetry of simultaneous charge conjugation and parity (CP) is indeed violated in nature. The 
CKM phase is the only source of this violation within SM. For in-depth coverage of the subjects 
on CKM mechanism and CP violation we refer the reader to the following references [6, 7, 8]. 

1.1.3 Need to go beyond SM 

Despite the unprecedented success of SM it is evident that it cannot be the final theory of elementary 
particles and their interactions. The discovery of dark matter (see for example [9, 10]) and neutrino 



^There are no dimension 4 operators in SM that would generate FCNC transitions. 
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oscillations [11] show that the particle content of SM is not adequate since there is no dark matter 
candidate and neutrinos are massless in SM. While the CP violation at low energies is very well 
described by the CKM mechanism, there are strong indications, that far more violation is needed 
to explain the high dominance of matter over anti-matter that we observe in the universe today 
[12, 13, 14]. 

What is more, SM does not try to describe gravitational interactions which become relevant 
at very high energies, of the Planck scale Ap = y^hc/G^ ~ 10^^ GeV. Assuming that SM could 
be a valid theory all the way to the Planck scale gives rise to a puzzling situation referred to as 
the "hierarchy problem" [15, 16]. The word hierarchy is related to the large separation between 
the electroweak and the Planck scale. Due to the fact that the dimensionality of the Higgs mass 
operator is 2, the radiative corrections to the mass dniH turn out to be quadratically proportional 
to the cutoff scale A 

ml = mT-j^,A' + ---. (1.11) 

Under the assumption that A ~ Ap, a great deal of fine-tuning (~ 30 orders) between the parameters 
entering the expression is necessary in order to obtain the appropriately low Higgs mass (of 
the electroweak order). 

The different views about its significance as a problem of the theory has earned the hierarchy 
problem to be labeled as somewhat controversial. Nevertheless it is at least an aesthetic issue 
and one that beyond Standard Model (BSM) theories (often in the literature referred to as UV 
completions) try to address and do so in many different manners. 

Following the aesthetic drive, the idea to unify interactions by embedding the SM product 
gauge group (1.1) into a single larger group comes naturally. The so called grand unified theories 
aim to do just that [17]. 

Because the flavor parameters, masses as well as mixing parameters described by the CKM 
elements, exhibit strong hierarchy one can not help but wonder if there is an underlaying symmetry 
manifested at energies above the electroweak scale that could explain it. SM provides no answer 
to this question which has become known as the flavor problem of SM. 

There is clearly more than enough reasons for theorists to explore the possibilities of new physics 
(NP) BSM and further to look for observables and processes which could help to discover NP and to 
discriminate between different NP scenarios. In many of them, the top quark provides a preferred 
search window due to its large coupling to the physics responsible for the electroweak symmetry 
breaking. A fascinating possibility, that we shall be exploring in our work, is that the top quark 
properties exhibit deviations from their predictions within the SM. 

In the next section we argue that among other places, also the top quark decays are good 
"hunting grounds" for NP. In particular decay modes and observables which are within the SM 
predicted to be highly suppressed are promising for detection of NP effects, since typically a non- 
zero measured signal of such quantities would, right-away, present a signal of physics beyond SM. 
Important thing to keep in mind is that in addition to directly probing top quark physics at LHC and 
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Tevatron, top quark properties can also be explored in lower energy phenomena of meson physics, 
where top quark appears as a virtual particle often having the leading role in rare processes. 

In our exploration of BSM top quark physics, we shall not be committing to particular NP 
models or frameworks. Rather we will take the effective field theory approach, adding to the SM 
effective Lagrangians with which will be able parametrize our ignorance about the BSM theory. 
We will then study different observables which might be sensitive to our additions. In principle our 
results may be applicable to a variety of BSM models. 

1.2 Top quark decays 

Top quark is the heaviest experimentally confirmed elementary particle. It was discovered at the 
Tevatron in 1995 [18, 19]. The two main features of the top quark are its large mass, experimentally 
measured to be [20] 

= 173.2 ±0.9 GeV, (1.12) 

and its decay width. In the SM top quark is predicted to decay almost exclusively through the 
charged weak current (1.9). The t — )> Wb channel, which we shall refer to as the main decay 
channel, is highly dominant due to the extreme hierarchy between the CKM matrix (1.10) elements 
of the third row \Vtd\, {Vtsl <C \Vtb\. Branching ratios of top quark decays are always normalized to 
the main decay channel. The tree-level decay width computed at leading order (LO) in QCD can 
be written as 

r(*-l^^) = |V„|^^f^i^^;lii^~1.5GeV. (1.13) 

where x = mw/rrit, m\Y = 80.1 GeV, the mass of b quark has been neglected and \Vtb\ = 1. 
Numerical value of g is related to the Fermi constant through Gpjypl = 5^/(8m^) = 1.167 GeV~'^ 
[21]. Due to its large mass, the average life time of top quark is an order of magnitude below 
the typical hadronization time scale, causing it to decay before forming bound states [22] making 
theoretical treatment free of non-perturbative QCD effects. 

While the production of top quarks is a very interesting area of research as well, especially in 
the case of ti production, where we are witnessing a persisting anomaly in the forward-backward 
asymmetry from Tevatron [23, 24] that has stimulated various theoretical attempts to reconcile 
it [25, 26], we shall be concentrating in this work on two aspects of top quark decays that we describe 
below within the framework of SM showing that they might be interesting for NP observations and 
constraints. On one hand we will study the main decay channel exploring the helicities of W bosons 
produced through the top quark decay. The information on what fraction of Ws produced in the 
decays have certain helicities allows us to directly probe the structure of the tWb interaction and its 
potential deviations from the SM form (1.9). On the other hand we shall consider the possibility 
of observing FCNC decays of the top quark. The branching ratios for these decays are highly 
suppressed within the SM and any observation of such a process would signal a presence of NP. 
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Both analysis, which are given in sections 2.3 and 3.5 wih be conducted at next-to-leading order 
(NLO) in QCD. 

We should note that in the last five years, the precision of the experimental top quark mass 
determination has been gradually improving and the central value (1-12) of the top quark mass 
has been continuously changing. As a consequence we will, in this work, encounter a few different 
values for top quark mass being used in the numerical analysis since the work presented here spans 
over four years. Deviations are however small and variations within these values have no significant 
effect on the results presented. 

1.2.1 Helicity fractions in the main decay channel 

Since W boson is a spin 1 particle, we can split the decay width of the top quark's main decay 
channel into three parts depending on which of the three helicity states the produced W is in 

T{t^Wb) = TL + T+ + r_, (1.14) 

where L stands for longitudinal, while + and — denote positive and negative transverse helicities 
respectively. We further define the helicity fractions as 

•^L,+,- = ^^, (1.15) 

telling us what fraction of W bosons produced in top quark decays have certain helicity. The 
main reason why helicity fractions are interesting for NP searches is that they are sensitive to the 
structure of the tWb vertex governing the decay. 

On the computational side, there is more than one way to extract a certain helicity of the final 
state vector boson. We shall be making use of the covariant helicity projectors [27], which are 
particularly useful when computing loop diagrams for QCD corrections. To define them we write 
down the squared matrix element for the t — >• Wb decay as 

\M\^ = H,,e^^{q,X)e*''{q,X). (1.16) 

€^{q,X) are the polarization vectors of the W fields with A = 1,2,3 labeling their basis and q 
denotes the momentum of the W. We have put everything else into the H^^. When going from the 
squared matrix element to the decay width we can, even when considering particular helicity final 
state, perform the summation over the polarizations of the W boson X^;^ e^((7, A)e* (g. A) replacing 
it with appropriate helicity projector given in Tab. 1.2. In the last column we show the LO helicity 
fractions in the limit of the massless b quark. We can see that within SM helicity fraction is 
in the presented approximation. The suppression is not hard to understand and is illustrated in 
Fig. 1.1. If we consider b quark to be massless and produced in the weak interaction, which strictly 
involve left-handed components of the fermionic fields, its helicity has to be negative, since for 
massless fermions helicity and chirality coincide. From a simple consideration of spin conservation, 
the situation where W boson would have a positive helicity is not possible. 
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Helicity 


Projector: Yl\ 




, A) W 


SM LO J='i with mb = 


Unpolarizec 


: r 


P^^ = 




m'i,. 






Asymmetric: 


Tf 


P^^ = 


mt |q| PtaQ/B 






Longitudinal: 


Tl 




'^W 1 
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Table 1.2: Covariant projectors extracting different helicities of a final state massive vector boson 
in a three-body decay. Presented projectors are for t — )• Wb decay, where pt is the momentum of 
the decayed top quark, and q is the momentum of the W. The three- vector length |q| is assumed 
in the top quark rest frame and eoi23 = 1- Appropriate projector Ff'^ is to replace the sum over 
polarization vector basis depending on what helicity state we wish to project out. Last column 
shows the tree-level SM helicity fractions at LO and in the limit of the massless b quark. 

Figure 1.1: An illustration of t — >■ Wb decay in the rest frame of the 
top quark where the limit nif, = is taken. Wide arrows represent 
third component of the spin with respect to the horizontal axes, while 
narrow arrows represent the direction of momentum. Because helicity 
and chirality of the massless b quark coincide, its helicity is always 
negative. Third pictures represents a situation that is not possible since 
top quark would have to have spin greater than 1/2 to accommodate 
the spin conservation, indicating the helicity suppression of 

This simple picture is altered if we consider the mass of b quark or the process involves more 
than just three particles, which is the case once higher order quantum corrections to the decay are 
considered. All of these effects have been analyzed within the SM including m;, and finite width of 
top quark effects, NLO QCD and electroweak corrections as well as next-to-next-to leading order 
(NNLO) QCD corrections [28, 29, 30, 31]. We summarize the theoretical prediction to be 

J-SM ^ 0.687(5) , Tf^ = 0.0017(1) . (1.17) 

Even with inclusion of these corrections J-+ remains highly suppressed, and a measurement of 
the positive helicity fraction of the per-cent order would undoubtedly signal the presence of NP. 
How these predictions get altered by the presence of NP governing the t — )■ Wb decay, is the subject 
of section 3.5. 

Helicity fractions are experimentally accessible through the measurement of angular distribu- 
tions of the leptonic final state to which the W boson decays t — )• Wb — )• bii' as indicated in Fig. 1.2. 
The extraction is based on the distribution of decays over cos 9* , where 9* is the angle between the 
momentum of the charged lepton in the W boson rest frame and the momentum of the W boson 
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Figure 1.2: Leptonic final state of the top quark's main decay channel used for the extraction of 
helicity fractions. Angle 9* is defined as the angle between the direction of the lepton in the W 
rest frame and the direction of the W in the rest frame of the decayed top quark. 



in the top quark rest frame [32] 

dcost/* 4 8 8 

At the moment the most precise measurements of the helicity fraction still come form the 
combined DO and CDF analysis [33] 

jrCDF+DO = 0.722 ±0.081, = -0.033 ± 0.046 , (1.19) 

which are for now showing agreement with the SM predictions (1.17). Helicity fractions are also 
being measured at the LHC [34] and the sensitivity expected to be reached is [35] 

a{T+) = ±0.002 , a{TL) = ±0.02 , (1.20) 

where cr{J^i) denotes the predicted absolute error on the measurement of the helicity fraction with 
10 fb~^ of accumulated data. 

1.2.2 FCNC top quark decays 

As mentioned in section 1.1.2 there are no FCNC decays possible at tree- level in the SM. They can 
however occur at one-loop level where through two insertions of charged flavor changing currents 
we can obtain a neutral flavor change. 

Top quark FCNC decays that we will be interested in t — )• qV, where V = ^, Z, g and q = u,c, 
proceed through the so-called penguin diagrams shown in Fig. 1.3. If we choose to work in the 
unitary gauge for the electroweak interactions, there are indeed only two diagrams to consider. By 
computing these diagrams we obtain the decay widths of the following form 

^it^qV)^\Y,>^^^f^''H^^)\\ (1-21) 

i 

where i denotes the flavor of the down-type quark with mass rrii running in the loops and Xi = 
mf/m'^. The loop functions /^-^^ are specific for the V gauge boson in the final state. Further we 
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Figure 1.3: Feynman diagrams contributing to FCNC decays t — )• (77, Z,g in the unitary gauge for 
the weak interactions. Quarks running in the loops are down- type and their flavor is denoted by i. 



have defined 

X'f'^ = V:,V,,. (1.22) 
Note that due to unitarity of the CKM matrix, the following equation holds 

+ Ai«^ + a['^ = , (1.23) 

since there is an orthogonal relation between different rows of CKM matrix (1.10). Consequently, 
when computing the decay widths (1.21), the sum over the flavors enables us to drop any term 
in /(^) that is not dependent on Xj. This is a property of loop functions often encountered when 
dealing with one-loop induced FCNC processes. Neglecting the ma^s and making use of unitarity 
relation (1.23), the branching ratios for these decays are proportional to 

BT[t^qV]^\V,k\'\f'^''Hx,)\\ (1.24) 

The suppression of the branching ratios is two fold. Firstly Xb <^ 1, so the loop functions give small 
contributions and secondly \Vgh\ ^ 1. The resulting values for t — )• cV are [36, 37] 

Br[t C7] ~ 10"^^ , Br[t cZ] ~ 10"^"^ , Br[t eg] ~ 10"^^ , (1.25) 

while the results for t — )• uV are an additional order of magnitude smaller. Within many BSM 
models, like Two Higgs Doublet Models, Minimal Super Symmetric Standard Model, models with 
up-type quark singlets, etc., the suppression of FCNC top quark decays can be lifted [37, 38, 39, 40]. 
It has been pointed out recently, that top quark FCNC phenomenology is crucial in constraining a 
wide class of NP scenarios where new flavor structures are present but can be aligned with the SM 
Yukawas in the down sector [41, 42, 43, 44]. 

Top quark FCNCs can be directly probed both in production and in decays of the top quark 
and all three FCNC decays presented here are being experimentally searched for. There has been 
no observation made so far, thus upper limits at 95% confidence level (C.L.) have been set, most 
stringent of which are 

Br[t ^ {u,q}'y] < {5.9, 32} x 10-^[45, 46] , (1.26) 

Br[t qZ] < 3.4 x 10"^ [47] , 
Br[t^{u,c}g] < {5.7, 27} x lO'^^S] . 
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With the increasing accumulation of data LHC wiU be able to probe branching ratios of lower 
orders [49], in particular in the case of no signal, ATLAS projects to improve the upper bounds 
(1.26) to 

Br[t ^ C7] < 10-5 , Br[i ^ cZ] < 10"^ . (1.27) 

1.3 Effective field theories 

The concept of effective field theories is highly applicable in high energy physics where we often 
encounter problems involving widely separated scales. 

On one hand effective theories are very useful when the underlying theory is unknown and 
allows us to parametrize its effects on the physics at lower energies in a systematic fashion. On the 
other hand, it is also useful when the underlying theory is known, since in general the full theory 
can be quite complicated and going to an effective theory simplifies matters greatly. In particular, 
going to an effective theory can manifest approximate symmetries that are not visible in the full 
theory and increased symmetry means increased predictive power. 

Furthermore, when the full theory contains several disparate scales m <^ M, perturbation 
theory can be poorly behaved as typically, when considering higher order quantum corrections, one 
generates logarithmic terms of the form log(m^/M^). When these logs are large, they need to be 
re-summed in a systematic fashion in order to keep perturbation theory under control. Working 
within an effective theory simplifies the summation of these logs. 

In this section we aim to briefly introduce the effective theory techniques that we shall be 
employing in our analysis. It is worth mentioning that the strength and applicability of effective 
theories in particle physics goes far beyond what we shall be presenting. For detailed explanations 
we refer the reader to the following pedagogical works [50, 51, 52, 53]. 

1.3.1 Operator product expansion 

The operator product expansion (OPE) [54] translates a time-ordered product of two operators to 
a series of local operators 

m(x)02(y) ^ j;Cp)(x-y)0,(x), (1.28) 

i 

where the spatial separation x — y is assumed to be small. Wilson coefficients, which are c- numbers, 
capture all the short distance x—y dependance. We can apply this expansion when we are computing 
amplitudes for different processes, which is typically done in the momentum space. In particular, 
when we encounter a Feynman diagram with a virtual heavy particle of mass M and we are 
interested in the external momenta p, where p <C M we can perform a Taylor expansion of the 
amplitude A in the parameter p/M. We can then ask ourselves what kind of effective Lagrangian 
which does not include the heavy field we would need to write down in order to be able to reproduce 
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the full theory result. This leads us to an expansion of the form [53] 



id) 1 




(1.29) 



i.d 



where the sum runs over d, representing the dimensionality of the local operators, and over i, the 
basis of operators with a given dimensionality. Typically the basis consists of more than just one 
operator, and each operator comes accompanied with its own Wilson coefficient, which is in this 
equation assumed to be dimensionless. OPE reveals a very important point, that contributions of 
higher dimensional operators come suppressed with higher powers of the high scale M, representing 
the scale of the physics that has been "integrated out" . 

The procedure of integrating out the heavy fields reduces the number of dynamical fields in the 
Lagrangian. If the underlying theory is known, we are able to match the effective and full theory, 
thus obtaining the Wilson coefficients. If, on the other hand, the full theory is not known, the 
matching procedure can not be performed, but we can still rely on the OPE to write down an ap- 
propriate basis of operators and analyze how they affect certain observables and most importantly, 
truncate the series at a certain dimensionality of the operators, knowing that higher dimensional 
operators are accompanied with higher suppression as indicated in (1.29). 

1.3.2 Running of the Wilson coefficients 

As already mentioned, computation of quantum corrections within a theory containing two scales 
m, M (which we will refer to as the full theory) will typically introduce the following form of 
logarithmic terms in the amplitude 



where a and h in general denote a product of different couplings. We assume a to come from LO 
diagrams, while the h term comes from some NLO corrections. For the NLO part we have written 
out only the logarithmic term, making our analysis a leading-log approximation. Finally, (Oj) is the 
matrix element of a certain operator. If the two scales are widely separated m <^ M, the logarithm 
of their ratio can be big and we might encounter a problem, since even if b is small, comprised 
of parameters appropriate for perturbative expansion, the logarithm creates a potentially large 
enhancement, rendering the perturbative expansion at least questionable. 

The best way to resum the large logs is to employ the effective theory approach. We compute 
the amplitude for the same process and to the same order in perturbation theory in the effective 
theory £eff = CiOi, from which the heavy degree of freedom (having mass M) has been integrated 
out. The amplitude will be UV divergent and the factor in front of the divergence will exactly 
match the factor in front of the large log in the full theory 




(1.30) 



AeS = --- + Ci(l + -{--log 
V a e 



m 



2 



2 




7 + log 47r . 



e 



e 
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Here we have chosen the dimensional regularization [55, 56, 57] of the UV divergence, working 
in d = 4 — e dimensions, which necessitated an introduction of an arbitrary scale fj, and 7 is the 
Euler constant. Performing the matching procedure between (1.30) and (1-31) order-by-order in 
perturbative expansion, thus gives us the Wilson coefficient 

Ci = a + 6(^- + log-^j + ••• , 

which is UV divergent. We can renormalize it using the modified minimal subtraction (MS) renor- 
malization scheme, very appropriate for the renormalization group applications (for detailed dis- 
cussion see [52]). Notice, that had we performed the matching only at LO, the extracted Wilson 
coefficient would have been the same as the leading-log coefficient, with the matching scale is set 
to n = M. 

The renormalization of the effective theory is said to involve operator renormalization, stating 
that either the operators or the Wilson coefficients in the Lagrangian are bare objects, denoted by 
(0), and need to be renormalized by the appropriate renormalization matrices denoted by Z 

Bare operators: QOf ^ = CiZijOj , (1.32) 

Bare W. coefficients: cf = dZflOj = dZfJ ZJ^O^^'^ . (1.33) 

Comparing the last expression of the second line with the first expression of the first line we find 
the relation among the two renormalization matrices 

2"=^ = Z-'^ . (1.34) 

The extraction of Z in the minimal subtraction scheme is achieved by finding the UV divergent 
parts of the effective theory amplitude. Operators are composed of fields and other parameters 
which themselves need to be renormalized. Usually we consider these objects appearing in the 
operators to be renormalized, so when looking for the operator renormalization matrix, this needs 
to be taken in consideration (for details see Ref. [52]). To obtain renormalization group equations 
(RGE) for the Wilson coefficients we use the fact, that the bare quantities cannot depend on 



d ,^(0) ^ d 



d 



H—(Dl">=0 =^ fi—0 = -Z-'fi—ZO = -^0, (1.35) 
d// dfi I d/i 



d/i . 



T 



C = j^C. (1.36) 



We have defined the anomalous dimension matrix 7 which governs the running of Wilson coefficients. 
It is dependent on some perturbative coupling, which we generally denote g, and is in the cases 
that we will be considering always QCD coupling constant gs- Important thing to note is that the 
coupling constant is itself dependent on /i and its running is determined by its beta-function 

M^ = /3(9), (1.37) 
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where QCD beta-function is known to 4 loops [58]. Taking that into account, the solution to the 
differential equation (1.36) can be expressed in the following way 



We can now depending on which scale we want to evaluate the matrix element of the operator at, 
run the Wilson coefficient to that scale and by doing that re-sum the large logarithms promoting 
our result to renormalization group improved perturbation theory prediction. 

1.4 Top quark in meson physics 

Top quark plays an important role also in the physics of energies lower than its mass. In such 
processes there is not enough energy to produce an on-shell top quark, rather it appears as a 
virtual particle in loop diagrams. Similarly as the existence of the charm quark was predicted by 
the Glashow-Iliopoulos-Maiani mechanism [59] before its experimental discovery, there was a strong 
belief in the existence of the top quark well before its discovery at the Tevatron. What is more, 
due to good theoretical and experimental understanding of rare processes in kaon and B meson 
physics, its mass was well estimated [60, 61, 62, 63]. 

As pointed out in Ref. [41] when searching for BSM physics in top quark sector, which we have 
set out to do, one should also consider the meson physics observables and the indirect effects that 
NP in top quark physics might cause. 

Employing the OPE and effective theory techniques presented in section 1.3, our study of 
such indirect effects is reduced to finding the Wilson coefficients of lower energy theory where the 
top quark and the heavy vector bosons have been integrated out. Our modification of the SM 
Lagrangian will impact only physics of scales at which QCD is perturbative, and will therefore 
all be contained in the Wilson coefficients. Once we shall compute the Wilson coefficients, the 
procedure of obtaining low energy meson observables will be exactly the same as in SM since no 
modification to physics of low energies is made. 

In this section we go to some detail in explaining the matching procedure at LO in QCD for 
|Ai?| = 2 and |Ai?| = 1 processes, in which top quark turns out to play an important role. The 
same computational approaches that we introduce here for the SM case shall be employed when 
we consider NP manifestation in charged quark currents in chapter 3. 

1.4.1 |Ai?| = 2 transitions 

Mixing between Bq and Bq mesons, where q stands for either d or s down type quarks, is a \AB\ = 2 
process since a 6 o 6 transition occurs. It is a FCNC process highly suppressed in the SM and 
sensitive to NP effects. For pedagogical description of theoretical treatment as well as insight into 
the experimental aspects of meson mixing we refer the reader to the following references [6, 7, 64]. 





(1.38) 
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The Bq and Bq states are flavor eigenstates and they oscillate between each other. Within the 
Wigner-Weisskopf approximation, the oscillation is governed by the Schrodiner equation 

where M*^ and are Hermitian mass and decay matrices. The physical eigenstates \Bl) ,\Bh), 
having masses and T^^ decay widths, are obtained by diagonalizing the M'^ — ir''/2 matrix. 

The oscillation (1.39) between the flavor eigenstates involves three physical quantities 



Mj 

12 



\M1,\, |r?2|, ,^^ = arg(-^), (1.40) 



which are the off-diagonal mass and width matrix terms and the CP phase respectively. On the 
computational side, Mj^2 is obtained from the dispersive part of the transition amplitude between 
the meson and anti-meson, 

M!, = ^^{Bq\n[^^^^='\BqU.,, (1.41) 

where is the mass of the B meson. On the other hand is obtained from the absorptive 

part of the same matrix element, but we shall not be considering it further. 'H^^^^~'^ = 
is the effective Hamiltonian governing the |Ai?| = 2 transitions and in the SM we have 

^|AB|=2 ^ -^|^(F;^1/,,)2Ci(/.)0? , Of = [q.j^b,] [q,j,b,] . (1.42) 

In order to find the Wilson coefficient at high scale of the W boson and top quark mass, which we 
denote as nw, we need to perform the matching at LO in QCD. 
Since the quark content of B and B mesons can be written as 

Bq^bq, Bq^bq, (1.43) 

we need two q and two b field operators in order to contract the final and initial states and obtain 
a non-zero matrix element. On the full theory side, this necessitates a fourth order perturbative 
insertion of charged current interactions given in Eq. (1.9). For simplicity we shall make use of the 
unitary gauge for weak interactions, eliminating the would-be Goldstone fields from the theory. At 
the lowest order the mixing of B mesons proceeds through the so-called box Feynman diagrams 
which are presented on the left side of Fig. 1.4. 

In computation of the box diagrams the momenta of external quarks are set to zero, since 
the operator we are matching to contains no derivatives or masses [53]. This leaves us with all 
the propagators in the loop having same momentum and greatly simplifies the integration. Under 
these assumptions both diagrams contribute equally and the obtained amplitude reads 

. 4 u,c.t 

^fuii = 'I 2 [qlLYh:iL][q2Ll^ML\ + [giL7''&4L] feiTA^^ai] \hq) V*gVaV,hV*qF{x„x,) , (1.44) 
D47r . 
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Figure 1.4: Feynman diagrams for Bq — )• Bq transition at LO in QCD for the full theory (left) and 
in the effective theory (right). Quark running in the loops are up-type and i,j denote their flavor. 



where indices 1, ... ,4 on the quark fields label the contractions of the external fields. F{xi,Xj) is 
an Inami-Lim function first presented in Ref. [65] with Xij = mfj/my^. When computed in the 
unitary gauge it is UV divergent 

F{x,,Xj) = [-_-\og^)^^ + •••, (1.45) 

where we have written out only the manifestation of the UV divergence within dimensional reg- 
ularization, denoting the rest with dots. What renders the amplitude (1.44) finite are the two 
summations over the flavor of the up-type quarks running in the loops. In a similar fashion as in 
Eq. (1.22), we define 

>^? = y:,y^b, (1.46) 

where now the role of up-type and down-type quarks are interchanged and the unitarity of CKM 
matrix yields, in analogy to Eq. (1.23) the relation 

a('?) + a(«) + a1'')=0, (1.47) 

which tells us that we can safely drop any term in F(xj, Xj) that is not dependent on both Xi and 
Xj. Consequently, the terms presented in Eq. (1.45) get cancelled. 

Neglecting the masses of light up quarks x^j = Xc = and making use of (1-47) and for brevity 
suppressing the {q) superscript in A^-''^ we arrive at the simplified result 

^fuu = '^^{k\ [qiLl^hL][q2Ll,b,L] + [qiLl''b,L][q2Ll,hL] mX^tSl'^ixt) , (1.48) 

where 

Sl^ix,) = F{x„ X,) - 2F(x„ 0) + F(0, 0) = "^^^(x/- 1^^ + ' ^'-^'^ 

This results explicitly confirms what we have been stating about the dominant role of the top quark 
in the mixing process. To complete the matching procedure we have to calculate the amplitude for 
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the same process using effective tfieory (1.42), where we have just one local operator to be inserted 
at first order of perturbation. We can choose any of the two b fields in the operator to contract 
the b quark final state, giving us a factor of 2. We then have a further choice of contracting the q 
fields of the operator with the external q quark fields, corresponding to the two Feynman diagrams 
given on the right side of Fig. 1.4, obtaining the result 

Aes = i^^|^A?C7i(6g| [qili'^^sl] feL7/.fe4L] + [qiL7''b4L][q2L-/^.bsL] \bq) . (1.50) 

Comparing (1.4S) and (1.50) we find for the Wilson coefficient 

C!'^{fiw) = Sl'^{xt). (1.51) 

We have demonstrated the matching at LO in QCD, however higher order perturbative QCD 
corrections are not negligible in such processes and can be captured in rescaling factor r]B [64, 66] 

^SM.NLO ^ (^SM,LO^^ ^ = 0.8393 ± 0.0034 . (1.52) 

Furthermore, considering NLO QCD corrections and employing RG methods the coefficient can 
be run down to lower energy scales at which different non-perturbative methods can be used 
to evaluate the matrix element of the Of operator. Just for completeness we show the typical 
parametrization of the matrix element following Ref. [64] 

{B,\Ol{f,b)\B,) = \mIjIBbM) , (1.53) 

where and Bb^ are nonperturbative parameters, the decay constant and the bag parameter 
respectively. 

It is apparent that considering NP to effect the charged quark currents, especially the top 
quark interactions, will effect the mixing amplitude. On the computational level it will lead to new 
contributions on the full theory side resulting in a change of the Wilson coefficient at the weak 
scale. 

1.4.2 |Ai?| = 1 transitions 

In this section we consider another type of rare FCNC processes, namely the radiative decays of 
the neutral B mesons. On quark level the transitions that we will be interested in are 

6-^s7, b^sg, b^st^r, b ^ svi> , (1-54) 

giving rise to different decays on the hadronic level. Note that while we could again consider the 
light down-type quark to be either d or s we shall commit to the case of s, because the experimental 
sensitivities for \AB\ = 1 processes that we will be studying are better when the final state quark is 
the s quark. The results of the matching can however be applied to d case by simple sod change. 
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Very similarly to the FCNC decays of the top quark described in section 1.2.2, within the 
SM transitions given in Eq. (1.54) can not proceed through tree- level diagrams, but through two 
insertions of charged current interactions at one-loop level. The effective Lagrangian of the low 
energy effective theory, from which the heavy vector bosons and the top quark have been integrated 
out and is adequate for description of processes (1.54), can be written as 



+ 



AGf 



10 
j=3 



GiO, + G,pO 



where Aj stands for A^*'' defined in Eq. (1.46). The relevant operators read 



(1.55) 



08 



errii, 
(4^ 

Osinb 
(4^)2 



{sLa^^^T%R)Gl,, 



O9 



(47r)2 

e2 
(47r)2 

(4^)2 



(sl7^6l)(^V)> 
(sl7^&l)(^^7^(1-7')^) 



(1.56) 



Here are SU{3)c generators in fundamental representation and cr^i/ = i/2[7^,7^]. The electro- 
magnetic and gluonic field strength tensors are defined as 



9s fi 



abc 



(1.57) 



where fahc are the S?7(3)c structure constants. Since they are not that crucial for our analysis, we 
omit the definition of the remaining four-quark operators 03,...,6i which can be found for example 
in Ref. [67]. 

In the procedure of matching the SM to the effective theory (1.55), again virtual top quarks 
turn out to play the dominant role. Representative diagrams to be computed are given in Fig. 1.5. 
In the first diagram the photon and/or gluon are considered to be on-shell, contributing to C7 
and Cg respectively. The remaining two penguin diagrams contribute to Cg, Cio and C^y. Here 
the gauge bosons are off-shell coupling further to the lepton pair. We note that Cio coefficient 
can only receive contributions from Z mediated process due to the purely vectorial structure of 
the photon- lepton coupling. In addition to the penguin diagrams, Cg, Cio and C^y also receive 
contributions from box diagrams. 

Calculation of the box diagrams proceeds in a similar fashion as described in section 1.4.1 for 
the mixing case. To compute all the Wilson coefficients coming from the penguin diagrams, we 
have to consider the photon and gluon to be in general off-shell and perform an expansion of the 
amplitudes to the second order in external momenta, neglecting throughout rrig and rn^ terms. For 
the off-shell massive Z boson the computation is somewhat simplified since its momentum can be 
neglected compared to its mass. 
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Figure 1.5: Representative penguin and box Feynman diagrams to be calculated in the matching 
of SM at LO in QCD to the effective Lagrangian given in Eq. (1.55). Quarks running in the loops 
are up-type and i denotes their flavor. In addition to the presented penguin diagrams one needs 
to consider also diagrams where the gauge bosons are emitted from the W boson (except for the 
gluon) or the external quark legs. Diagrams with would-be Goldstone bosons which are present 
when working in general gauge for the weak interactions are not presented. 



The amplitudes for processes can be written as 



AfMb ^ s^+r )box = iXi^j^^Boix,)[sLl''bL][iLl^.iL] , (1.58a) 



V2 (47r)2 s2 
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Bo{x,)[sLi''bL][i>Ll^iyL], (1.58b) 



Aiuiiib — > sZ) ~ \X 



AGf 



2 cw 



Ain\\{b S7) = lA 



4Gf 



Ahi\\{b sg) = lA, 



V2 (4^)2 
AGf 9. 



T^m^ — Co(x,)[sL7"6L]^o(fc), (l-58c) 
Do{xi)[sL{k'-i'' - fc"|^)foL] + mf,Z?[,(x,)[siia"^fc^6fl]l A^{k) , (1.58d) 



£^o(^0[sL(fc'7" - k"l^)T%L] + mbE'o{x,)[sLi<j"^kpT%u] Gl{k) . (1.58e) 



^/2 (4^)2 

The first two amplitudes are related to the box diagrams, while the rest are obtained from the 
penguin diagrams. We have introduced the abbreviations cw and .sw which denote the cosine and 
sine of the Weinberg mixing angle. Further, k denotes the momentum of the gauge boson, which 
can be neglected in Eq. (1.58c) as argued above. 

The loop functions Bq, . . . ,Eq depend on Xi = mf/m^^, where rm is the mass of the up- type 
quark running in the loops of the diagrams. The subscript stands to remind that the process 
was calculated at LO in QCD. The analytical expressions for the loop functions are given in the 
Appendix B.3. We note, that again due to the CKM unitarity relation (1.47) all the Xi independent 
terms encountered in the computation are dropped. 

Functions Bq, Cq and Dq remain dependent on the choice of the gauge for weak interactions. 
The gauge independence is recovered once all the contributions to a particular physical final state 
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are summed up. In particular, diagrams with off-shell photons and Z bosons contribute to the 
leptonic final states as the bosonic field gets contracted into a propagator and further coupled 
to the leptons. Combining these contributions with those coming from box diagrams gives the 
following gauge independent combinations 



^fuii(& ^ sl+l ) = iA^ 



4G/ 



V2 (47r)2 



4=2 



+ Co{x,) + Do{x^)) [sLj^bL] [Ij^i] (1.59a) 



^full(^ svv) = iX 



/2 (47r) 



2BoiXi) + Coix^) r u, 1 r- /I \ 1 



4s 



(1.59b) 



w 



To complete the matching procedure we have to compare the full theory amplitudes, with the 
amplitudes computed within the effective theory (1..55). For the on-shell photon and gluon the 
matching is straightforward since only contributions on effective theory side are tree-level insertions 
of operators. Same holds for the neutrino final state and the axially coupled charged leptons 
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(1.60) 
(1.61) 
(1.62) 
(1.63) 



Since the matching was performed at LO in QCD, the expressions again apply to the high /i^i/ scale. 
For the extraction of Cg coefficient we have to, in addition to the trivial tree- level contribution of 
Og, take into account the one-loop contribution of O2 operator presented in Fig. 1.6. The amplitude 
on the effective theory side is 



.4Gf 



\/2 (47r) 



AtCg(/ivK) 



(1.64) 



w) 



+ log + log Xi + 1 



The 2/e UV divergence is removed using the MS renormalization'\ The constant term is charac- 
teristic to the naive dimensional regularization, which has been employed in the computation [52]. 
Because C2{fJ-w) = 1) we can write down, comparing (1.59a) and (1-64) the final expression for Cg 

^This means that O2 and Og mix under QED renormalization. 
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Figure 1.6: The one-loop contri- 
bution of O2 operator to the h — )■ 
sl'^l~ on the effective theory side. 



As in the case of | Ai?| = 2 process, all of the Wilson coefficients 
obtain higher order perturbative QCD corrections [68, 69, 70]. 

Consideration thereof introduces renormalization mixing among many of the operators involved in 
I AS I = 1 processes. Application of RG methods allows us to establish the QCD running of the 
Wilson coefficients [71, 72, 73, 74] which can then be, re-summing the large logarithms, run down 
to appropriately low scales, where the operator matrix elements can be evaluated thus making it 
possible to compare theoretical predictions with the experimental measurements of decay rates. 



1.5 The main strategy 

Having introduced the two phenomena of top quark physics that are interesting for NP searches 
and presenting also the importance of top quark contributions in B physics, we devote the last 
section of the introductory chapter to introduce our main strategy of specifying and analyzing the 
deviations from SM physics in the top quark sector. We closely follow the reasoning presented in 
Ref. [41], where the case of NP generating FCNC top quark decays was considered. We apply the 
same approach to the case of charged currents as well^. 

The main idea is illustrated in Fig. 1.7 and it starts with the assumption that there exists some 
NP at the energy scale A, which is much higher than the electroweak energy scale. As this NP is 
integrated out, it generates operators at the electroweak scale (denoted ^t)-, which consist of SM 
fields only, and are invariant under the SM gauge group (1.1). Making use of the highly stressed 
property of OPE (1.29), that the effects of higher dimensional operators generated at high scale A 
come suppressed with higher powers of 1/A, we avoid committing to a particular UV completion 
of the SM, rather working in the framework of an effective theory, described by the Lagrangian 

CoK = CsM + ^Y.^'^' + ^-^^ + O(VA^) , (1.66) 

i 

where Cqm is the SM part, and Qi are dimension-six operators with the aforementioned properties. 
The set of operators forming Ces are chosen such that they generate tqZ, tqj, tqg vertices, where 

*We note that this concept is not hnked only to NP in top quark physics which we are considering here. The 
indirect constraints on physics including new heavy degrees of freedom are often important and need to be considered, 
see for example [75]. 
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Integrating out t,W, ... Integrating out NP 




Figure 1.7: Illustration of the effective theory approach employed in our analysis. First step 
presents integrating out NP particles with mass well above the electroweak scale. The second step 
is integrating SM degrees of freedom above the scale of b quark mass to be taken when analyzing 
effects in B physics. 



q = u,c or modify the SM tWb vertex, depending on whether we are interested in NP in FCNC top 
quark decays (chapter 2) or the main decay channel (chapter 3). By parametrizing the appropriate 
vertices in most general way and analyzing the consequences of such modifications on top quark 
observables, we can quantify the effects of the effective operators in vCgfr on the top quark physics 
side. 

If however, we want to establish the effects of the same operators on the B physics side as 
well, we need to further match our effective theory (1.66) to the low energy effective Lagrangians 
responsible for |Ai?| = 2 and \AB\ = 1 processes given in Eq. (1.42) and Eq. (1-55) respectively, as 
is illustrated by the second second arrow in Fig. 1.7. By doing that, we gain the access to a variety 
of observables in B physics and come across an interesting interplay of top and bottom physics. 

We should note that, since the weak scale matching of NP contributions will be done at LO 
in QCD, there is an ambiguity of the order as(mj)/47r and a residual scheme dependence, when 
performing the RGE evolution at next-to-leading log, which we shall be employing. However as 
corrections to the matching are in general model dependent and thus beyond the scope of our 
effective theory approach (c.f. [7fi] for a more extensive discussion on this point). 



22 



Chapter 2 

NP in Top Decays: Neutral currents 



2.1 Introduction 

As we have pointed out in section 1.2.2 of the introductory chapter, SM predicts highly suppressed 
flavor changing neutral current processes of the top quark 

t^qV, V=Z,-f,g, q = c,u, 

while NP beyond the SM in many cases hfts this suppression. 

For the case of t — ?• qZ, 7 FCNC top quark decays, the eff'ective theory approach that we have 
described in section 1.5 has been used in Ref. [41], where the authors considered the constraints 
on operators that generate FCNC top quark decays from B physics observables. They found that 
contributions of some dimension six SU{2)l gauge invariant operators are not yet constrained by B 
physics data to such extent that the consequential top quark FCNC decays could not be observed 
at LHC, if the predicted sensitivities summarized in Eq. (1.27) are reached. On the other hand, 
gluonic operators governing the t ^ qg decays are not constrained by such indirect considerations 
and can, at NLO in QCD, contribute to t — )• qZ^ 7. In this chapter we therefore focus our attention 
on the top quark physics side only, analyzing the decay rates of FCNC top quark decays mediated by 
effective operators generating most general FCNC effective vertices. The correspondence between 
the Wilson coefficients of our operators and those of the SU{2) l invariant operators used in Ref. [41] 
is given in Appendix D.2. 

In the first part of the chapter, which is based on our published work [77, 78], we analyze 
the two-body decays t — )• qZ,^ at NLO in QCD. In Ref. [79] it was found that t ^ qg decay 
receives almost 20% enhancement from NLO QCD contributions while corrections to the t — )• qZ, 7 
branching ratios are much smaller. However, the authors of [79] only considered a subset of all 
possible FCNC operators mediating t — )• qV decays at leading order and furthermore neglected the 
mixing of the operators induced by QCD corrections. In the case of t — )• 57 decay in particular 
the QCD corrections generate a nontrivial photon spectrum and the correct process under study is 
actually t — t- qg'j. Experimental signal selection for this mode is usually based on kinematical cuts. 
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significantly affecting the spectrum. The vahdity of theoretical estimates based on the completely 
inclusive total rate should thus be reexamined. Finally, renormalization effects induced by the 
running of the operators from the NP scale A to the top quark scale are potentially much larger 
than the finite matrix element corrections. Although these effects are not needed when bounding 
individual effective FCNC couplings from individual null measurements, they become instrumental 
for interpreting a possible positive signal and relating the effective description to concrete NP 
models. 

The second part of the chapter is devoted to the study of t — )• q£~^£~ decays and is based on 
our published work [80]. The FCNC vertex is governed by the same effective Lagrangian as in the 
first part, but the neutral gauge boson is further coupled through SM interaction with the pair 
of charged leptons. The basic goal of this analysis is the identification of possible discriminating 
effects of different NP models in top FCNCs, by considering different types of observables, which 
become attainable due to the larger phase space of the final state. 



2.2 Framework 

In writing the effective Lagrangian that will generate the tZq, t^q and tgq vertices of the most 
general form, we rely on the notation of Ref. [37, 81]. Hermitian conjugate and chirality flipped 
operators are implicitly contained in the Lagrangian and contributing to the relevant decay modes 



eflf 



A2 



<ot + 



A2 



LR^LR 



'LR^LR 



+ {L ^ R) + h.c. . 



To explain the notation, operators considered are 



O 



L,R 



azZn 



— eF 



qL,Rl''tL,R 
qL,RCJ^''tR^L 



o 



LR,RL 



9zZ. 



qL,R(r^''tR^L 



(2.1) 



(2.2) 



0%RL = 9sG% qL,Ra>^''TatR,L 



and gz = 2e/ sin 20^. In addition to Ff^i, and G^^ that we have defined in Eq. (L57), we have 
introduced the derivative part of the Z boson field strength tensor 



Z 



duZa + 



(2.3) 



neglecting the terms with more than one vector field, since they are not relevant for our analysis. 
Finally v = 246 GeV is the electroweak condensate and A is the effective scale of NP. In the 
remainder of the chapter, since there is no mixing between chirality flipped operators we shorten 
the notation, setting a and b to stand for either a^, biR or ajj, bRL. The Feynman rules for t — )■ qV 
vertices generated by operators (2.2) are given in the Appendix A.L 

Note that in principle, additional, four-fermion operators might be induced at the high scale 
which will also give contributions to t — )• qV processes, however these are necessarily as suppressed. 
On the other hand, such contributions can be more directly constrained via e.g. single top pro- 
duction measurements and we neglect their effects in the present study. Throughout this chapter 
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we will be neglecting the mass of the final state (c, u) quark. Furthermore, when considering NLO 
QCD corrections we will regulate UV as well as IR divergences by working in d = 4 + e dimensions. 
This kind of approach necessitates performing phase space integration in d dimensions. We sug- 
gest the reader interested in dimensional regularization of IR divergences to consult the following 
references [82, 83, 84]. 

2.3 Two-body t — > qV decays 

In this section we present the results for the NLO QCD corrections to the complete set of FCNC 
operators which mediate t — )• qZ,^ decays already at the leading order (2.2). 

We start of by considering the virtual one-loop correction and the renormalization of UV diver- 
gences for both decay channels. We present the RGE effects linking the values of Wilson coefficients 
at the top quark scale to those at higher NP scale. Contributions from gluonic dipole operators 
are also taken into account. Next we turn our attention to the finite part of virtual corrections - 
the matrix element corrections, as well as the corresponding bremsstrahlung rates. For the t — )• ^7 
channel we also study the relevance of kinematical cuts on the photon energy and the angle between 
the photon and the jet stemming from the final state quark. We present our results in analytical 
form and also give numerical values to estimate the significance of NLO contributions. 

2.3.1 Operator renormalization and RGE 

We assume the effective a, b couplings are defined near the top quark mass scale at which we evaluate 
virtual matrix element corrections and a^. A translation to a higher scale matching is governed 
by anomalous dimensions of the effective operators and can be performed consistently using RGE 
methods. To employ this mechanism we need to examine the UV divergencies generated by the 
NLO QCD corrections. 




Figure 2.1: Feynman diagrams for one- loop virtual corrections to t — )• qZ,j decays. Squares mark 
the insertion of effective operators given in Eq. (2.2) and crosses the additional points from which 
Z or 7 can be emitted. 



In addition to the diagrams presented in Fig. 2.1, the first two diagrams should be accompanied 
by one particle reducible diagrams with SM gluon corrections attached to the external legs. These 
diagrams are taken into account by setting the quark fields forming the operators and 
to be renormalized q — ?■ y^~Z^q, t — )• y/Ztt. Since the final state light quark q is considered to be 
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massless, the corresponding field renormalization differs from tfie one of tlie initial top quark. Using 
the on-shell renormalization conditions [851 we obtain 
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where /i is the renormalization scale parameter, Cp = 4/3 is the color factor and separate track 
is kept of UV and IR divergences. Finally F denotes the Euler's gamma function. Including the 
tree- level LO diagrams, the quark field renormalization and the diagrams presented in Fig. 2.1 
yields the following amplitudes 
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where the complete expressions of form factors Fy are given in the Appendix B.l. We were able 
to crosscheck our expressions with those found in the literature. Namely, Eqs. (B.3-B.6) agree 
with the corresponding expressions given in Ref. [SO] for the B — )• Xs£'^i~ decay mediated by a 
virtual photon after taking into account that the dipole operator in [86] includes a mass parameter 
which necessitates additional mass renormalization. On the other hand, the two form factors 
for the photon case (B.l, B.2) are obtained from the corresponding Z form factors in the limit 
where the mass of the Z boson is sent to zero. To some extent we were also able to crosscheck 
the gluon operator induced form factors F^g and F^^ given in Eqs. (B.7, B.8). Namely, we find 
numerical agreement of the form factor's vector component with the corresponding expressions 
given in Ref. [87]. The crosscheck is only possible in the vector part, since the SM photon coupling 
appearing in [87] has no axial component. 

We note that F^ , F^^, F^ and F^^ contain UV divergences that necessitate additional operator 
renormalization which we carry out in the MS scheme obtaining the following renormalization 
factors 
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where Q = 2/3 is the electric charge of the up- type quarks and v is defined in Eq. (B.9) of the 
Appendix. The RG running is governed by the anomalous dimensions of the operators. Since 
operators do not have anomalous dimensions we assemble the remaining six operators into 

two vectors 



0, = {0],0f -^^^ 



. Of? 



RL, LR , 



(2.11) 



which do not mix with each other under QCD renormalization. The corresponding one-loop anoma- 
lous dimension matrix is the same for both chiralities and reads 
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2^ 



Cf 


8Cf/3 





Cf 



Cf{S 






5Cf - 2Ca 



(2.12) 



where Ca_ = 3. We note that to compute the last entry on the diagonal of the matrix (2.12) we need 
to consider virtual corrections to t ^ qg process mediated by Of operator. We have performed 
this calculation to crosscheck it with the well known result found in the literature (see for example 
Ref. [52]), we however refrain from explicitly showing the details of the calculation. 

Depending on the nature of new physics which generates the dipole operators at the scale A, 
the relevant LR operators might explicitly include a factor of the top mass. By redefinition of 
operators 

Olr = {mt/v)OLR , 

their running is altered by the additional mass renormalization Zm (found for example in Ref. [52]), 
which can be taken into account by adding QCf to the diagonal entries of ^lr given in Eq. (2.12). As 
we shall demonstrate, this effect is numerically not important for the interesting range of couplings 
and scales, which can be probed at the Tevatron and the LHC. 

We are interested in particular in the mixing of the gluonic dipole contribution into the photonic 
and Z dipole operators. For the case with no explicit top mass effect, LR and RL operators receive 
identical corrections and the effective couplings at the top mass scale read 



16 



o2 



(2.13a) 
(2.13b) 



where fit is the top mass scale, r] = as{A) / asi^it)-, i^i = 4/3/3o, ^2 = 2/3/3o and /3o is part of the one- 
loop QCD beta function (found for example in Ref. [52]). Assuming that no new colored degrees 
of freedom appear below the UV matching scale which would modify the QCD beta function, it 
evaluates to /3o = 7 above the top mass scale. If we include the top mass running in the RGE of 
LR operators, then ki^2 are modified to ki = 16/3/3o, K2 = 14/3/3o- 

We illustrate the effect of the RGE running in Fig. 2.2 where we plot 



when 67'^ (A) = 0, 



(2.14) 
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Figure 2. 2: The ratio of \b]'^ {fit)/bf{A)\ function of A, when bj' (A) = and the top mass 
scale ~ 200 GeV. Sohd hnes represent the case with no exphcitjop mass effect, while the dashed 
line corresponds to the Wilson coefficient of O^i? operator. The is not shown as its deviation 
from 5£p is unnoticeable on the plot. 



This shows how much b]'^ {/j-t) can be generated at the top mass scale fit — 200 GeV, due to the 
QCD mixing of the operators and the presence of the gluonic dipole operator at UV scale A. We 
see that for NP matching scales above 2 TeV the induced contributions to bj are around 10% of the 
bf in the UV. On the other hand, due to cancelations in the RGEs for the bf , these receive much 
smaller corrections (below 1 % for the interesting range). Including the top mass renormalization 
reduces the induced corrections to the coupling, however for UV scales of a couple of TeV or 
below, this effect is negligible. 



2.3.2 Matrix element corrections 



To consistently describe rare top decays at NLO in as one has to take into account finite QCD 
loop corrections to the matrix elements {qj\ Oi \t) and {qZ\ Oi \t) evaluated at the top mass scale as 
well as single gluon bremsstrahlung corrections, which cancel the associated infrared and collinear 
divergencies in the decay rates. The total FCNC top quark decay width to Z boson or a photon 
governed by the effective Lagrangian given in Eq. (2.1) is therefore, at NLO in QCD, the sum of 
r(t — )■ g.Z', 7) and V{t — )■ qgZ,j) decay rates, where the two-body final state decay width includes 
the virtual QCD corrections. 

Contributions due to the Oj^^^^j^ have already been computed in Ref. [79]. Here we expend 
the analysis to the operator basis given in Eq. (2.2), including results for Ofj^ current operators 
as well as for the admixture of the gluonic dipole operators Of^n . The final results that we are 
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after can therefore be parametrized in the fohowing way 
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where V = Z,^ and a''' = 0. Note that ^ appearing in the second and third row of Eq. (2.15) 

correspond to contributions from the gluonic operator and are therefore absent in the LO result, 
emerging only at order. 



Tree level expressions 

At the tree- level we only have F^, T^''^ and contributions, which we write in 4 + e dimensions 
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where = m^/mj. 
Virtual corrections 

The one-loop virtual QCD corrections to the decay amplitudes have already been presented in 
section 2.3.1, where the UV divergences were renormalized. This leaves us with UV finite form 
factors appearing in Eqs. (2.5, 2.6), which however remain IR divergent and the divergences are 
carried over to the expressions for t — t- qV NLO decay widths, for which the complete expressions 
are given in the Appendix C.1.1. Here we only outline their form 
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and stress that F^^"'^]^ all posses IR divergences, while F^ii^„ and f ^iil are finite. 
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Bremsstrahlung Contributions 

The relevant Feynman diagrams contributing to t — )• qgZ, 7 bremsstrahlung processes are given in 
Fig. 2.3. At the level of the decay width these diagrams give contributions of the same order in 
as as the one-loop virtual corrections presented above. Soft and collinear IR divergences emerge 
in the phase space integration and have to cancel the divergences present in r^''^"'* once we sum 
the two contributions. Computation of T{t — >• qgZ,j) from the first two diagrams of Fig. 2.3 



Of 



Figure 2.3: Feynman diagrams for t — )• qgZ, 7 bremsstrahlung process. Squares mark the insertion 
of effective operators given in Eq. (2.2) and crosses the additional points from which the gluon (in 
the first two diagrams) or Z, 7 (in the last diagram) can be emitted. 



gives contributions presented in Eqs. (C.11-C.15) that indeed render the sum of three-body and 
two-body final state decay width IR finite 
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We were able to crosscheck our results given in Eqs. (2.18) with the corresponding calculation done 
for a virtual photon contributing to the B — t- Xsi^i~ spectrum [88]. After taking into account 
the different dipole operator renormalization condition in [88] (including mass renormalization) 
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we find complete agreement with their results. was also cross-checked with the corresponding 
calculation of the t — >• Wb decay width at NLO in QCD [89]. Finally, we have compared our 
expression with the results given by Zhang et al. in Ref. [79]. In the limit — )• our results agree 
with those given in [79] , but we find disagreement in the rz dependence. After our first publication 
of these results in [77], we were made aware of a new paper in preparation by the same authors, 
which has now been published [90] and therein a corrected result for is given that coincides with 
ours. 

The remaining bremsstrahlung contributions are induced by the gluonic dipole operator. What 
needs to be pointed out here is that while final result (2.15) for F^ is finite, F"^ remains IR divergent. 
The divergences appear in F^ (squared contribution of third diagram of Fig. 2.3) and are not 
canceled by any of the virtual corrections we have considered. To cancel them we would have to 
consider the decay width for t ^ qg governed by the gluonic operator and include the one-loop 
virtual QED corrections. The corresponding Feynman diagram is shown in Fig. 2.4. 

The t — qjg decay process involves three (one almost) 
massless particles in the final state. Virtual matrix element 
corrections contribute only at the soft gluon endpoint {Eg = 0) 
and result in non-vanishing b'^b^ interference contributions. 
They involve IR divergencies which are in term canceled by 
the real gluon emission contributions. These also produce non- 
vanishing IftS'p contributions, and create a non-trivial photon Figure 2.4: Feynman diagram for 
spectrum involving both soft and colhnear divergencies. The one-loop virtual QED correction 
later appear whenever a photon or a gluon is emitted collinear for t ^ qg decay governed by 
to the light quark jet. An analogous situation is encountered Of^j^j^^. 
in the B — )• decay measured at the i?-factories. However, 

there the photon energy in the B meson frame can be reconstructed and a hard cut (E^^^) on it 
removes the soft photon divergence. The cut also ensures that the B — )• Xsg process contributing at 
the end-point £J-y = is suppressed. On the other hand, in present calculations the collinear diver- 
gencies are simply regulated by a non-zero strange quark mass, resulting in a moderate \og{ms/mh) 
contributions to the rate. The situation at the Tevatron and the LHC is considerably different. 
The initial top quark boost is not known and the reconstruction of the decay is based on trigger- 
ing on isolated hard photons with a very loose cut on the photon energy (a typical value being 
E^ > 10 GeV in the lab frame [91]). Isolation criteria are usually specified in terms of a jet veto 
cone AR = \J Ar/^ + A(/»2 where A?7 is the difference in pseudorapidity and A0 the difference in 
azimuthal angle between the photon and nearest charged track. Typical values are Ai? > (0.2 — 0.4) 
[49]. 

Rather than including QED corrected t ^ qg rate, we render the F^ finite by modeling the 
non-trivial cuts in the top quark frame with a cut on the projection of the photon direction onto 
the direction of any of the two jets {5rj = 1 — • pj /E-fEj), where j = g,q labels the gluon and 
light quark jet respectively. The effects of the different cuts on the decay Dalitz plot are shown in 
the left graph of Fig. 2.5. Since at this order there are no photon collinear divergencies associated 
with the gluon jet, the 5rg cut around the gluon jet has a numerically negligible effect on the rate. 
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On the other hand the corresponding cut on the charm jet - photon separation does not completely 
remove the divergencies in the spectrum. However, they become integrable. The combined effect 
is that the contribution due to the gluonic dipole operator can be enhanced compared to the case 
of B X,7. 

We present the full analytical formulae for the t — )• q^g and t — )■ qZg decay rates including the 
effects of kinematical cuts for the former channel in Appendix C.1.2. 

Numerical analysis 
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Figure 2.5: Left: The t — t- q^g Dalitz plot. Contours of constant photon and gluon infrared and 
collinear divergent contributions are drawn in red (dot-dashed) and blue (dashed) lines respectively. 
The collinear divergencies appear at the horizontal and vertical boundaries of the phase-space, while 
the IR divergencies sit in the top and right corners. The cuts on the photon energy correspond 
to vertical lines, the cuts on the gluonic jet energy to horizontal lines. Full green lines correspond 
to cuts on the jet veto cone around the photon. Right: Relative size of corrections to the 
Br(t — )• (77) at representative ranges of 5rc = 5r and -E™*. Contours of constant correction values 
are plotted for 6^ = (gray, dotted), = (red) and = —IP (blue, dashed). 

In all the numerical analysis of this section we use the following values for the parameters 

mvi/ = 80.4 GeV , mt = 172.3 GeV, mz = 91.2GeV, (2.19) 

fj, = mt , a s (mt) = 0.107 , sin^ = 0.231 . 

Turning first to t — )• 97 decay we show in the right graph of Fig. 2.5 the induced correction to the 
tree- level Br(t — t- (77) for representative ranges of 5r and i?™*. We observe, that the contribution 
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of can be of the order of 10 — 15% of the total measured rate, depending on the relative sizes 
and phases of Of,'^ and on the particular experimental cuts employed. Consequently, a bound 
on Br(t — )• ^7) can, depending on the experimental cuts, probe both b^''^ couplings. In order to 
illustrate our point, we plot the ratio of radiative rates T{t — )• qj)/T{t — )• qg), both computed at 
NLO in QCD versus the ratio of the relevant effective FCNC dipole couplings \b'^ /b^l in Fig. 2.6. 
The NLO T{t — t- qg) result is taken from Ref. [79]. We show the correlation for two representative 




0.0 0.2 0.4 0.6 0.8 



Figure 2.6: The ratio of radiative rates T{t — )• q^)/T{t — )• qg) versus the absolute ratio of the 
relevant effective FCNC couplings \b'^ Two representative choices of experimental kinematic 
cuts are shown. The two regions represent the possible spread due to the unknown relative phase 
between and W couplings, while the lines correspond to maximal positive (full) and negative 
(dashed) interference V^b^ . See text for details. 



choices of experimental kinematic cuts for the t ^ q^y decay. The vertical spread of the bands is 
due to the variation of the relative phase between 6^ and b^ couplings. We also display the two 
interesting limits where the Wb^ interference is maximal positive (zero relative phase) and negative 
(relative phase vr). We see that apart from the narrow region around \b'^ /b^\ ~ 0.2, where the two 
contributions may be fine-tuned and conspire to diminish the total t ^ q^y rate, the two radiative 
rates are well correlated. In particular, depending on the kinematical cuts employed, there is a 
natural lower bound on ratio of decay rates, valid outside of the fine-tuned region. Finally, for 
\b'^' /b^\ > 0.6 the correlation becomes practically insensitive to the particular experimental cuts 
employed and also the unknown relative phase between b"' and b^ couplings. 

For the t — t- qZ decay channel we present some numerical values to estimate the significance of 
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QCD corrections. In particular we parametrize the decay width given in Eq. (2.15) as 



167r^^lA4' ' 



47r 



CpVa 



H \uZ\2 



A4 



2^6 + -rCpyb 

47r 
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A4 



in 



v^mt 



2Re{a''*bS}-^CFyag - 2lin{a''*bsy^CFyag 
47r 47r 



A4 



2Re{b^*b^}^CFybg - 2lm{b''*ba}^CFybg 
47r 47r 



(2.20) 



Here Xj stand for the tree-level contributions, while yi , jji denote the corresponding QCD correc- 
tions. Numerical values of the coefficients are given in Tab. 2.1. We see that corrections due to the 
gluon dipole operator are an order of magnitude smaller (except Ug, which is even more suppressed) 
than corrections to the Z operators themselves and have opposite sign. 



Xb = 


2.36 


— 


1.44 


^ab — 


1.55 






Vh = 


-17.90 


Va = 


-10.68 


Vab = 


-10.52 


Va -- 


= 0.0103 


Vhg = 


= 3.41 


Hag = 


-- 2.80 


Vbg = 


2.29 


Uag 


= 1.50 



Table 2.1: Numerical values of coefficient functions appearing in Eq. (2.20). 



Next we investigate the relative change of the decay rates and branching ratios when going 
from LO to NLO in QCD. The results are presented in Tab. 2.2. We see that the change in the 





=b3 = 


=b9 = 


= b^,b3 = 


= 0, = b9 


=0,5^= ¥ 


pNLO y-pLO 


0.92 


0.91 


0.92 


0.95 


0.94 


BrNLO/BrLO 


1.001 


0.999 


1.003 


1.032 


1.022 



Table 2.2: Numerical values of r^LO/pLO ^^^^ Br^L^(t qZ)/Br^'^{t qZ) for certain values 
and relations between Wilson coefficients. 



decay width is of the order 10%. There is a severe cancellation between the QCD corrections to 
r(t — )■ cZ) and the main decay channel r(t — )• bW). This cancellation causes the change of the 
branching ratio to be only at the per-mille level when ¥ is set to zero. In the case when only 
operators Ofj^ are considered this cancellation is anticipated since the NLO correction to is 
of the same form as the correction to the rate of the main decay channel. If we treat b quarks as 
massless, exact cancellation is avoided only due to the difference in the masses of Z and W bosons. 
It is more surprising that similar cancellation is obtained also when only the dipole Z operator is 
considered. However, setting ¥ = or ¥ = b^ , the impact of QCD corrections is increased by an 
order of magnitude and reaches a few percent. 
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2.3.3 Summary 

To summarize, we have presented a study of t — )• qZ, 7 decays mediated by the effective operators 
given in Eq. (2.2) at NLO in QCD. We found that QCD corrections can induce sizable mixing of 
the relevant operators, both through their renormalization scale running as well as in the form of 
finite matrix element corrections. These effects are found to be relatively small for the t — )• qZ 
decay, but can be of the order 10% in the t ^ qj channel, depending on the kinematical cuts 
employed. The accurate interpretation of experimental bounds on radiative top processes in terms 
of effective FCNC operators requires the knowledge of the experimental cuts involved and can be 
used to probe Of, r contributions indirectly. 
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2.4 Three-body t — > qi^i decays 

This section is devoted to the study of t — )• qi^i~' with the basic goal of identifying discriminating 
effects of different NP models in top FCNCs which can be approached by the experimental study. 
Exploring the three-body decay channel brings about two main advantages. First one is the larger 
phase-space which offers more observables to be considered - in particular the angular asymmetries 
among the final state lepton and jet directions. The second advantage is that the three-body final 
state that we are to consider is common to both Z and 7 mediated FCNC decays. Since some 
BSM models predict observable FCNC top quark decays in both Z and 7 channels the interference 
effects in the common three-body channel is something worth exploring. 

Since the standard forward-backward asymmetry for the leptons vanishes in the photon medi- 
ated decays, we consider another asymmetry which we call the left-right asymmetry and is associ- 
ated with the lepton angular distribution in the lepton-quark rest frame (see section 2.4.2). This 
asymmetry is nontrivial also in the photon mediated decays. We explore the ranges of values for 
these two asymmetries in t — )• qi~^£~ decays mediated by both Z boson and the photon. We also 
consider the interference effects as we expect them to significantly affect the ranges of the asym- 
metries. Our results can serve as a starting point for more elaborate investigations of experimental 
sensitivity to the proposed observables including QCD corrections, proper jet fragmentation and 
showering and the impact of experimental cuts and detector effects 

2.4.1 Effective Lagrangian 

The effective Lagrangian governing the FCNC top quark decay is 
considered to be the same as given in Eq. (2.1). We shall assume 
the FCNC mediating gauge boson (Z boson or the photon) to fur- 
ther couple with a pair of charged leptons. Note that the gluonic 
operator will not play a role in the analysis of this section since the 
gluon does not couple to leptons. For completeness we present here 
the part of SM Lagrangian that couples charged leptons with the 
photon and the Z boson 

Ci = gzZ^ [cR iR-f^'iR + CL IlI^l] + eA^l^l , (2.21) 

where the Z couplings are cr = siv? 9w and cl = —{cos29w/2). 
The Feynman diagram of the process that we shall be studying is given in Fig. 2.7. Since top quark 
is massive enough for the produced Z boson to be on-shell we shall use the standard Breit-Wigner 
formula to accommodate for its finite decay width by replacing (s — m?z) ~^ ~ + '^'^z^z) 
in the denominator of the propagator, where s denotes the square of the intermediate Z boson 
momentum and Tz is its total decay width. 



t q 




Figure 2.7: t ql+i' 
Feynman diagram, where 
the FCNC operators 

^f,R' ^LR,RL^ are given 
in Eq. (2.2). 
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2.4.2 Observables 

We consider scenarios where detection of a NP signal in the FCNC decay channel t — )• qi~^£~ 
could be most easily complemented by other observables in the same decay mode. This would 
allow distinguishing between different possible effective amplitude contributions and thus different 
underlying NP models. We neglect kinematical effects of lepton masses and the light quark jet 
invariant mass, as these are expected to yield immeasurably small effects in the kinematical phase- 
space set by the large top quark mass. 

We start with the double-differential decay rate dT/(duds), where s = ?tT'£+£_ is the invariant 
mass of the lepton pair and u = is the invariant mass of the final state quark (jet) and the 

lepton of positive charge . Integrating this decay rate over one of the kinematical variables, we 
obtain the partially integrated decay rate distributions (dP/du, dP/ds), while the full decay rate 
(P) is obtained after completely integrating these distributions. The branching ratio is obtained by 
normalizing the decay width to the width of the main decay channel (l.f3). 

The differential decay rate distribution can also be decomposed in terms of two independent 
angles, as defined in Fig. 2.8. In the rest-frame Zj = cos 9j measures the relative direction 




Figure 2.8: Definition of two angles relevant to our analysis. The arrows represent the three 
momenta of particles. Left diagram corresponds to the lepton pair rest-frame and the right diagram 
to the rest-frame of the positively charged lepton and the light quark. 



between the negatively charged lepton and the light quark jet. Conversely, in the rest-frame of 
the positively charged lepton and the quark jet, we can define zi = cos 6i to measure the relative 
directions between the two leptons. In terms of these variables, we can define two asymmetries 
{i = j, £) as 

Ai = ~ , (2.22) 

where we have denoted Pz^go as the integrated decay rates with an upper or lower cut on one of the 
Zi variables. We can then identify Aj = Ap-Q as the commonly known forward-backward asymmetry 
(FBA) and in addition define Ai = ^lr as the left-right asymmetry (LRA). The two angles and 
the asymmetries they define are related via a simple permutation of final state momentum labels 
between the quark jet and the positively charged lepton, and consequently via a u o s interchange. 
Since the asymmetries as defined in Eq. (2.22) are normalized to the decay rate, they represent 
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independent observables with no spurious correlations to the branching ratio. On the other hand, 
correlations among the two asymmetries are of course present and indicative of the particular NP 
operator structures contributing to the decay. 



2.4.3 Signatures 

Next we study the signatures of various possible contributions to the t — )• q£~^i~ decay using the 
integrated FBA and LRA observables defined in the previous section. Before exploring individual 
mediation cases a general remark is in order. Since all the effective operators of our basis (2.2) come 
suppressed with an undetermined NP cut-off scale, the actual values of the effective couplings (a, 
b) are unphysical (can always be shifted with a different choice of the cutoff scale). The total decay 
rate determines the overall magnitude of the physical product of the couplings with the cut-off 
scale. On the other hand relative sizes or ratios of couplings (independent of the cut-off) determine 
the magnitude of the asymmetries. The extremal cases are then naturally represented when certain 
(combinations of) couplings are set to zero - often the case in concrete NP model implementations. 



Photon mediation 



As pointed out in section 2.3.2, the direct detection of energetic photons is considered to be the 
prime strategy in the search for photon mediated FCNC top quark decays. However the t — )■ qi~^i~ 
channel, where the photon is coupled to the charged lepton pair can serve as an additional handle. 
Due to the infrared pole in the di- lepton invariant mass distribution we introduce a low s = m'j /m^ 
cut denoted Smin = and present the total decay width as its function 
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(2.24) 
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LR\ 



+ 



(2.23) 



The physical cut is of course at e = 4m|. 
We have also define an auxiliary variable 
summarizing the relevant NP parameter depen- 
dencies. Function depends only on the di- 
lepton invariant mass cutoff Smin and is pre- 
sented in Eq. (D.2) of the Appendix. 

While the FBA vanishes identically, due to 
the purely vectorial coupling of the photon to 
the leptons, the LRA can be written in the fol- 
lowing way 




^/7 [GeV] 

Figure 2.9: The dependence of the photon- 
mediated LRA on the low di-lepton invariant 
mass cut e. 



where function (7^ is presented in Eq. (D.3) of the Appendix and also depends only on the cut of 
the di-lepton invariant mass. Consequently LRA does not depend on the effective dipole couplings 
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in any way, however showing a non-trivial dependence on the low s cut, which we plot in Fig. 2.9. 
We see that the value of the integrated LRA is highly sensitive to the cut having value 1 in the 
limiting case of ^/e — )• 0, decreasing with the cut set higher and even exhibiting a change of sign 
between 60 and 80 GeV. 



Z mediation 

Current search strategies for t — )• qZ decays actually consider t — )• q£^i^ decay channel, where 
the charged lepton pair is to be identified as the decay product of the Z boson. This is achieved 
by imposing a cut on the invariant lepton mass around the Z mass to reduce backgrounds. As 
long as such cuts are loose compared to the width of the Z, we do not expect them to affect our 
observables. The decay width and the two asymmetries can be written as 



r 



4 4 



fAA + fsB + fcC , (2.25a) 



- ^"-'^fAA + fsB + fcC^ 

- fAA + fsB + fcC ' ('-''^^ 
where the parts depending on the NP generated FCNC top quark couplings have been separated 
from the parts that depend just on the SM parameters and the phase space integration. The NP 
dependent parameters introduced are 

_,2 ilZ |2 , |2 m^\h^ \^ - \h^ |2 

2 2 
mi Re{bfX* + ^lL4*} . ^ _ mi Re{bfX*-^gL4.*} r 

^ — 7, ^+ 1 7 — 7, ^~ ) 

V 2 V 2 

where L± = ^ sin^ 9w it | cos'^ 26yir are the factors coming from the charged lepton couplings to the 
Z boson governed by the Lagrangian given in Eq. (2.21). 

On the other hand the remaining parameters fi and gi do not depend on the effective FCNC 
couplings and are presented in the form of integrals in Eqs. (D.6, D.7) of the Appendix. 

Performing a random sweep over the values of FCNC couplings that give the same value of the 
total FCNC decay width (2.25a), we explore the possible ranges and correlations between the two 
asymmetries (2.25b, 2.25c) in Fig. 2.10. On the same plot we also project the limits, where only 
dipole or only current interactions of the Z contribute. In Ref. [41] strong indirect limits were re- 
ported on the left-handed FCNC couplings of the Z coming from low energy observables. Therefore 
we also superimpose the possible predictions for the two asymmetries when these couplings are set 
to zero. 

We observe that the LRA can be used to distinguish between dipole and current FCNC couplings 
of the Z, while the FBA can distinguish the chiralities of the couplings. 
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Figure 2.10: The correlation of FBA and LRA in Z mediated decay. The gray area (sohd border) 
represents decays with aU possible current and dipole Z FCNC couplings. The orange area (dotted 
border) corresponds to decays with a£ set to zero, while the solid and dashed lines represent decays 
with only current and only dipole couplings respectively. The following numerical values have been 
used mz,t = 91.2 , 171.2 GeV, Tz = 2.5 GeV and sin^ Ow = 0.231. 



Interference of photon and Z mediation 

Several NP models predict comparable decay rates for t — )■ qZ, 7. This may in turn lead to a 
situation, where an experimental search using a common final state may be more promising than 
dedicated searches in each channel separately. In addition, the asymmetries in t — )• ql~^£~ may shed 
additional light on the specific couplings involved. The decay rate in this case depends again on 
the di-lepton invariant mass cutoff 



where the pure photon contribution F'^ is given in Eq. (2.23) while the pure Z and interference 
contributions can be written as 




(2.27) 
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Here Wi,...,W4 are the newly introduced NP dependent constants containing both Z and 7 
effective FCNC couplings 



Wi 



mt 1 



Wo 



m} 1 



(2.29) 
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The two asymmetries can be expressed as the following fractions 

- 4/3 + 47) + fw(2{W2 - Wi) + Wa- W^) 
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The newly introduced functions fi and gi now depend on the Z boson parameters as well as the 
di-lepton invariant mass cutoff e. They are presented in Eqs. (D.8) of the Appendix. 
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Figure 2.11: The correlation of FBA and LRA in Z and 7 mediated decay. The gray area (solid 
border) represents decays with all possible Z and 7 FCNC couplings. The blue (dotted border) area 
corresponds to decays with only current Z FCNC couplings. For comparison, the orange (dashed 
border) area represents Z mediated decays. Numerical values used are the same as in Fig. 2.10 
with addition of \/e = 40 GeV. 



Performing a random sweep across both 7 and Z NP FCNC couplings that give the same decay 
width (2.27) we explore the possible correlation between the FBA and the LRA in this scenario. We 
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present the results with a fixed cut on y/e set to 40 GeV in Fig. 2.11. 
for the case when only the current FCNC Z couplings contribute, 
interference effects can produce a larger LRA compared to the case 

In Ref. [41] upper bounds on coefficients accompanying the op- 
erators responsible for FCNC t — )• cZ and t — )• 07 are presented. 
Using transcription formulae presented in Eqs. (D.IO) of the Ap- 
pendix, we can evaluate FBA and LRA associated with these upper 
bounds. The numerical values are presented in Tab. 2.3. 

These values serve just for illustration that nonzero values of 
asymmetries can indeed be obtained. They do not represent any 
kind of upper bounds for asymmetries. There is no reason to think 
that the highest allowed values of coefficients (D.IO) are to give 
the largest possible asymmetries which are complicated functions 
of these coefficients. 



We also present possible points 
We observe that in principle 
of pure Z mediation. 





only Z 


Z and 7 


^FB 


0.045 


0.035 




0.206 


0.226 



Table 2.3: Values of FBA and 
LRA for highest allowed coef- 
ficients given by Fox et al. in 
Ref. [41], = 40 GeV. 



2.4.4 Summary 

We have considered the top quark decay t — )• ql^l^ as a probe of BSM physics manifested in 
FCNC top quark transitions. In addition to the branching ratio, we have defined two angular 
asymmetries which can serve to further discriminate between different NP scenarios. Comparing 
possible contributions to the decay mode via SM 7 and Z mediation we can draw the following 
conclusions: large values of FBA (|^fb| 3> 0.1) cannot be accounted for in decay modes mediated 
by 7 bosons as long as we assume these bosons to have SM couplings to the charged leptons. 
We have shown in Fig. 2.10, that the ApB G [-0.12,0.12]. 

A measured point in (AfBj^lr) plane could exclude models with only current or only dipole 
FCNC couplings of Z if it were located off the solid and dashed black lines in Fig. 2.10. Treating 
the Z and photon as indistinguishable mediators expands the allowed LRA region to larger positive 
values. 

Current experimental sensitivity studies look at the two-body decay modes t — t- cZ, 7 [49]. Our 
analysis may be applicable to the potential measurement of t — )• cZ at ATLAS since they will be 
identifying the Z boson through its decay to a lepton pair. Angular asymmetries of this pair and 
the remaining hard jet could provide additional information on the tcZ FCNC vertex. The t ^ cy 
decay is generically characterized by a single high pT photon. Current search strategies for this 
FCNC include the detection of this photon, and not its eventual decay to a lepton pair. 

In order to fully explore our decay mode, one would need to relax or modify certain criteria 
used by current search strategies to reduce SM backgrounds. In addition, the reconstruction of 
the LRA might require top quark charge tagging. In principle our results are applicable also to 
the purely hadronic decay modes, where the two leptons are replaced by 6-tagged jets for example, 
however in this case the asymmetries are compromised by the lack of knowledge of the sign of the 
b quark charges. 
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NP in Top Decays: Charged currents 



3.1 Introduction 

In this chapter we analyze the possible deviations from the SM coupling strength and structure 
of the charged quark currents. Following the strategy outlined in section 1.5 of the introductory 
chapter we explore on one hand the implications of such deviations on the low energy B physics 
and on the other hand delve into the possible consequences to be observed in top quark decays, 
particularly focusing on the main decay channel and the tWb interaction. Results presented in this 
chapter are based on our published work [92, 93, 94]. 

Unlike in the case of FCNC top quark decays, the complete analysis of low energy constraints 
for model independent anomalous charged quark currents has not been performed yet. In particular 
analysis of anomalous tWd,s couplings on B meson mixing has been attempted in Ref. [95, 96], 
using however a subset of all possible effective tWd, s operators. Furthermore effects in 6 — )■ 57 
transition have been analyzed in Ref. [97]. Again, only a subset of the operators that we shall be 
considering has been used and no other AB = 1 transitions have been ruminated. 

We therefore first set out to upgrade the analysis of indirect implications available in the 
literature focusing on observables in AB = 2 and AB = 1 processes presented in sections 1.4.1 and 
1.4.2 respectively. The experimentally well measured low energy observables which mostly agree 
with the SM predictions provide constraints on the anomalous charged quark current structures. 
Since the top quark with its possible anomalous interactions is in these processes a virtual particle 
we label such constraints on top quark physics as indirect. 

Only then do we turn to the direct top quark physics where the structure of tWb vertex can 
be probed by analyzing the main decay channel of on-shell produced top quarks at Tevatron and 
LHC. The derived indirect constraints give us an idea of how much room there is left for deviation 
away from SM predictions given that the NP can be parametrized in the form of anomalous tWb 
couplings. 

We briefly comment that there is another class of observables in which virtual top quarks play an 
important role called electroweak precision observables (EWPO). Namely, anomalous tWb vertices 
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can effect the S, T, U oblique parameters and the Z ^ bb decays (see for example [98, 99]). Following 
our publications of the results that are presented in this work, the analysis of anomalous top quark 
couplings on EPWO was given in Refs. [100, 101]. Although the operator basis parametrizing NP 
used by the authors of Refs. [100, 101] does not coincide with the one used in this work, especially 
in regards to flavor structure, there are some common operators for which we can confront the 
obtained indirect bounds to find they are comparable and compatible. 

This chapter is structured as follows. The first section is devoted to the formulation of the 
effective theory. Specifying our framework in terms of effective operator basis we proceed to analyze 
the effect in B physics, giving detailed study of AB = 2 transitions, namely B meson mixing and 
AB = 1 transitions comprised of FCNC B meson decays. In the last section we focus on the top 
quark physics and the effects of anomalous tWb couplings on helicity fractions in its main decay 
channel. Here we perform our analysis at NLO in QCD, which seems sensible, since as pointed out 
in section 1.2.1, NLO QCD corrections prove crucial in consideration of helicity fractions within 
the SM. We show that the direct bounds are in a sense complementary to the indirect bounds and 
reveal a nice interplay between top and bottom physics. 

3.2 Effective Lagrangian 

Following the strategy outlined in section 1.5, our first objective is to specify the operator basis 
Qi of our interest, namely the determination of gauge and fiavor structure of the dimension six 
operators appearing in the effective Lagrangian (1.G6). To analyze the effects in B physics, we will 
have to further perform the second step illustrated in Fig. 1.7. Much like in the SM case presented 
in sections 1.4.1 and 1.4.2, we will have to integrate out the SM degrees of freedom with masses 
greater than that of the b quark, matching £eff to low energy effective Lagrangians (1.42, 1.55). By 
doing that we shall gain access to different B physics observables of interest and the possibility to 
see how predictions are affected by NP. 

3.2.1 Gauge structure 

Our operator basis consists of all dimension-six operators invariant under the SM gauge group that 
generate charged current quark interactions with the W. The possible gauge structures we can use 



Here Q stands for the quark SU{2)l doublet u, d are the up- and down-type quark SU{2)i singlets 
respectively, are the SU{2)l Pauli matrices. In addition we have the covariant derivative and 



are [102] 




(3.1a) 
(3.1b) 
(3.1c) 
(3.1d) 
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field strength definitions 

D, = d, + i^WX + -4^,Y, (3.2) 

and finally 4>u,d are the up- and down-type Higgs fields (in the SM (j)u = cj) = ir^(/)^). The operators 
in Eq. (3.1) are written with quark fields in the interaction basis being flavor universal. 



3.2.2 Flavor structure 

On the flavor side, we restrict the structure of the operators to be consistent with Minimal Flavor 
Violation (MFV) hypothesis [103, 104, 105], which postulates that even in the presence of NP 
operators, Yukawa matrices present in the SM remain the sole source of flavor violation. 

The way to implement this concept is to make the operators formally invariant under the SM 
flavor group (1.4) and insist that the only symmetry breaking spurionic fields in the theory 
are the up and down quark Yukawa matrices Yu^d, introduced in Eq. (1.3), formally transforming 
as (3, 3, 1) and (3, 1, 3) respectively. 

The four distinct quark bilinears appearing in Eq. (3.1) have different transformation proper- 
ties under Q^^: ud, QQ, Qu and Qd transforming as (1,3,3), (1®8, 1,1), (3,3,1) and (3,1,3) 
respectively. Prom them we can construct the most general invariant structures as 

uY^AudYdd, QAqqQ , QAquYuU, QAqdYdd, (3.3) 

where Axy are arbitrary polynomials of YuYu and/or Y^yJ, transforming as (1 ® 8, 1, 1). 

In order to identify the relevant flavor structures in terms of physical parameters, we can 
without the loss of generality consider Y^^^ condensate values in the down basis in which (1^) is 
diagonal (1.7) 

{Yd) ~ diag(0, 0, mb)/vd , {Y^) ~ 0diag(O, 0, mt)/v^ . (3.4) 

Here V is the SM CKM matrix (1.10) and we have introduced separate up- and down-type Higgs 
condensates Vu,d- We have also neglected the masses of first two generation quarks, which is the 
approximation we will be using throughout this chapter. 

Further, once we assume electroweak symmetry breaking (see Eq. (1.8) and text around it) we 
rewrite the quark fields in the mass eigenbasis. Making the fiavor indices and chirality explicit 
these fields are 

Qi = (ykiULk,dLi)^ , UiR, diR. (3.5) 

We consider first the simplest case of linear MFV where Axy is such that powers of Yd in (3.3) 
do not exceed 1, and powers of Y^ do not exceed 2. Obtained flavor structures are given in flrst two 
columns of Tab. 3.1. Following [106], the generalization of the above discussion to MFV scenarios 
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where large bottom Yukawa effects can be important is straight forward. We implement it by 
raising the highest allowed power of Yd appearing in (3.3) to 2. This gives us additional flavor 
structures presented in the last two columns of Tab. 3.1. We can see from the form of YuYu and 
Y^yJ with flavor indices explicitly written out 

2 2 

(YuY^h = "^yuVtj , {YdY^h = "i^si^sj , (3-6) 

that apart from these two forms, the only additional flavor structure that can be obtained is 

(YuY^YdYi),, = ^^V;,6s,. (3.7) 

This means that our list given in Tab. 3.1 is exhausting as forms of A more complex than presented 
therein generate no new flavor structures, rather just give higher powered overall factors. 





A = l 


A — V 


^ = YdXl 


A — YyYuYdY^ 


uY^AudYdd 




"^QrV^^tjQ, 


^Q3Q3 


'^'^QiV.tVttQs 


qaqqQ 


QiQi 


QAquYuU 


f^Q^VutR 




^^QsV^pR 




QAgdYdd 




-^^^^^tiOR 







Table 3.1: MFV consistent flavor structures of the quark bilinear parts of dimension-six operators 
that generate charged current quark interactions with the W . First two columns represent the 
simplest case of linear MFV with the highest allowed powers of Y^ and Yd set to 2 and 1 respectively 
while the last two columns represent MFV scenarios where large bottom Yukawa effects can be 
important, so the highest allowed Yd power is raised to 2. 

3.2.3 Final operator basis 

Before writing down the final set of operators that we will be analyzing we first put the obtained 
structures under further scrutiny and make some modifications. 

• Since QiQi is completely flavor universal, when coupled to the W it would modify the effective 
Fermi constant as extracted from charged quark currents compared to the muon lifetime. 
Existing tight constraints on such deviations [107] do not allow for significant effects in B 
meson or top quark phenomenology and we do not consider this structure in our analysis. 

• On the other hand, QiV^*VtjQj potentially leads to large tree-level FCNCs in the down quark 
sector when coupled to the Z. Explicitly, both the singlet and triplet operators (3.1b) include 
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the following term 



2cw 



Vt*Vt,[dL^rdLj\Z^^l, (3.. 



where (pQ is the lower SU{2)]^ component of (pd, having zero electric charge. After it acquires a 
VEV, such terms generate the aforementioned FCNC couplings. We therefore consider as our 
operator the linear combination of the singlet and the triplet, choosing the relative negative 
sign between them therby getting rid of the FCNCs in the bottom sector. 

• Similarly, QiV^*Vtbbji includes the following term 

- Vt* [dua'^'bR] MswFf^u + cwZf^u) , (3.9) 

which, once (po acquires a VEV, generates tree-level FCNC b — )• sj, Z transitions. These 
are already tightly constrained hy B ^ Xg^ and B — )• Xsi'^i~ [108]. Since all the charged 
current mediating SU{2) i invariant operators of dimension six or less containing such a flavor 
structure do necessarily involve either the Z or the photon, we drop this structure from our 
subsequent analysis. 

Taking these considerations into account, we obtain the set of seven effective dimension six operators 
invariant under the SM gauge group and consistent with MFV hypotheses that involve charged 
quark currents 

Qrr = Vtb[iBn^bR]{(t)iiD^4>d), (3.10a) 
Qll = [Q'3rVQ3](</'!,T'^i^M'Ad)-[Q37^Q3]Wi^M'/'rf)> (3-lOb) 
Q!ll = [QsrVQsK'/'i.r'^il^M'Ad) -[Q37^Q3](</':^i^M'/'rf)> (3-lOc) 

Q!'ll = [Qs^VQaK^l^r'^i^/.'/'d) -[Q37^Q3](0:^i^M'/'rf)> (3-lOd) 

Qlri = [Q'sT^a^^'tnjKW^,, (3.10e) 

Q'LRt = [Q3rV^''tR](/>„T^;,, (3.10f) 

QiRb = [QarV^'^^Rl^rfT^;,, (3.10g) 



where we have introduced 



Q'3 = Q^V,* = {tL,V;,d,Lf. (3.11) 



We note, that the final set of operators coincides with those considered in the B — t- XsJ analysis 
of anomalous tWb couplings [97] expanded by the three primed operators which originate from the 
structures given by the last two columns of Tab. 3.1, corresponding to higher order down- Yukawa 
MFV^. 



^In the final form of the operators we do not explicitly write out the rrit/vu and rrib/vd factors present in Tab. 3.1, 
technically shifting them to the Wilson coefficients. 
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Furthermore we do not make the operators hermitian, hence effects of operators Qj are accom- 
panied by C* and will be kept track of separately. 

Notice that starting with the most general MFV construction we are led to a set of operators, 
where largest deviations in charged quark currents are expected to involve the third generation (a 
notable exception being the flavor universal QiQi structure present already in the SM, which we 
have dropped from our analysis). 

As a consequence of our effective theory approach, operators given in Eq. (3.10) modify 
not only the tWb vertex but withhold a much richer flavor and gauge structure. Since our aim is 
to analyze the effects of these operators in B physics we make two notes regarding the additional 
effects Qi might cause that we shall not be pursuing. 

• Qll and Qim also modify tW.s and tWd vertices. Consequently, they also contribute to 

mixing at one-loop. However, their contributions to neutral kaon and as well as B 
meson oscillations turn out to be universal and purely real (see discussion below Eq. (3.18)) 
so they cannot increase the e and e' predictions. 

• QhRb and Q'l^i^ also modify uWh and cWb vertices. This could interfere with Vcb and Vub ex- 
traction from semileptonic B decays. Since these quantities are crucial for the reconstruction 
of the CKM matrix in MFV models, a consistent analysis of these operators would require a 
modified CKM unitarity fit, which is beyond the scope of this work. 

The Feynman rules for all the vertices generated by Qi that are relevant for our analysis are 
presented in the Appendix A. 2. Since we shall be working in the general R^ gauge, beneficing the 
check of the ^ dependence cancelation in the final results, we will have to consider also the would-be 
Goldstone bosons in our calculations. 

3.3 |A5| = 2 transitions 

Recently, possible NP effects in the Bg — Bg, mixing amplitudes {q = d, s) have received considerable 
attention (c.f. [64] and references within). In particular within the SM, the B^ — B^ mass difference 
and the time-dependent CP asymmetry in B^ — )• J/ipKg are strongly correlated with the branching 
ratio Br(i?"'' — )• r+z/). The most recent global analysis point to a disagreement of this correlation 
with direct measurements at the level of 2.9 standard deviations [6 I]. 

Similarly in the Bs sector, the measured CP-asymmetries by the Tevatron experiments, namely 
in Bg — )• J/^(f> and in di-muonic inclusive decays when combined, deviate from the SM prediction 
for the CP violating phase (ps in Bs — Bs mixing by 3.3 standard deviations [(34]. This indication of 
NP effect is however weakened by the latest LHCb result of the (ps phase inferred from combined 
analysis of Bs — s- J/ip(j) and Bs — s- J/ 1^/(980) channels [109], showing agreement with SM prediction 
therefore eliminating the possibility of large NP contributions. 

In this section we analyze the effects of our operators (3.10) on the Bg — Bg mixing amplitudes. 
We first perform the matching to the low energy Lagrangian, where we consider only diagrams with 
one Qi operator insertions, resulting in first order corrections in 1/A^ expansion. Consideration 
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of only single NP operator insertion is a good approximation given the small size of observed 
deviations in the CP-conserving Bq mixing observables from SM predictions. However, we have 
also computed higher order insertions and checked explicitly that they do not change our conclusions 
of the numerical analysis presented in section 3.3.2. 

To obtain the constraints on NP contributions we rely on the recent global CKM and Bq mixing 
fits given in Refs. [G l] and [110]. 

3.3.1 Matching 

In order to study the effects of our operators (3.10) on the matrix elements relevant in Bq — Bq 
mixing, we normalize them, following [64], to the SM values given in Eq. (1.40) by writing 

M!, = ^{B'q\n,s\B^q)d,sp = Mli^Aq , (3.12) 

where the deviation of parameter Ag from 1 quantifies NP contributions. Proceeding in a similar 
fashion as in section 1.4, where we have analyzed the mixing amplitudes in the SM, we now match 
our effective theory to the low energy effective theory relevant for |Ai?| = 2 transitions which is 
governed by the Lagrangian 

Afr = -^^{Vt,v:qfjzc^{^^)o^. (3.13) 

i=l 

Compared to the low energy effective Lagrangian (1.42), to which we were matching the pure SM 
contributions, effective Lagrangian of (3.13) contains four additional operators [111] 

Oi = [d%bl] [d^X] ' = [d%b^] [462] , (3.14) 

of = [462] [^6^] , Oi = [463 [^6^] , 

which need to be included since non-SM chirality structures are present in our operator basis. We 
have explicitly written out the a, f3 color indices. 

In the matching procedure the W boson and the top quark are integrated out by computing 
the box diagrams such as the one depicted in Fig. 3.1, which now contain anomalous couplings, 
from insertions of operators Qj. The box diagrams with anomalous couplings appearing in the 
bottom-right corner instead the top-left and the crossed diagrams with internal quark and boson 
lines exchanged are completely symmetric and need not be computed separately. 

We note that working in general gauge for weak interactions brings about new anomalous 

interactions of Would-be goldstone bosons generated by Q^^'^ and Qrr operators. What is more, 
in the general R^ gauge, operators Qll and Q'^^ contribute to mixing amplitudes also through the 
triangle diagrams shown on the righthand side of Fig. 3.1. 

By simple consideration of the chirality structure of the diagrams we find that one insertion 
of operators Qrr and QiRb give contributions suppressed by the down quark masses. This is 
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Figure 3.1: Feynman diagrams for Bq — )• Bq transitions with one insertion of Qj operators, labeled 
with a square. The zigzag lines represent W gauge bosons or would-be Goldstone scalars (p. Quarks 
running in the loop are up- type quarks. Left: Box diagram to which all Qj contribute. Right: 
Triangle diagrams generated only by Q^'l operators. 




Figure 3.2: Schematic consideration of Qrr and Qim insertions in the mixing diagrams. Crosses 
and qs on the top quark propagators mark which part of the propagator is needed due to specific 
chirality demands of the vertices, q at the anomalous vertices mark that the Feynman rules for 
those vertices include the loop momentum q. 

illustrated in Fig. 3.2, which shows explicitly the chiralities of quark fields and the number of times 
the loop integration momentum q appears in the mixing diagrams. Because we neglect the down 
quark masses, no chirality flip of these fields is possible, therefore their chirality is determined by 
the interaction vertices. The up (top) quarks in the loop however can experience a helicity flip 
or instead pick up the loop integration momentum q from the propagator. What is more some 
anomalous couplings also include momentum q in the Feynman rule. If we sum the number of 
times that q appears in the diagram to be odd, such diagram gives a zero contribution due to 
the symmetry of the integration over d'^q and the fact that the only momentum appearing in the 
diagram is q (see section 1.4.1) 

/ r9Urrr2 21 =0^ n = 0,l,2,.... (3.15) 

This is what happens in all one-insertion diagrams for Qrr and QiRb, which are presented in 
Fig. 3.2, meaning that these two operators contribute only upon two insertions and will therefore 
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not be considered further in this section. In turn the only relevant four quark operator of the 
effective Lagrangian (3.13) for the matching is Of and we can write 

Ci = Cf^ + 6Ci, (3.16) 

where Cf^ is given in Eq. (1.51) and 6Ci, which is identical for q = d and q = s case is 

5Ci = Re[KLL]S^'^{xt) + Re[KLm]S^''\xt) + 4^ S^'^'^'^xt) + K'^mS^''" {xt) , (3.17) 

where xt = ml / uiy^ and we have defined 

_ ^LL (') _ ^LRt (n ■|o^ 

" A^V2Gf ' " A^Gf ■ 

The 5'o(x() loop functions are given in Eqs. (B.ll) of the Appendix. Their gauge independence has 
been checked by the cancelation of all .^-dependent terms. There are two features of the obtained 
results that we want to point out. 

• First one is that, Sq^ and Sq^" contributions turn out to be UV-divergent. The divergences 
originate from the triangle diagrams with two would-be Goldstone bosons running in the 
loop'^. We renormalize them using the MS prescription, leading to remnant renormalization 
scale dependent terms of the form xt log / /i^ . Because of this ultraviolet renormalization, 
it would be inconsistent to assume that no other operators but those in Eq. (3.10) comprise 
the dimension-six part of the Lagrangian (1.66). In particular, on dimensional grounds it 
is easy to verify that the appropriate MFV consistent counter-terms are generated by the 
four-quark operators of the form 



giving the change in Wilson coefficient 



Appearance of xt in the expression above is crucial to match the xt factor accompanying the 
UV divergences in Qll and Q'[]^ contributions to SCi. 

Generic tree- level contributions of this kind to 5Ci have been analyzed in detail in [112] 
although not in the context of radiative corrections but as standalone dimension-six AF = 2 



^Had we been working in the unitary gauge, these diagrams would not exist. The UV divergences would however 
still emerge from the box diagrams, since one Qll or Q'll insertions demand one of the up-type quarks to be the 
top quark. This means that rather than having two summations in Eq. (1.44) as we did in the SM case, we have just 
one, which leaves the results UV divergent. 
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effective operators adhering to MFV - we will not consider them further. It is however 
important to keep in mind that our derived bounds on presented in the next section 
assume that the dominant NP effects at the /x ~ mt scale are represented by a single Qi 
insertion. 

• The second feature to be pointed out is that only real parts of ull and klui enter Eq. (3.17) 
and thus cannot introduce a new CP violating phase. On a computational level, this is due 
to the fact that these operators always contribute to the mixing amplitudes in hermitian 
conjugate pairs. In particular, a box diagram with an insertion of operator Qll or Qlhi in 
the upper-left corner is accompanied by a diagram with insertion of Q|^^ and Q|^^^ in the 
upper-right corner, since as was pointed out in section 3.2.3, these operators also effect tWd 
and tWs vertices. Both diagrams give the same result but one with HLL,LRt and the other 
with K^L LRt pre-factors respectively, resulting in the appearance oiRe[KLL,LRt] in the sum of 
both contributions. Similarly, the triangular diagrams generated by Qll are also generated 

by qL- 

This inability to introduce new phases can also be understood more generally already at the 
operator level. Namely as shown in Ref. [113], a necessary condition for new flavor violating 
structures yx to introduce new sources of CP violation in quark transitions is that 

Tr{yx[{Y^Y^),{YaY^)])^0. (3.21) 

In MFV models (where yx is built out of Y^ and Y^ ) this condition can only be met if yx 
contains products of both Y^ and Y^. In our analysis this is true for all operators except Qll 
and QiRt- 

3.3.2 Semi-numerical formula 

In order to evaluate the hadronic matrix elements of the operators, we have to evolve the Wilson 
coefficient (3.16) from the matching scale at the top quark mass to the low energy scale at the 
bottom quark mass. Next-to-leading log (NLL) running for the SM Wilson coefficient Cf^ in the 
MS (NDR) scheme is [114] 

Cf^irub) = OMOCf^imt). (3.22) 

Because we are relying on results with consistent MS renormalization procedures, we need to use 
the MS quark masses mt = rnt{mt), mf, = rni,{m.{j). Following the reasoning outlined in the last 
paragraph of section 1.5 we assume the same running also for the 5Ci part. Under this assumption 
the effects of running between the change 5Ci and the SM contribution cancel out and the parameter 
for quantification of NP in mixing amplitudes can then be written as 



2.57Re[Kii] - lMRe[KLR.t] + 2.00 k'^l - 1-29 4^ - 0.77 , (3.23) 
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where the parameters Ki are understood to be evaluated at the high matching scale rrit. In order 
to be consistent with the global analysis of Ref. [64], on which we shall rely in the next section, we 
have used the numerical values for masses and other parameters as specified therein. 

3.3.3 Bounds on NP contributions 

Let us first assume that the anomalous couplings Ki are real. Using Eq. (3.23), we consider one 
Ki{fi = mt) at the time to be non-zero. The assumption of real Ki makes our NP fall under the 
"scenario 11" of [64], for which the global analysis gives 



A = 0.90 



+0.07 




+0.31 


-0.07 




-0.10 



here the bracketed intervals represent the la and 2cr C.L. intervals around the central value, which 
we can use to obtain 95% C.L. bounds on Kj. We present our results in Tab. 3.2. Compared to 



95% C.L. 




Table 3.2: 95% C.L. allowed intervals for Ki which are considered to be real and analyzed with only 
one being different than zero at a time. Accompanying graph shows A as function of one Ki. The 
horizontal line represents the central fitted value and the orange bands the la and 2a regions for 
A as obtained in the global analysis of [64]. 95% C.L. limits on the parameters Kj, are obtained by 
looking at where the functions cross the outer orange regions. 



existing B — )• constraints for couplings kll and KLRt given in Ref. [97], we find our bounds on 
kll to be comparable, while bounds on Kim are improved. 

Relaxing the assumption of real Kj, contributions of the primed operators can introduce new 
CP violating phases if the anomalous couplings have nonzero imaginary components. The analysis 
of such general complex contributions to A fall under the "scenario IH" of [64] . As we have already 
mentioned at the beginning of this chapter, the best fitted values for NP when new phases are 
considered given in [64] are highly influenced by the extraction of the (j)s, the mixing angle in the 
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Figure 3.3: 68% (dashed) and 95% (solid) C.L. allowed regions for and k'^^^ in the complex 
plain obtained from Eq. (3.23) and the results of Ref. [110]. 



Bg system, from Tevatron data. One can imagine the fits to change quite significantly on the side 
of imaginary component of A, if the newest measurement of from LHCb [109] is taken into 
consideration since it is in agreement with the SM prediction. This is confirmed by comparing the 
graphical results from Ref. [64] and the recent update given in Ref. [110], which we use to present 
the 68% and 95% C.L. regions in the complex plane for and k^/j^ shown in Fig. 3.3. We can 
still observe the preference of non-zero imaginary components by the global fits, significance of 
which is however reduced compared to the situation prior to the LHCb measurement [93]. 

3.3.4 Summary 

To summarize, we have matched our effective Lagrangian (1.66) to the low energy effective La- 
grangian (3.13) and analyzed the impact of the effective operators Qi on |Ai?| = 2 mixing to first 
order in 1/A^, allowing one Qi insertion in mixing diagrams. 

We have shown that operators Qrr and QhRh do not contribute upon one insertion, and while 
operators Qll and Qliu can not contribute any new phases to the mixing, the primed operators 
Sff and Q'i^j can. 

Effects of Qi turn out to be the same for and Bg systems, so that they can be parametrized in 
terms of one parameter A. Following the global analysis of [64, 110] we were able to put constraints 
on Kj. In particular, first assuming the Wilson coefficients to be real, were able to obtain for 
them the 95% C.L. allowed intervals given in Tab. 3.2, which compared to the 6 — )• 57 constraint 
given in Ref. [97] prove to be competitive for ki^i^ and improved for kl^. For the three primed 
operators, which can contribute new phases, we have obtained the 95% C.L. allowed regions in the 
corresponding complex plains, presented in Fig. 3.3. 
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3.4 |A5| = 1 transitions 

In this section we turn our analysis of the rare \AB\ = 1 processes introduced in section 1.4.2. 
After performing the one-loop matching of our operator basis (3.10) onto the low energy effective 
Lagrangian (1.55), we obtain corrections to the relevant Wilson coefficients. We proceed by calcu- 
lating the effects in the inclusive B — )• Xs^ and B — )• Xsi~^i~ ■ In order to derive bounds on both the 
real and imaginary parts of the appropriate Wilson coefficients we include the experimental results 
not only for the decay rates but also for the CP asymmetry in — t- Xgj. After performing a global 
fit of the Wilson coefficients, we derive predictions for several rare B meson processes: Bg — )• fj,^fi~, 
the forward-backward asymmetry in B ^ K*£^i~ and the branching ratios for B — )■ K^*^v9. 

3.4.1 Matching 

The procedure of matching closely resembles the one described for the SM case in section 1.4.2. We 
are again interested in NP contributions to the observables at order 1/A^ and thus only consider 
single operator insertions. 

Generic penguin and box diagrams with anomalous couplings are shown in Fig. 3.4 and Fig. 3.5, 
where again due to the commitment to R^ gauge, we are faced with diagrams with would-be 
Goldstone bosons. Exact diagrams for a specific Qi can be reconstructed using Feynman rules 
given in the Appendix A. 2. 




Figure 3.4: Types of Feynman diagrams encountered when computing h — ?■ sV transitions, where V 
stands for 7, Z, g. Dotted lines represent would-be Goldstone bosons, crosses mark additional points 
where V can be emitted in one-particle-reducible diagrams and square represents an anomalous 
coupling. Gluon emission is only possible from quark lines and with the SM coupling. Quarks 
running in the loops are up-type. 

As the result of the matching procedure we obtain deviations from the SM values for the Wilson 
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Figure 3.5: Box Feynman diagrams contributing to 6 — t- si'^i~ and b — sviy transitions. Diagrams 
with would-be Goldstone bosons are absent, since the leptons are treated as massless. Quark 
running in the loop is up- type. 



coefficients, Cj = Cf^ + 6Ci, which we parametrize as 

<5C,(m) = + (3.25) 
j 

where index j runs over the operator basis (3.10), // is the matching scale and kj are rescaled 
Wilson coefficients defined as'^ 



,('/') _ ^LL ^ _ Crr _ CLRb (/) _ CLfif (o 9c^ 

A2V2G^' ''^^ - A22V2G^ ' '^^^^ " A^G^ ' "^^^ " A^G^ ' ^^"'^^ 

Separate track is kept of Qi and Q| contributions that are quantified by functions /• and 
respectively of which the analytic expressions are given in the Appendix B.3. 

We note, that the matching procedure for operator Q'^ stands out compared to the other 
operators. The charged current structure of this operator resembles that of the SM operator 
Qn'^T'^QiW^. Consequently, when considering Q'j^j;^, Wilson coefficients Ci^...^8 are changed in a 
trivial way Ci{fi) = (1 + «;'^^(^))Cf ^(/u). The change of the remaining Wilson coefficients Cg^iCi/p, 
matching of which involves the Z boson, is however not of this form. 

As in the case of \AB\ = 2 process, some of the diagrams in Fig. 3.4 are UV divergent. We 
remove these divergences using the MS prescription leading to remnant log(m^/;U^) terms. We 
shall quantify the matching scale dependence of our results and consequently their sensitivity to 
the UV completion of the effective theory by varying the scale between fj. = 2mw and fj, = m^y. 
Since UV renormalization is necessary, again our operator basis needs to be extended to include 
operators that can serve as the appropriate counter-terms. Within the employed MFV framework 
examples of these operators read 

• = [Qa^'^AQdY^T'^d] cPdW^, , Q^' = [Qi^AqqQ] [h^i] , (3.27) 
• = [QYt-A'qqQ] [<pyiD,cp,] + [Q-z^A'qqQ] [4>^,iD^M . 

The operator Q^ *' produces a counter-term for divergences in 5Ci , while Q^' 3 provide counter- 
terms for divergent parts of SCg^io^up. The operator Qg*- generates a tree-level bZs vertex. The sets 



'For easier reading we repeat the definition of k^'^'' and fi:^^^ previously defined in Eq. (-i. !•'"') 
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of flavor matrices needed to match the structures of divergencies generated by the various operators 
in (3.10) are 

Agd = YuY^, (3.28) 
Aqq = Y^Y^,YuY^YdYj, 

^QQ = YuY^ , {Y^Y^f , YuY^Y.y} , {YuY^fY.Yj , Y^Y^Y^Y^YuY^ . 

Just as in the analysis of NP in |Ai?| = 2 processes, we will selectively set contributions of certain 
operators to be nonzero in our numerical analysis. Consequently in the following, we will drop 
the implicit (tree-level) contributions of the operators in (3.27) to SCi, as these have been already 
investigated and constrained in the existing literature [108]. 

All f^^\fi^'^ are found to be ^ independent and a crosscheck with results from Ref. [97] is 
possible for some of them. We confirm their original results for all the operators except QiRbi 
while an updated version of [97] confirms our result also for this operator. 





SM 




* 




k'* 


'^"ll 


//* 


KRR 


KLRb 


K-LRt 


'^LRt 


^'lri 


k'* 


h 


-0.19 


0.45 


0.45 


-0.19 





0.45 





-45.3 


85.5 


-0.13 


-0.17 


-0.15 


-0.17 


h 


-0.095 


0.24 


0.24 


-0.095 





0.48 





-20.2 


54.5 


0.15 


0.05 





0.05 


h 


1.34 


-1.11 


-1.11 


1.35 


0.09 


-1.11 


0.009 








0.64 


0.64 


0.009 


0.64 


fio 


-4.16 


1.48 


1.48 


-4.28 


-0.12 


1.48 


-0.12 








-2.41 


-2.41 





-2.41 


fuu 


-6.52 


2.38 


2.38 


-6.63 


-0.12 


2.38 


-0.12 








-4.25 


-4.25 





-4.25 



Table 3.3: Numerical values of functions fl and at /i = 2m\[r- Numerical values used for the 
input parameters are rnt{2mw) = 165.0 GeV, = 0.231, mw = 80.4 GeV, m;,(2mvy) = 2.9 GeV, 
[Vtfep = 1. All fi values correspond to matching at LO in QCD. 

To quantify the effects of our seven operators on Wilson coefficients (3.25) we present the 
numerical values of evaluated at = 2mw in Tab. 3.3. We see that the contributions of the 
operator Qll and qJ^^^ are identical in all cases, which means that kll can not induce new CP 
violating phases in the Wilson coefficients. Likewise, QiRt contributions to Cg^io.^p are Hermitian 
but this operator can induce a new CP violating phase in Cj^^. Finally at order 1/ , operators Qrr 
and QiRb which contain right-handed down quarks only contribute to Cy^g- These contributions 
are however very significant, since they appear enhanced as mt/rrib (B.23, B.24) due to the lifting 
of the chiral suppression, as already pointed out in Ref. [97]. 

3.4.2 Bounds on anomalous couplings 

Having computed 6Ci in terms of kj, we turn our attention to observables affected by such shifts. In 
particular at order 1/A^, the presently most constraining observables - the decay rates for B — )• Xg^ 
and B — )■ X.si'^i^ are mostly sensitive to the real parts of Ki [115]. While in general both B — )• X^^sl 
channels are complementary in their sensitivity to flavor violating NP contributions [11 G], within 
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MFV such effects are to a very good approximation universal and the smaUer theoretical and 
experimental uncertainties in the later mode make it favorable for our analysis. In order to bound 
imaginary parts of Kj, we consider the CP asymmetry in i? — )• Xg^. Finally, we compare and 
combine these bounds with the ones obtained from Bq — Bq oscillation observables given in Tab. 3.2. 

Real parts 

We consider the inclusive B — t- and B — t- Xsl!.^l~ branching ratios, for which the presently 
most precise experimental values have been compiled in [117, 118] 

Bi[B ^ X,7]s^>i.6 Gev = (3.55 ± 0.26) x 10"^ , (3.29) 
Br[^ ^ Xs/i+^"]iowg2 = (1-60 ± 0.50) x 10"^ . 

Because the SM contributions to Cj(/ife) and the corresponding operator matrix elements are 
mostly real [115], the linear terms in 6Ci, which stem from SM-NP interference contributions 
contribute mostly as Re[(5Cj]. These are the only terms contributing at order l/h?. Therefore, 
the bounds derived from these two observables are mostly sensitive to the real parts of Kj. Using 
results of [115], we have explicitly verified that the small Im[(5Cj] contributions to Br[i5 — t- Xgi'^i^] 
have negligible effect for all operators except QRR^LRb- However, even for these operators Im[Kj] 
are much more severely constrained by Ax^^y, discussed in the next section. Also, using known 
NLO B —7- Xg'j formulae [119], we have verified that Im[(5Cj] contributions to this decay rate are 
negligible. To analyze the effects of 6Ci on the two branching ratios, we therefore neglect the 
small Im[(5Cj] contributions and employ the semi-numerical formulae given in Ref. [120] with a few 
modifications that we specify below. 

• In [120] all predictions are given in terms of 5Ci at the scale fif, = 4.8 GeV. Since we wish to 
check how our results depend on the matching scale //, we express 5Ci{^b) using NNLO QCD 
running [71, 72, 73, 74] as 

5Cj{nb) = 0.627 (5C7(mvF), (3.30) 
dC-ji^ib) = 0.579 5C7(2mH/). 

On the other hand, Cg^io are only affected by EW running and their change with scale from 
2mw to mw is negligible. 

• The authors of Ref. [120] assumed 5C% = 0, which is not the case in our analysis. However, 
LO Cy and Cg (thus also SCj and 5C^) enter both observables in approximately the same 
combination (conventionally denoted as Cj^^ , c.f. [121]). Employing the known SM NNLO 
matching and RGE running formulae [71, 72, 73, 74] we can correct for this with a simple 
substitution in the expressions of Ref. [120] for the branching ratios 5C^ — )• 5C7 + 0.24 iJCs, 
where we have neglected the small difference between the matching conditions at = 2mw 
and /i = mw- We have verified that in this way we reproduce approximately the known 5C% 
dependencies in S ^ X,7 [122] and B X^^+f" [115]. 
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• As pointed out in the previous section, Q'^^ is to be treated differently than the other op- 
erators. Its effects in Oi,,..,8 can be seen as a shift in the CKM factor appearing Vf^V^* — t- 
(1 + K!Li)VtbV*g. Consequently SM predictions for these contributions simply get multiplied 
by the factor of |1 + f^'n}'^ and only dCg^iQ^yy need to be considered separatly. 

Taking all this into account and considering only one operator Qi to contribute at a time, we obtain 
the 95% C.L. bounds on Re[Kj] shown in Tab. 3.4. 

The first column shows bounds obtained from Bq — Bq mixing as analyzed in section 3.3, while 
the "combined" column corresponds to combined bounds from all three observables. For the later 
we also present the results when the matching scale is set to /u = myi/ to check the scale dependence 
of our results. We can see that the bounds obtained change significantly only in the case of Kim 
where lowering the scale to /i = mw loosens the bounds by almost an order of magnitude. We have 
also checked that the B — )• Xg'^ bounds agree nicely with those obtained in Ref. [97]. In the last 
column we present the lower bound on NP scale A obtained from the combined bounds and under 
the assumption that Wilson coefficients are of the order 1 at scale 2mw- 

The bounds on k and k lri, indeed turn out to be an order of magnitude more stringent then 
for the rest of the coefficients. This was anticipated by the numerical values given in Tab. 3.3 where 
very large effects of operators Qrr and Qim on C^ and Cs were observed. 



Imaginary parts 



We have shown in section 3.3 that imaginary parts of primed Wilson coefficients (3.18) can affect 
the CP violating phase in Bq — Bq mixing and nonzero values were found to be favored by the global 
fit of [110]. To constrain imaginary parts of the remaining four operators, which do not contribute 
with new phases in Bq — Bq mixing, we consider the direct CP asymmetry in B ^ Xs^ for which 
the current world average experimental value reads [117] 



A 



TjB ^ X,7) - TjB -> X,7) 
r{B ^ Xsj) + T{B ^ Xsj) 



-0.012 ±0.028. 



(3.31) 



Based on the recent analysis of this observable in Ref. [12,3] we obtain the following semi- 
numerical formula 



0.006 + 0.039(A5'7 - A' 



17) 



(3.32) 



+ 



Re[5C7{2mw)] 



0.008 + 0.051 (A5'7-A57 
0.012(A^7 - AIj) + 0.002 Re[6C8i2mw)] 

- 0.256 + 0.264A78 - 0.023A^7 - 2.7991^7 

- 0.668A?7 - O.OO5AT7 - 0.563A78 + 0.135 



lm[6C7{2mw)] 
lm[6Cs{2mw)] 



The estimated intervals for hadronic parameters A"7, AJ7 and A78 as specified in Ref. [123] dom- 
inate the theoretical uncertainty making it sufficient to use a LO QCD analysis in the perturbative 
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0.04 / 
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0.03 \ 
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A > 0.82 TeV 




nil 

U. ii 

-0.11 
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U.i ( 

-0.04 


U.oi 

-0.30 


0.11 / 

— U.UD \ 


0.10 \ 

— U.UD / 


A > 0.74 TeV 


'^"ll 


0.18 
-0.18 


0.06 
-0.22 


1.02 
-1.04 


0.08 / 
-0.17 I 


0.05 \ 
-0.15 J 


A > 0.60 TeV 


K-RR 




0.003 
-0.0006 


0.68 
-0.66 


0.003 / 
-0.0006 \ 


0.002 \ 
-0.0006 J 


A > 3.18 TeV 


f^LRb 




0.0003 
-0.001 


0.34 
-0.35 


0.0003 
-0.001 


f 0.003 \ 
I -0.01 ) 


A > 9.26 TeV 


ULRt 


0.13 
-0.14 


0.51 
-0.13 


0.38 
-0.37 


0.13 / 
-0.07 V 


0.12 \ 

-0.14 ; 


A > 0.81 TeV 


'^LRt 


0.29 
-0.29 


0.41 
-0.11 


0.75 
-0.73 


0.27 / 
-0.07 I 


0.25 \ 
-0.06 ) 


A > 0.56 TeV 



Combined 95% C.L. intervals for ji = 2m^, 




-0.1 0.0 0.1 0.2 



Table 3.4: Lower and upper 95% C.L. bounds on real parts of individual anomalous couplings kj 
for fi = 2mw and bracketed for fi = mw, where mt{mw) = 173.8 GeV and m\){mw) = 3.06 GeV 
have been used. *The B — )• Xsl'^l~ bounds on ^e[nRR^LRb\ are valid in the lm[K,RR^LRf)] = 
limit; see text for details. Last column shows estimated lower bounds for NP scale A assuming the 
combined bounds and Wilson coefficients Ci to be of the order one. Accompanying plot serves for 
visual comparison of the presently allowed intervals. 

regime. Thus, in addition to the numerical parameters specified in [123], we have used the LO 
QCD running for dCj^s in this observable. Performing a combined analysis of all considered bounds 
on the real and imaginary parts of individual Ki in which we marginalize the hadronic parameters 
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Current and projected 95% C.L. regions 

0.02 



0.01 



^ 0.00 

e 

-0.01 
-0.02 

-0.001 0.000 0.001 0.002 0.003 

R£[k] 

Figure 3.6: 95% C.L. allowed regions in the KjiR(2m\Y) and KLRb{2mw) complex plain (solid). The 
constraints are dominated by Bt[B — )• Xg^] and Ax^-y. Dashed and dotted lines represent bounds 
obtained under the projected sensitivity of Super-Belle [124, 125] measurement with 5 and 50 ab~^ 
of integrated luminosity. 



entering ^4x^7 within the allowed intervals, we can obtain the allowed regions in the complex plain 
of (Re[Kj], Im[Kj]) shown in Fig. 3.6. 

As already argued, the imaginary part of ull does not contribute to the 5Ci and thus remains 
unconstrained. It also turns out that due to the large hadronic uncertainties, the imaginary parts of 
K^^j, K^'^ remain largely unconstrained by Ax^y and Bq — Bg mixing observables still provide the 
strongest constraints (except for Im[K/,/jf] which again remains unconstrained). On the contrary, 
constraints on the imaginary parts of krr and KLRb reach per-cent level, as can be seen in Fig. 3.6, 
where we also illustrate the projected bounds for Super-Belle, assuming the measured central value 
to be the same as given in Eq. (3.31) and using the estimated Super-Belle accuracy given in Ref. [124, 
125]. Finally we note that in absence of the long-distance effects on NP contributions considered 
in [123], Ax^-y would exhibit an even greater sensitivity to the imaginary parts of SCt^s [126], thus 
we consider our derived bounds as conservative. 

3.4.3 Predictions 

Having derived bounds on anomalous kj couplings, it is interesting to study to what extent these 
can still affect other rare B decay observables. Analyzing one operator at a time we set the matching 
scale to /i = 2mw and consider kj to be real. 

We turn once more to the semi- numerical formulae given in Ref. [120], and first consider the 
branching ratio Br[i?s — )• for which CDF's latest analysis yields [127] 

4.6 X 10"^ < Br[Bs < 3.9 x 10"*^ , at 90% C.L. , (3.33) 



--- 5ab-' 
Krr 50 ab"' 
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while the LHCb collaboration reports on the upper limit [128] 

Bt[Bs < 4.5 X 10"^ , at 95% C.L. , (3.34) 

and a CMS, Atlas and LHCb combined limit has recently been made public [129] 

Bt[Bs /i+/x-] < 4.2 X 10"^ , at 95% C.L. . (3.35) 

In addition we explore the differential forward-backward asymmetry Afb(q^) in the Bd — )• K*£^£~ 
decay, for which the latest measurement of LHCb has recently been published in Ref. [130]. Finally, 
following Ref. [131] we analyze the allowed effects of kj on the branching ratios Br[i? — )• K^*''i/i>], 
which are expected to become experimentally accessible at the super-B factories [132]. 





Figure 3.7: Ranges of values for branching ratios obtained as anomalous couplings are varied 
within the 95% C.L. intervals given in Tab. 3.4. We also show the SM predictions (black) with 
1(7 theoretical uncertainty band (dotted) and for the muonic decay channel the lower end of the 
experimental 90% C.L. interval from [127], the 95% C.L. upper bound from LHCb [128] and the 
latest combined LHC upper bound [129]. 



We present our findings in Figs. 3.7 and 3.8. The effects of anomalous couplings kj on all 
branching ratios are similar. There is a slight tenancy of anomalous kj couplings to increase the 
predictions compared to the SM values at the level of the present theoretical uncertainties, with the 
exception of and klri, of which effects are negligible. In particular, none of the contributions 
can accommodate the recent CDF measurement of Br[i?s — )• /x"*"/!"] at the 90% C.L. , while the latest 
measurements from LHC are already starting to constrain k'j^j^ and k^^^ contributions and could 
in future become an important factor to be included in the combined bound analysis. Furthermore 
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we find that the forward-backward asymmetry Ape (9^) can still be somewhat effected by n'j^j^ 
and krr, for which we present the bands obtained when varied within the 95% C.L. intervals in 
Fig. 3.8. While not sensitive at the moment, in the near future, improved measurements by the 
LHCb experiment could possibly probe such effects. On the other hand, the contributions of other 
anomalous couplings all fall within the theoretical uncertainty bands around the SM predicted 
curve. 




1234 5 6 123456 



[GeV-J q- [GeV-J 

Figure 3.8: ^fb(9^) band obtained when varying real parts of k'^^ (left) and k^r (right) within 
the 95% C.L. interval given in Tab. 3.4. Also presented are the SM predicted central value (black) 
with la theoretical uncertainty band (dashed) and the latest measured points with experimental 
errors given in Ref. [130]. 



3.4.4 Summary 

We have investigated contributions of anomalous charged quark currents in flavor changing neutral 
current mediated \AB\ = 1 processes within an effective field theory framework assuming minimal 
flavor violation. We have determined the indirect bounds on the real and imaginary parts of the 
anomalous couplings. In particular, we are able for the first time to constrain the imaginary parts 
of Kjiji and KLjih already at order 1/A^. Taking into account the obtained bounds on real parts 
of Ki we have predicted the magnitude of effects that the operators considered might have on the 
branching ratio of the Bg — )• fJ'~^f^~ decay, the forward-backward asymmetry in S — )■ K*i~^£~, as 
well as the branching ratios of B ^ K^*^i/i> decays. The better knowledge of these (especially the 
potential further lowering of the Br[i35 — t- upper limit) and other recently proposed [133] 

observables in the future could further constrain some of the anomalous couplings. 
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3.5 Helicity fractions at NLO in QCD 

Having exhausted the impUcations of charged quark current operators (3.10) on the B meson 
mixing, radiative and rare semileptonic decays, we turn in this section, to the study of how the 
non-SM tWb interactions induced by these operators influence the W gauge boson hehcity fractions 
in unpolarized top quark decays at NLO in QCD. We aim to confront these effects with the indirect 
bounds obtained for the couphngs from B physics to see how much deviation from SM predictions 
in hehcity fractions may still be compatible with the low energy observations and further examine 
if perhaps LHC and Tevatron measurements might turnout to provide more stringent constraints 
on such NP. 



3.5.1 Framework 

Since in this section we shall be dealing with t — t- Wb decays exclusively, the tWb vertex and 
its deviation from the SM value and structure is the only charged quark interaction of interest. 
Following [81], we can obtain the most general parametrization thereof by considering the following 
effective Lagrangian 



r 9 T 



-/^{aLPl + anPR) - {bRLPi + bLRPR) qu tW^ , (3.36) 



where q is the momentum of the W boson and Pr^l = 1/2(1 ±7^) are the chirality projectors. Note 
that ai includes also the SM contribution ul = Vtb + SiaL. Extraction of the modified Feynman rule 
for tWh vertex is straight forward and is given in the Appendix A. 2. All of the operators (3.10) 
considered in the previous two sections generate these couplings with the following correspondence 

5aL = V^fc^f^'^* , aR = VtbK*RR , biR = -£^^t>LRt ' ^RL = ~^^^tb'^LRb > (3-37) 

where kj have been defined in Eq. (3.18). 

We shall parametrize the main decay channel of the top quark in the following way 

t^Wb 



^^Vr, (3.38) 



where i = L, +, — stands for longitudinal, transverse-plus and transverse-minus as introduced in 
section 1.2.1. 

The r* decay rates have already been studied to quite some extent in the existing literature. 
The tree-level analysis of the effective interactions (3.36) has been conducted in Ref. [134]. QCD 
corrections, however, have been studied only for the chirality conserving (SM) operators. As we 
have argued in chapter 1.2.1 QCD corrections are especially important for the observable J-+ since 
they allow to lift the helicity suppression present at the LO in the SM. Helicity suppression in 
this observable is also exhibited in the presence non SM dipole structure bLR of tWb vertex, which 
is especially interesting since it is least constrained by indirect bounds coming from B physics 
presented in Tab. 3.4. It might therefore have the potential to modify the t — >• Wb decay properties 
in an observable way. 
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3.5.2 Computation 

We compute the 0{as) corrections to the polarized rates F* in the mf, = hmit using the most 
general tWb interaction vertex extracted from Eq. (3.36). The appropriate Feynman one- loop 
and bremsstrahlung diagrams to be considered are presented in Fig. 3.9. We regulate UV and IR 
divergences by working in d = 4 + e dimensions. The renormalization procedure and the fusing 
of virtual and bremsstrahlung contributions to attain the cancelation of IR divergences closely 
resembles the procedure (for the t — )• qZ case) described in chapter 2.3.1, where we have analyzed 
the NLO QCD corrections for FCNC top quark decays. There are however some differences in the 
computation worth pointing out. 




Figure 3.9: Feynman diagrams for next-to-leading order QCD contributions. Square marks the in- 
sertion of the generally parametrized tWb interaction specified in Eq. (3.36) and cross the additional 
point from which the gluon can be emitted. 

First, since there are no contributions of gluonic NP operators, there are no QCD mixing effects 
to be considered. Furthermore, we need to make use of the covariant projector technique introduced 
in section 1.2.1 of the introductory chapter and summarized in Tab. 1.2, to project out the desired 
helicities of the W boson. This makes the analysis more involved. What is more, since one of the 
projectors includes an explicit Ca/^^yS tensor one can envision encountering problems with 7^ in d 
dimensions if naive dimensional regularization, whereby 7^ is assumed to anti-commute with other 
'j^ matrices in d = 4 + e dimensions as well, is used. To avoid conceivable problems we use the 
prescriptions based on 't Hooft-Veltman 7^ regularization that have been derived by S.A. Larin 
and are given in Ref. [135, 136] 

7^75 ^ (l-4a,)^e^.,.,.37"^7"'7"^ (3-39) 

where = Cpasj (47r) is needed because since the anticommutativity of 7^ is violated the standard 
properties of the axial current and Ward identities are also violated and need to be restored by 
additional renormalization (see Ref. [135] for details). 
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3.5.3 The decay rates 

In the mfe = limit, which we are employing, there is no mixing between chirality flipped operators 
and the decay rates can be written as 

r(^'+'-) = |ai|2r(^'+-) + |6Lfi|2r[^'+'-) + 2Re{ai6I^}rit'+'"^ + {L^R,+ ^-) ,(3.41) 

where considering the ur, bm pair can be accommodated by changing the role of transverse plus 
and transverse minus decay widths r+ o T~ as indicated by the bracketed term in Eq. (3.41). 
Analytical formulae for T^^^^ functions are given in the Appendix C.2. We have crosschecked T\ 
with the corresponding expressions given in [.31] and found agreement between the results. 





L 


+ 




unpolarized 




L 


+ 




-pi,LO 
J- a 

-pi,LO 
b 

-pj,LO 
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(l-x2)2 
4(1 -x2)2 
2(l-x2)2 




■pj,NLO /-pi.LO 
J- a 1 a 
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■pj,NLO /-pi.LO 
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0.90 
0.96 
0.93 
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0.93 
0.91 
0.92 


2x2 

2x^(1 - x'^ f 


2(1 -x2)2(2 + x2) 
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Table 3.5: Left: Tree-level decay widths for different W helicities and their sum, which gives the 
unpolarized width. All results are in the m,;, = limit and we have defined x = mw/mt- Right: 
Numerical values for rNLOypLO -^yj^^}^ following input parameters nit = 173.0 GeV, mw = 80.4 
GeV, Usinrit) = 0.108. Scale /i appearing in NLO expressions is set to /i = nit. In addition m;, = 4.8 
GeV. These values are used throughout the section for all numerical analysis. 



The LO (0(0^)) contributions to decay rates L^'^'^j are obtained with a tree-level calculation 
and are given on the left side of Tab. 3.5. Our results coincide with those given in [134], if the mass 
mb is set to zero. 

The change of FJ^ ^ going form LO to NLO QCD is presented on the right side of Tab. 3.5. 

Since in the = limit L^^^^ vanish, we use the full mi, dependence of the LO rate when dealing 
with W transverse-plus helicity. Effectively we neglect the 0{asnrn,) contributions. In Ref. [28] it 
has been shown, that these sub-leading contributions can scale as as(mb/mvi/)^ log(?7ib/mt)2 leading 
to a relative effect of a couple of percent compared to the size of 0{as) corrections in the nn, = 
limit. 

3.5.4 Effects on 

In Fig. 3.10 we present how each separate anomalous coupling, assumed to be real, affects 
Deviation of the left-handed current coupling from the SM value 5aL can not be probed with 
helicity fractions as long as no interference effects with other NP is considered, since its effects 
simply factor out in the decay widths. The impact of going from LO to NLO in QCD is presented 
in terms of the bands where the lower line corresponds to LO while the upper line presents the 
inflated NLO result. 
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A relatively narrow range in the size of anomalous couplings is shown since using Eq. (3.37) we 
translate the indirect constraints on anomalous couplings given in Tab. 3.4 to find the 95% C.L. 
allowed intervals are quite narrow 

- 0.0006 <aR< 0.003 , -0.0004 < bjiL < 0.0016 , -0.14(-0.29) < < 0.08 . (3.42) 

Two separate lower bounds on 6_l_r, which is substantially less constrained the other two anomalous 
couplings, correspond to which operator Qim or Q'un we assume the bLR to be generated from. 
The graph on the righthand side of Fig. 3.10 shows the J-+ dependence on biR in more detail, along 
with intervals given in Eq. (3.42). We see that the increase is substantial when going to NLO in 
QCD, but still leaves J-+ at the 1 — 2 per-mille level. 






^ 1 


NNLO SM 




u 






-J 


LO 




B-physics 
allowed 





-0.10 -0.05 0.00 0.05 



0.0 
^LR 



Figure 3.10: Dependence of on anomalous couplings which are considered to be real and non- 
zero one at the time. Left: Value of for aji (blue, dotted), bm (orange, dashed) and biji (black, 
solid). Lower and upper lines correspond to LO and NLO results respectively and cross marks the 
SM NNLO prediction. Right: Value of J-+ for Blr. Dashed line corresponds to LO results, while 
the solid line represents the NLO results. We also present the SM NNLO value along with its error 
bars given in Eq. (1.17) and the 95% C.L. allowed intervals for bLR given in Eq. (3.42). 

Since the indirect constraints on non-zero values of ur and bRi are very stringent they can not 
produce large both giving the maximal values of J-+ = 0.00133, which is within 1% of the SM 
prediction jtSM.nlo ^ 0.00132. 

3.5.5 Effects on J'l 

Analyzing a single real NP contribution at the time, leading QCD corrections decrease J-'l by 
approximately 1% in all cases. In in Fig. 3.11 we show the J-'l NLO dependance on the anomalous 
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couplings. Possible effects of ur and bjn are again severely constrained by indirect B physics 
considerations. On the other hand, we find that the most recent combined Tevatron measurement 
of Tl given Eq. (1.19) allows to put competitive bounds on Ilr compared to the indirect constraints 
given in Eq. (3.42). A detailed plot of dependance on Blji is given on the right graph of Fig. 3.11. 
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Figure 3.11: Dependence of J^l on anomalous couplings which are considered to be real and non- 
zero one at the time. Also shown is the central measured Tevatron value (dashed) and the 95% 
C.L. interval as well as the expected ATLAS 95% C.L. interval put on top of the Tevatron central 
value. Left: Value of for an (blue, dotted), Brl (orange, dashed) and Blr (black, solid) at 
NLO in QCD. Right: Dependance on b^R and the 95% C.L. allowed intervals for biR given in 
Eq. (3.42). 

We see that at present the indirect constraints are a bit better, however if the projected sensi- 
tivity is reached, the direct bounds from J^l could turn out to be more stringent. 

3.5.6 Comparison with direct constraints 

As we have shown, there is an interesting interplay between direct and indirect constrains when 
considering the anomalous charged quark currents of the top quark. In this subsection we would 
like to stress this point further, going a bit beyond the scope of this work, by including another 
important top quark process that is influenced by anomalous tWb vertices, namely the single top 
quark production which proceeds through weak interactions. 

Not going into details of the subject, we only comment that measurement of the single top quark 
production cross section [137] and its agreement with the SM predicted value serve to constrain NP 
contributions affecting the cross section [138]. Having an additional sensitive observable at disposal 
one can consider pairs of NP operators altering the tWb vertex contributing simultaneously and 
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obtain 95% C.L. allowed regions in the corresponding NP parameter planes. This was performed in 
Ref. [139] using the single top production cross section and the helicity fractions as the constraining 
observables. 




Figure 3.12: 95% C.L. allowed regions in different planes. The gray bands represent the 

allowed regions from direct Tevatron constraints given in Ref. [139] and Ki are assumed to be real. 
Left: krr - kll (solid), k'j^j^ (dashed), hi'[j^ (dotted) plane. Matching scale is set to /i = 2mw. 
Right: Kx^jih - KLRt (solid), k'^^^^ (dashed) plane. Matching scale is set to /x = 2my/ (narrow 
regions) and /i = raw (wider regions). 

We compare the regions presented there with those that we can obtain using our indirect con- 
straints from B physics in Fig. 3.12. The comparison nicely summarizes the interplay of direct and 
indirect constraints on anomalous tWb interactions and shows that they are in a way complemen- 
tary. On both plots the gray area represents the 95% C.L. allowed regions obtained in Ref. [139]. 
They appear as bands because the direct constraints in krr ^'iid KLm, directions are much weaker 
than the indirect. On the other hand we can see that Tevatron constraints on k^^'^ and k^^^^ are 
comparable and in some cases more stringent than the indirect. Having analyzed the helicity frac- 
tion constrains in detail we can deduce that the single top production contributions to the direct 
constraints presented in Fig. 3.12 are significant, improving the direct bounds considerably. 

3.5.7 Summary 

We have analyzed the decay of an unpolarized top quark to a bottom quark and a polarized W 
boson as mediated by the most general effective tWh vertex at 0{as). We have shown, that within 
this approach the helicity fraction Tj^ can reach maximum values of the order of 2 per-mille in 
the presence of a non-SM biR contributions. Leading QCD effects increase the contributions hiR 
substantially owing to the helicity suppression of the LO result. 
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Indirect constraints coming from B physics already severely restrict the contributions of anoma- 
lous tWb couplings. In particular, considering only real contributions of a single anomalous coupling 
at a time, all considered anomalous couplings except b^R are constrained to yield within 2% of 
the SM prediction. Even in the presence of the much less constrained bLR contributions, a potential 
determination of J-+ significantly deviating from the SM prediction, at the projected sensitivity of 
the LHC experiments [3''>], could not be explained within such framework. Based on the existing 
SM calculations of higher order QCD and electroweak corrections [29, 30], we do not expect such 
corrections to significantly affect our conclusions. 

Finally, with increased precision of the J-'l and the single top quark production cross-section 
measurements at the Tevatron and the LHC the direct bounds on and 6aL are expected to 
regnant over the indirect. 
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Concluding remarks 



As we are approaching the end of Tevatron era and have well entered into the exciting era of the 
LHC, the hunt for BSM physics is in full effect. Searching for new particles is by no means the 
only way that LHC can produce new answers and questions. One of the area that we hope it 
will shed some light on, is the flavor problem of the SM. There is no question that this is where 
top quark, with its high mass, plays an outstanding role. Since LHC can be considered a true 
top quark factory, top quark physics is for the first time being probed at high precision. Precise 
determination of top quarks parameters and interactions could serve as a window to observations 
of physics beyond SM. In this work we have considered different aspects of top quark decays and 
how NP, which we have parametrized using effective theory approach, could affect them. 

On one hand we have investigated the effects of perturbative NLO QCD corrections on different 
decay rates of the top quark, something to be considered when dealing with quarks and being 
confronted with measurements of ever increasing precision. In particular, we have investigated the 
branching ratios of t — > q'y, Z decays and different kinematical asymmetries in subsequent three 
body decays t — > qt^t' as well as the main decay channel of the top quark t — > Wh^ paying special 
attention to helicity fractions as observables sensitive to the structure of the tWh vertex. 

On the other hand, we tried to accent the importance of considering the effects in well measured 
observables of meson physics, whenever deviations from SM in top quark physics are present. The 
dominant role that top quark plays in rare processes of meson physics, where it appears as a virtual 
particle, should always be kept in mind. While the analysis of indirect constraints for operators 
generating FCNC top quark decays has already been performed and can be found in literature, a 
comprehensive analysis of indirect constraints on operators generating anomalous charged currents 
with top quark has not and is an essential result of our work. As we have shown, the precise 
measurements of different "top quark sensitive" observables in |Ai?| = 2 and |Ai?| = 1 processes 
put constraints on NP. The significance of some indirect constraints is not expected to be met by 
the direct constraints from the LHC data. 

Whether any deviation from SM predictions in top quark decays is to be observed or not, the 
future measurements arc expected to play an important role in the flavor aspects of constructing 
and constraining BSM models. 
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Appendix A 



Feynman rules 



In this Appendix we present the Feynman rules involving different types of anomalous couplings 
generated by effective Lagrangians defined various sections of chapters 2 and 3. 

A.l Neutral currents 

In Tab. A.l we present Feynman rules for FCNC top quark transitions governed by the Lagrangian 
given in Eq. (2.1). 



Table A.l: Feynman rules for tVq FCNC vertices, is the momentum of the outgoing gauge boson 
and Pr^l = (1 i 7^)/2 are the chirality projectors. 
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A.2. Charged currents 



FEYNMAN RULES 



A. 2 Charged currents 

Feynman rules for vertices generated by operators given in Eq. (3.10) that turn out to be relevant 
for our analysis are shown in Tab. A. 2. We use the following abbreviations, labeling flavor with i,j 

vr = HB^hihi , vr^—vr, (A.l) 

sw 

Cw ~ *W 1 / * r , I* ! ,11* "^ibVtj \ 
VL - -f -VL-7, {'^LLS3r + KllSsj + ULL-rr^} , 

9R = -I^LRb , 

9L = --l^lRthi - l^'hRthihj , 



icy" qv 



9r,l 



Pr,l 




.i2_T/. 



{vR,L;VR,Lh'' 



kxV»;^. Pr,l 



Sw 



2mw 



-^^i3l^i-^R,L)PR,L 




.i9_T/. .p 



Pr,l 





H', 




'K <7^f^tmVtn^PL 



ie KLRbhmhn 2mw 




-16; 



1 „a 



TPl 




Table A. 2: Feynman rules for the anomalous vertices. Indicies i,j and m, n label quark flavor. 
Feynman rule for tWb vertex obtained from Eq. (3.36) is given in Tab. A. 3. 




^73 



7^aR,L 



—^OLR,RL 



Pr,l 



Table A. 3: Feynman rule for the general parametrization of tWb vertex, is the momentum of 
the outgoing gauge boson. 
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Appendix B 

Loop functions 



In this Appendix we present analytic expressions for various loop functions obtained in the calcu- 
lation of one-loop amplitudes in FCNC processes. 



B.l FCNC top decay form factors 

Here we present the form factors defined in Eqs. (2.5, 2.5). Expressions are given in d = 4 + e 
dimensions regularizing both UV as well as IR divergences denoted euv and eiR respectively. Further 
we define 

a = K^)T(l-6/2)/(47rr/2, 
and for the UV divergent form factors add the counter terms denoted by 5. The form factors are 



F^ = 



CiR, eiR euv 
2 



euv 
4 

4 

^IR 



- 11 + -TT^ - 27ri 
3 

5 - 41og(l - rz) 
eiR 

5 - 41og(l - rz) 



21og2(l - rz) + 31og(l - rz) ~ 2U2{rz) - 6 
2 



eiR 



euv 



-21og"(l-rz)+41og(l-r2) 



2Li2(rz)-6 



b 1 



'4mt log(l - rz) , 

- — ^log(l-rz), 
rrit 2rz 



B.l) 

B.2) 

B.3) 
B.4) 

B.5) 
B.6) 
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B.2. \AB\ = 2 loop functions 



LOOP FUNCTIONS 



ag 



mt 



V + a 



+ (v-a) 
- C, 
+ {v-a) 



rz{'^-rz){l + rz) , 2rz{'i-rz), 1 - 7rz + 3r| 
h \4 h 



2r. 



i^-rz) 



(1 - rzY 



i^-rz) 



l-rz 



log rz 



2v + {v + a){h-2) 



I - rz 

where we have defined 



log(rz) - ITT ^77:; -^^5 /i + j-. 



i^-rz) 



2(l~rz)3 



{l-rz)' 



(B.7) 



(B.8) 



rz = rnz/m^, v = - 2sni9wQ ■, a = T3 



(B.9) 



For the up-type quarks Q = 2/3 and T3 = 1/2. Further we have introduced auxiliary functions /i 
and /2 for shorter notation 



/i 

/2 



2. i^arctan 



rz V V 4 - ?'z 

— 2Li2(rz — 1) + 2arctan f —-i —\ arctan ( — W 

\3-rz \ rz ) \2-rz\ 

2Re{Li2((l-rz) 



(B.IO) 



4- rz 



2 l-rz 



rz 



rz 



rz 



W{^-rz)rz))} 



B.2 |A5| = 2 loop functions 

Here we present SM as well as NP loop functions for |AS| = 2 processes defined in Eqs. (3.17). 
Functions Sq^^'^ are UV divergent and the forms presented here are MS renormalized. 



Sl^ix,) 



OMS 



Sl^'^'ixt) 



2'^° A{xt - 1)2 



2 - llxt + ^) j_ ^xflogxt 



2{xt - 1)3 ' 



LLtt, 



OMS^ 



+ 



xt (xf^ + lOxt + 1) 

2(xt-l)2 
(xf — Sxf + 12xj 



+ xt log 



4) log 



xt(xf + 1) 2x1 log 

(xt - 1)2 + {xt - 1)3 



(B.lla) 
(B.llb) 



(B.llc) 
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LOOP FUNCTIONS B.3. \AB\ = 1 loop functions 



B.3 |A5| = 1 loop functions 

Below we present loop functions obtained in the calculation of \AB\ = 1 processes within the SM 
in the general gauge defined in Eqs. (1.58). 



Bo{x 
Bo{x 
Co{x 
Do{x 

Eo{x 



E'o{x 



Function 



X X log x 

2(x- 1) ~ 2{x - If 

2x 2x log X 
+ 



1 



(p{x,C) 



x-1 ' [x-iy 2 

x{x'^ -7x + 6) x{3x + 2) 



2(x - 1)^ 



2(x - 1) 



logx - (p{x,C) 



4, x2 l92;-25 x^(-5x^ + 2x + 6) ^ 

- log X H ^ H ^ z log X + 6(x, $) 



+ - 7a; x2(3x-2) 



12(x - 1)=^ 



2(x-l)4 



logx, 



2 x(x2 + llx-18) x2(4x2 - 16x + 15)logx 

3 ^^"^ 12(x-l)3 6(x-l)4 
x(x2 — 5x — 2) ^ 3x2 ^ 



4(x - 1) 



2(x- 1)4 ■ 



x2(^ - l)(x^ + 7x - S^logx x(x^-7x + 60 x^(6x + _ 7^) ^ 

^""'^^ ~ 4(x-i)2(x-e)2 4(x- i)(x-e) 4(e- i)(x-e)2 



(B.12) 
(B.13) 
(B.14) 
(B.15) 
(B.16) 
(B.17) 
(B.18) 

(B.19) 



captures all the ^ dependance and has the property liuix^o (p{x,^) = lim^^i0(x,^) = 0. It is 
obvious that the following linear combinations are gauge independent 



2So(x) - Co(x) , 2So(x) + Co(x) , Coix) + Do{x) . 



(B.20) 



Next we present analytical expressions for functions f^''^ and fj^''^ defined in Eq. (3.25). For 
shorter notation we further decompose 



(i) 



g 



1 



4s2 



J 10 



1 



4s2 



Functions containing explicit dependance have been renormalized using MS scheme. Below we 
give all nonzero contributions 



(LL) 



f(LL) _ ALL,,) _ 22x^ - 153x2 + 159x - 46 ^ 3x' 



2x^ 



[LL,) 



8x^ + 5x^ - 7x 3x^ - 2a;2 



4(.T - ly 



logx, 



2A{x- 1) 



4(a;- 1) 



4 log a;, 



(B.21) 
(B.22) 



91 



B.3. \AB\ = 1 loop functions 



LOOP FUNCTIONS 



^{RR) ^ 

ALRb) _ 

J7 — 

(LRt) _ 



{LRt) 



{LRU) 



fi 



{LL) 



ARR) 
ALRb) 

h 

ALRt) 

Is 

ALRt) 

h 

giLL) 

g(LLO 

-giLLO 
niLLll) 



-{LL») 



ALRt) 



mt 

mt 

mw 
rub 
mt 

mw 



-5x2 _^ _ 20 2x- 



12(x- 1)2 



2(.T - 1)3 



6x-^ - 31.i 



log a; 

- 19.r 



logx 



1 , m^ 

o log— ^ 
8 



+ 



12(x-l)2 6(x-l)3 
-9x3 _^ g32,2 _ Qi^^ iQ 32.4 _ i2x^ - 9x2 _^ 20x - 8 



+ 



_ ALRti) _ mt 

— J7 — 

mw 



48(x - 1)3 
-Sx^ + 17x2 -Ax -4 2x - 3x2 



24(x - ly 



logx 



mt 

mw ' 



\Vti,\' 



24(x - 1)3 

X x2 log X 



+ 



4(x - ly 



logx 



_ f(LL) _ ALLii) _ 



: 4(x-l)3. 

5x3 _ 92,2 _^ 30a, _ 8 3x2 ^ 



24(x - 1)3 



4(x - 1)4 



{LL/) _ — x^ + 5x2 + 2x 3a;2logx 

4(x- 1)4 ' 



mt 
mb 

mw 
rub 
mt 

mw 



(X- 1)3 
^2 



X — 4 ^ 3x log X 



4(x-l)2 ' 2(x-l)3 
x2 + 5x ^ 2x^ - 6x2 _|_ ^ 



4(x- 1)2 
r3x2 - 13x 



2(x- 1)3 
4 5x- 2 



fi 



{LRU) 



{x-lf 
mt 
mw 



4(x- 1)4 



logx 
logx 



5x 



2 3xlog(x) 
8(x- 1)3 ^ 4(x- 1)4 



^(ii') = (_a;__)log 



m 



w 



+ 



250x3 - 384x2 + 39x + 77 
108(x - 1)3 



+ 



+ 



-18x5 + 48x4 - 102x3 + 135x2 - 68x + 8 



18(x- 1)4 



4 X m 
9-2)^°^7i2 



^ 125x3 - 253x2 + 138x 



logx, 
- 16 



36(x- 1)3 
9x5 _^ 12x4 - 48x3 + 99x2 _ 593. + g 



18(x- 1)4 
= /i(^^')=/^f'^ = -flog 



logx-|14i,|'2 



+ 



-93;" + 12.1; ^ - 48,]--' + 99.i;2 



+ 



18(x- 1)4 
h^'^'^"^=~fir'^ = \Vtb\\-l\og 



+ -{l-Xogx-\Vtb?), 

250x3 - 384x2 + 39x + 77 
108(x - 1)3 

l0gX-|yt6| -logx, 



= g{LRt) ^ 

mw 



X \ 
/.^-2^°S"J' 

99x3 ^ ;^35^2 ^ 25x - 50 24x3 - 45x2 + 17x + 2 , 

H T^, log X 



72(x - 1)3 



12(x - 1)4 
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B.3. \AB\ = 1 loop functions 



JLRti) _ '^t 



g(LRti) 



j^(LRt) 



mw 
mt 
mw 



\Vtb\ 



a;2 + 3x - 2 x - 2x'^ 



logx 



?(.T-1)2 4(,T-1)3 

54.r'' + 59.1- + 35.1- - 34 15a;'' - 27.i;- + lO.i: + 1 



hiLL) 



36(;e - 1)3 
x"^ + Sx-2 X - 2x2 



6{x-lY 



logx 



(X-1)2 4(X-1)3 



logx 



llx - 5 -2x3 + x^ - 2x 



X , mt 

2 - 



3x 



w 

2 ' 4(x-l) 

x3 + X^ + X 



+ 



2(x- 1)2 



logx, 



2(x-l) 2(x-l)2 



iogx-|y*,p^, 



TO 



llx -5 x^+x'^+x 



4(x-l) 2(x-l)2 



logx - \Vtb\'^ -logx, 



+ 



3x 3x log X 

2(x - 1) ^ 2(x- 1)2 
3x log X 



(x-1) (X-1)2 ' 

3x 3x log X 



(X- 1) (X- 1)2 ' 

3x log X 



+ 



(x-1) (X-1)2' 



j^(LRt) ^ ^{LRtl) ^ ;j(iiJt) ^ 



3 3x log X 

x-1 ~ (X- 1)2 ■ 



(B.41) 
(B.42) 
(B.43) 
(B.44) 
(B.45) 
(B.46) 
(B.47) 
(B.48) 
(B.49) 
(B.50) 
(B.51) 
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Appendix C 

Decay Widths 



In this Appendix we present various decay widths obtained in calculating rates at NLO in QCD. For 
the FCNC processes covered in chapter 2 we separately present the decay widths for t — )• qZ, t ^ q'j 
including the virtual NLO QCD corrections and the decay widths for bremsstrahlung processes of 
t — >• qgZ, t — )• qg^. For the main decay channel analysis given in chapter 3 we present only the 
combined virtual and bremsstrahlung decay rates. 



C.l FCNC decays 



C.l.l Virtual corrections 

Here we present the decay widths defined in Eq. (2.17) with Cp = 4/3 

8 6 



p7,virt. 
^ 6 



p7,virt. 



-p7,virt. 



nZ.virt. 



a(o) 



b 

An 

(0) 



^ + 7--+21og^ 

efn eiR 3 



CfQ 



- 11 

- 2tt 



^ pZ(0) 



-161og(l - rz) 



l+2rz 



161og2(l-r2) 



2(5 + 8rz) 7r2 2(6 + Tr^) 

log(l - rz) 



l + 2rz 3 
8 -161og(l-rz)- 



8 

2+rz 



+ 10 



■^IR 



16 log- (1 - + log(l - - + 21og (-j 

2 + 3 V 



4Li2(rz) 



2(7 + 6rz) 
2 + rz 



(C.l) 

(C.2) 
(C.3) 
(C.4) 

(C.5) 



-4Li2(rz) 
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C.l. FCNC decays 



DECAY WIDTHS 



ah 



nZ,virt. 



ag 



nZ,virt. 



bg 



'-bg 



-.Z(O) 
ab 



"IR 



-161og(l - rz) 
eiR 



22 
3 



161og^(l - rz) ^-^-^.^ log(l -rz)-^-j- ^U,{rz) 



3rz 



2v log J 



1 



2 a(2-rz)+f;(l-2r^) , a,, , 

3 1 - 3 



2{)log( -L]+{v-a) 



rz 



rz 



Uv 

T 

log(rz) 



4/2 



(l-rz)2(2 + rz)J 



a(4 + rz -r|) -'0(3 + rz)rz Al + ^rz 



6 



(l-rz)(2 + rz) 
rl°)^C;^(«-a)(-7r), 



rz 



rz 



.(0) Os 

47r" 



(C.6) 
(C.7) 

(C.8) 

(C.9) 
(C.IO) 



C.1.2 Bremsstrahlung 

Below we give the t — ?• qgZ, 7 bremsstrahlung decay rates, where for the photon channel Eqs. (C.19, 
C.20, C.21) we include the kinematical cuts. For the purpose of shorter notation we define x = 6rc 
and y = 2E'i^^ /mt. 



p7.brcms. 

^ a 

T-iZ.brems. 



-iZ.brcms. 





8 




25 

^ y 




" 8 

L^IR 


4' + ^'^^log(l , 
l + 2rz 


^z(o) as 


" 8 


L<^IR 


4 „ ^ log 1 7 

2 + rz 


pZ(o) as 


" 8 

L^IR 


44 

^ log(l rz) - 



6 

eiR 



37 



161og(l-rz)- i4 



2rz 



6 



eiR 



(C.ll) 
(C.12) 

+ 161og2(l - rz) ~ 41og(rz) log(l - rz) (C.13) 



4(1 -rz- 2rl)rz , „2 



(l-rz)2(l + 2rz) 
161og(l - rz) + ^ - 10 



log(rz) -TT-" -4Li2(rz) + 



7 + rz 
(1 -rz)(l + 2rz) 



10 



2+rz 



eiR 

4(2 - 2rz ~ r|)rz 
(l~rz)2(2 + rz) 



+ 161og2(l - rz) - 41og(rz) log(l - rz) (C.14) 



log(rz) -TT^ -4Li2(rz) 



8rz 



43 



22 



^"^^"^^^ ''^^ ^ + 161og2(l - rz) - 41og(rz) log(l - rz 

eiR 

4(3 - 2rz)rz 



(l-rz)(2 + rz) 3 

(C.15) 



3(l-rz)2 



log(r2)-7r2-4Li2(r2)- 



3(1 -rz) 



47 

y 
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C.l. FCNC decays 



-iZ,brems. 
ag 



+ 



-iZ,brems. 

bg - 



,Z(0) 
ab 



2v 



-(3 - 4rz + r|) + log(rz)(l - rz - r|) - Li2(l - rz) 



rz-2 



)) 



3(l-rz)2 47r 

rz \/(4 - rz)rz ( arctan ( . / ^ ) + arctan ( ^ — 

' \ \\A-rz'' ^^{j^-rzYz 

+ 2Re{Li2 (^(1 - rz)(2 - rz - 1^(4 - rz)rz) } 

+ 2a ^(3-8rz + 5r|) + ^log(rz)(-2-7rz + 2r|)+Li2(l-r 

rz-2 



(C.16) 



(3 - rz)v'(4-rz)rz( arctan ^) + arctan 
2Re|Li2 (^(1 - rz)(2 - rz - i\/(4 - rz)rz) } 



.Z(0) 



2t) 



(1 - rz)22(2 + rz) 47r 
4(1 - rz)-\/(4 - rz)rz(^ arctan J 4"^) + arctan 



3(1 - rz)(-25 + 2rz - r|) - 4rz log(rz)(l + rz) 



(C.17) 



rz-2 



\/(4 - rz)rz 



)) -4Li2(l-rz) 



+ 8Re{Li2(^(l - rz)(2 - rz - 1^(4 - rz)rz) } 



+ 2a 



9 - rz(2 + 7rz) + rz log(rz)(8 + 5rz) + 4Li2(l - rz) 

+ 2(2 - rz)V{^-rz)rz( arctan (j /^ ) + arctan f 

V VV4-rz/ ^V(4-'"z)»'z^^ 

- 8Re{Li2(^(l - rz)(2 - rz - 1^(4 - rz)rz) } 
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■p7,brems. 
^6 



■p7,brems. 
'-bg 



+ 



.Z(0) 



+ 6 



(1 - rz)(77 -rz- 4.rl) + 31og(rz)(10 - 4rz - 9r|) 

rz-2 



(1 - rz)22(2 + rz) 47r 

(20 + lOrz - 3r|) (^arctan (^^ ^ ^ + arctan 



(C.18) 



4-rz 
+ 48Re|Li2( 



^(4 - rz)rz 



)) +121og2(rz) 



1 i 

2 + 2 



4-rz- 



(1 - 70rz + 38r| - 5r|) + 31og(rz)(2 + 46rz - 9r| + 41og(rz)) 



— 6(20 — 3rz)-\/ (4 — rz)rz^ arctan ^— ^ — ^ + arctan ^ 



rz - 2 



V(4-r-z)rz 



)) 



+ 48Re 
+ av 



{Li.( 



1 i /4-rz 

2 + 2 



rz 

„2 



)-Li.(!| + 1^(437^)} 



7 + 22rz - 15r| - log(rz)(6 - 5r| + 41og(rz)) 
2(2 + rz) V(4 - rz)rz ( arctan ( 

16Re|Li2( 



rz 
A-rz 



^ + arctan ^ 



^(4 - rz)rz 



)) 



■1 i 
^2 + 2 



4-rz 



rz 



)-Li.(^ + 1^(4^7^)} 



1-y 



+ l-7r^ 



^ y(l-y)(2y-l) 
2 - 



+ y+-(2-y)(l-y) 



- 16 



- 6 log 



- 21og'(l -y) + if + 2y- 10) log(l - y) - 



2x^ - 24i + 32 



log(. 



yx 



2-x 



2 - xy{2 - y) 



+ 2y) log ( 



2 - xy{2 - y) 
2-yx 



) 



p7(0) ^ r) 



(l-y)(2-x)(j/i2 -2yx-2f + 8) 2?: 



x'^{2 — yx) + 3 

2 - xy{2 - y) 



(C.19) 



^(^-^)^°^( (i:,?^:g) )-^^°^(^)^°^( 

^(f ^ - 4x + 4) log + 4(Li. (I) - Li.(^) - Li. ( 



)+--©-^(23|^) 



(C.20) 

xy{2-y)j 



8 arctan (^—^===^ (^2/x — 1 — arctan(-\/2/x — 1^ 

{Li(^ (2 - i - iV(2-i;)i)) - Li2 (2 - - iy ^(2 - x)x)) } 



8Re 
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C.2. Main decay channel 



p7 

g 



j.i{o)as 2 



{1 - y){2 ~ x){3yx^ - ixy ~ 8x + 16) 2tt^ 



x^{2 — xy) 



(C.21) 



+ {A-2x + 41og (I)) log(^) + (3 - y)(l - y) log (^^^) 



(2 - x){x^ -x^ + Qx-^) log + 4(Li2 (f ) - Li2 (|) - U,{y)) 



C.2 Main decay channel 

Here we present analytical formulae for all nine T^'^^^ appearing in Eq. (3.41) to 0{as) order and 
in the rrih = limit. Note that in this section x = mw/mt 



(l-a;2)2 as 



2a;2 



4tt 



(l-a;2)(5 + 47x2 -4x4) 271^ 1 + Sx^ + 2x'' 3(1- 

2x2 3 ^2 J.2 



log(l - x^) (C.22) 



2(l-x)2(2-x + 6x2 + .x3) , , ^, M ^ 2( 1 + x) ^ (2 + X + Gx^ - X^) , , 

^ '-^ ^ log X log 1 - X - ^ '-^ ^ '- log X log 1 + X 

x^ x-^ 

2(l-x)2(4 + 3x + 8x2 + x3)^. ^ ^ 2(l + x)2(4-3x + 8x2-x3)^. ^ ^ ^^^^ ^ 2m / ^ 
^ — ^ ^Li2(x) - ^ ^ ^Li2(-x) + 16(1 + 2x2) jQg^^<) 



x^ 



— 

id — 



2x^(1 -x^)^ + -^Cf 
An 



X 
2^2 



2x^(1 - x2)(21 - x^) + ^4x^(1 + x2)(3 - x^) + 4x^(1 - log (^) 



16x^(3 + 3x2 _ jQg(^^) _ 4^-^ ^2^2(^2 + x^) log(l - x^) - 8x(l - x)2(3 + 3x2 ^ 2^3) log(x) log(l - x) 
8x(l + x)2(3 + 3x2 _ 2x^) log(x) log(l + x) - 8x(l - x)2(3 + 2x + 7x2 _^ 4x^)U2{x) 



^ ab 



8x(l + x)2(3 - 2x + 7x2 _ 42,3)Li2(-x) 

(l-x2)(l + llx2) 



(C.23) 



(1-X2)2 + ^C. 
2(l-x2)2(l + 2x2) 



27r2 
3 



(l-7x2+2x*) + (l-x2)2log(4) 



log(l - x2) - 4x2(7 - x2) log(x) - 4(1 - x)2(l + 5x + 2x2) log(x) log(l - x) 



4(1 + x)2(l - 5x + 2x2) log(x) log(l + x) - 4(1 - x)2(3 + 9x + 4x2)Li2(x) 



4(1 + x)2(3 - 9x + 4x2)Li2(-x) 
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r+ 



47r 



Cf 



1 TT^ 

-(1 - a;) (25 + 5a; + 9x'^ +x^) + —{7 + Qx^ - 2x^) - 2(5 - Ix^ + 2x'^) log(l + x) 



(C.24) 



+ 



2(5 + Ix^ - 2x^) log(a;) - 
(l + a;)2(5 + 7a;2-4a;3) 



(l-x)2(5 + 7a;2+4a;3) 



log(a;) log(l — x) 



(l-x)2(5 + 7x2+4x3) 



lA-iix) 



, , , 5 + 10a; + 12a;2 + 30a;3-a;^-12x^^. , , 
log(a;) log(l + xj H Li2(-x) 



r+ = 



47r 



Cf 



-x(l - x)(30 + 3x + 7x2 _ 22.3 _ 2a;4) _ 47^23,4 _ _ q^2 _^ ^^4^ j^g^-^ _^ 



r+ 



+ 8x^(1 + 5x2) jyg(-^>) _ _ ^^2(^4 + 5a; + Gx^ + x^) log(x) log(l - x) 

- 4(1 + x)2(4 - 5x + 6x2 _ jQg^^) jQg(^ + x) - 4(1 - x)2(4 + 5x + 6x2 ^ x^^U2{x) 

- 4(4 + 3x - 16x2 + 6x3 ^ ]^g^4 _ ^5-jLi2(_a:) 

27r2 



Cf 



2x(l - x)(15 - llx) + ^x2(5 - 2x2) - 2(13 - 16x2 + 3x^) log(l + x) 



- 2(1 - x)2(5 + 7x + 4x2) log(x) log(l - x) - 2(1 + x)2(5 - 7x + 4x2) log(x) log(l + x) 

+ 2x2(1 + 3x2) ^Qg(^^ _ 2(^1 _ a.)2^5 + 7x + 4x2)Li2(x) - 2(3 + 3x - 31x2 + x^ + 12x'*)Li2(-x) 



2\2 



r- = (1-0^2) 



47r 



-Cp 



1 7r2 

- (1 - x) (13 + 33x - 7x2 + x^) + — (3 + 4x2 - ^^^-^ 



^.r ^2 . 4m /^ 2(l-x2)2(l + 2x2) , , 2(1 - x2)(l - 4x2) , ,^ 

2(5 + 7x2 _ 22.4^ log^j.) _ _^ y ^ iog(i _ J.) _ _^ A i iQg(j ^ ^) 

(l-x)2(5 + 7x2+4x3) , ^ , (l + x)2(5 + 7x2-4x3) , 

log(x) log(l - x) ^ log(x) log(l + x) 

(l-x)2(5 + 3x)(l + x + 4x2)^. , , 5 + 2x + 12x2 + 6x3 -x^ -4x5,. , 
^ — ^Li2(a;) H \ai(-x) 

X X 



4(1-X2)2 + |C. 



-(1 - x)(16 - 14x + 22x2 + 18x3 - 3x^ - Sx^) 
o 



- —4(4 + x"*) + 8x2(1 + 5x2) ^Qg^^^ _ 24(1 _ ^,2^2 j^g^-^ _ 3.) _^ 3(1 _ 2,2)(2 - ^2) log(l + x) 

- 4(1 - x)2(4 + 5x + 6x2 ^ ^3^ log(x) log(l - x) - 4(1 + x)2(4 - 5x + 6x2 _ ^3^ log(x) log(l + x) 



4(1 - x)2(12 + 21x + 14x2 + x3)Li2(x) - 4(12 + 3x + 6x3 _ x^^ju^^{_x) + 8(1 - x2)2 log ( 



(C.25) 



(C.26) 



(C.27) 



(C.28) 
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2\2 



" ab 



4tt 

2(l-x2)2(l+5x2) 



277 



2(1 - a;)(9 - 6x + - bx^) - — (5 + 2a;*) + 2a;^(l + Sa;^) log(a;) 



(C.29) 



log(l - x) 



2(1 -x2)(l- 9x2-2x4) 



log(l + x) 



- 2(1 - xf{f> + 7x + 4a;2) log(a;) log(l - a;) - 2(1 + xf{f> - 7a; + Ax^) \og{x) log(l + a;) 

2 

- 2(1 - xf{Vi + 23a; + 12x^)U2{x) - 2(15 + 3a; + Zx^ + x^ + Ax^)U2{-x) + 2(1 - x^f log ("y) 
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Appendix D 

Three-body FCNC top decays 



In this Appendix we present some details of the t — )• ql'^i analysis presented in section 2.4 of 
chapter 2. 

D.l Analytical formulae 

Below we give the complete analytic formulae for the partial differential decay rate distributions in 
terms of our chosen kinematical variables and the expression for functions appearing in FBA and 
LRA. Mostly they are given in unevaluated integral form, as analytic integration, though possible 
in most cases, yields very long expressions. 

D.1.1 Photon mediation 

The double and single differential decay widths are given as 



d-uds 167r3 A4 si 
dr^ _ mt g\v^ (1 - 



dr''' mt gif^ „ 1 
= — ^^^^R., X - s 




(D.l) 



ds 167r3 A4 



Functions and defined in Eqs. (2.23, 2.24) are 



1 r 



+ 9x - 61og(x) - 8 , 
2 

+ 3x - 2x^ + rc^ - - log(4x) . 



(D.2) 



9 L 



13 



(D.3) 



18 



D.l. 2 Z mediation 



We define 
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The double and single differential decay widths are given as 



dr^ 

dsdu 



4 4 

nit gz^ 



1 



A + a, , A — a, ,^ .s 

{l-s-u){s + u) + ^ — {l-u)u (D.4) 



+ + I3)us{u + s) + {B- (3)s{l - s - u){l - u) + {C + 7)5(1 - u - s) + (C - j)us 
dr^ mt ffl^;^ {s-lf 



ds 



^{2s + l) + js{s + 2) + Cs 



(D.5) 



For the sake of shorter notation we first define 



n 



r2 



:i-n) 



[(1 - z)(l - - 2mz]2 + 7|(1 - z^il - n)2 ' 



then present fi functions defined in Eqs. (2.25) in the from of the following integrals 



/•I 1 

fA = dsri-(l + 2s) 



fc= dsris, 
10 



1 



/b= / dsri-{2s + s' 
Jo -J 

/a/37 = ~o/c • 



(D.6) 



The gi functions present in LRA expressions defined in Eqs. (2.25) are more complicated due to 
the fact that the angular variable appears in the resonant factor of the matrix element. So for the 
sake of brevity we define additional functions Gt in which the u integration is performed 



gi 
Ga 
Gb 
Gc 

GajSy 



1 

dzGi- dzGi, i = A, B, C, a/37 , 

J-i 

/ du r2 (1 + 5w + 2uz — z'^ + iiz^) , 
Jo 

/ dn r2 4(1 - M + 2ii2 - 2uz - z^ + 3uz^ - 2u'^z^) , 
Jo 

/ du r2 4(1 + M — 2uz — z'^ + uz^) , 
Jo 

1 du r2 (1 — 3?x + 2uz — z'^ + uz^) . 
Jo 



(D.7) 
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D.1.3 Interference between Z and photon mediation 



The interference contribution between the Z and the photon to the double differential decay rate 
is 



dr 



int 



4 4 



dsdu 



Re< 



s-mz - isjz 
(s — m^)2 + s'^Jz 



167r3 A4 

2Wi{l - s - u){l - u) + 2Wiu{u + s) + W3(l - s - u) + W^u 



In all further computations we neglect the imaginary part in the propagator's numerator 7^ ~ 0.02. 
This means that Re acts only on the model dependent constants . . . , fl and g\ are the 
same as fi and gi^ except that the integration limits are altered due to the di-lepton invariant 
mass cutoff e. In fx the s integration is now in the [e/m^, 1] region, in gx the intervals for z are 
[0, 1 — 2e/mf] and [—1, 0], and for the u in Gx functions u G [0, 1 — 3^]. We further define 



[(l-^)(l-z)-2mz](l-n) 
[(1 -z){\-u)- 2mz? + 7|(1 - zY{l 



u 



The new fi and gi functions defined in Eqs. (2.28) are 

/•I 1 

fwi2 = / dsri{s -mz)-{s + 2) 

Jelnii "J 



fw 
9i 



1 

/ dsri{s-mz)- 

Je/m1 2 



l-2e/m2 



dz Gi 



dz Gi , 



(D.8a) 
(D.Sb) 
(D.Sc) 
(D.Sd) 



Gwi 



Gws 



1-z 



dursil-ufil + z), Gw2 



2e/m| 



dur3u{l + u — z + zu) , (D.9) 



dnrs -(1 + n + z - zu) , Gw4 



dn r3 -u . 
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D.2 Matching to the parametrization of Fox et al. 

Here we present the conversion of £cfr presented in Ref. [41] to the form in Eq. (2.1). Fox et al. 
give a complete set of dimension six operators that generate a tcZ or tcj vertex 



tRW, 



flU 1 



O'rl = 91 
OIr = 91 



Q3rQ2\ [hHD^H) - {D^H)^H 
Q2CF^^H tRB^y , 
Q^a>'^H\cRB^,, 



OIr = itRl'^CR W{D^H) - (D^HyH 



for notational details see Ref. [41]. We have suppressed the addition of h.c. for every operator. 
Keeping only FCNC parts and the VEV of the Higgs field we obtain 



o 



V 



OIr = 9^W'^, 



CL(T>"'tR\ , a 

OIr = ^[94- 9 B>.] [tRl^CR . 
Finally our coupling constants can be expressed as 



OlL = \[9Al-g'B,\ 



V 

RL ^ 9'-^B^^ \cLcr^"'tR 



tLCT^'CR 



^Ir 



"rl 



1 

2 

'-'RR 



C]^^ cos2 Ow - C^Rj. sin^ Ow 



RL ' 



V2 

C^R cos'^ ew - Clj^sin'^ Ow 



V2 



^RL + '-'RL 



V2 
V2 



'-'LR + '-'LR 



(D.lOa) 
(D.lOb) 
(D.lOc) 
(D.lOd) 
(D.lOe) 
(D.lOf) 
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Poglavje 5 

Razsirjen povzetek 



5.1 Uvod 

Glavni cilj teoreticne fizike osnovnih delcev je razumevanje, opisovanje in napovedovanje pojavov, 
ki se odvijajo na najmanjsih eksperimentalno dosegljivih razdaljah. Glavno matematicno orodje 
pri tern je kvantna teorija polja. 

Stalni napredki na teoreticnem in eksperimentalnem podrocju so privedli do oblikovanja teorije 
znane pod imenom Standardni Model (SM), katerega glavni koncepti so bill zasnovani v sestdesetih 
letih prejsnjega stoletja [1]. 

Glavni karakteristiki SM sta enostavnost in izjemna prediktivna moc. Zadnja desetletja je zaz- 
namovalo testiranje narazlicnejsih napovedi te teorije z bogato paleto sofisticiranih eksperimentov, 
najodmevnejsi od katerih je Veliki Hadronski Trkalnik (LHC). 

Kot renormalizabilno teorijo osnovnih delcev in interakcij SM zaznamuje umeritvena grupa, pod 
katero je invariantna, in mehanizem zloma elektro-sibke simetrije s pomocjo Higgsovega bozona 

SUi3)cX SUi2)LxU{l)Y ^ SU{3)cXU{1)q, Q = Y + T^, (5.1) 

kar povsem doloca vsebino vektorskih in skalarnih polj. Fermionska polja, ki jih delimo na kvarke 
in leptone, pa se v SM pojavijo v treh ponovitvah (druzinah) istih reprezentacij umeritvene grupe, 
za katere recemo, da imajo razlicne okuse. Baza v kateri so umeritvene interakcije diagonalne se 
razlikuje od baze masnih lastnih stanj, ki izvirajo iz interakcij skih clenov tipa Yukawa. Pravimo, 
da je Yukawa sektor SM edini vir fizike okusa, ki v kvarkovskem sektorju privede do nabitih tokov, 
ki spreminjajo okus 

^cc = [u^L^d^L] V,,W+ - -| [d.LruiL] V:,W- . (5.2) 
Mesanje med druzinami v nabitih tokovih opisuje unitarna matrika Cabbibo-Kobayashi-Maskawa 
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(CKM) 

Vub 

V = \ Vcd Vcs Vcb I . (5.3) 

Vtd Vts Vtb 

Kljub veliki uspesnosti SM kot teorije velja splosno sprejeto prepricanje, da SM ni dokoncna 
teorija osnovnih delcev in njihovih interakcij. Nenazadnje ne opisuje kvantne gravitacije, ki postane 
pomembna pri energijah reda Planck-ove skale Ap ~ 10^^ GeV. Tudi delcna vsebina SM ne zadostuje 
za opis vsega, kar smo do sedaj opazili v naravi, saj SM nikakor ne pojasni obstoja temne snovi 
in energije [9, 10]. Odkritje nevtrinskih oscilacij nedvomno potrjuje obstoj nevtrinskih mas, ki jih 

V SM ne poznamo. Poleg tega sta tu se dve veliki konceptualni uganki, ki nam dajeta misliti, da 
mora koncna teorija biti se bolj dovrsena. Hierarhicni problem [15, IG] izpostavlja tezko razumljivo 
veliko razliko med elektro-sibko in Planckovo skalo. Problem fizike okusa pa izraza, da SM ne zna 
razloziti parametrov fizike okusa in njihove ocitno hierarhicne ureditve. 

V dobi LHC lahko prvic preucujemo fiziko kvarka top z veliko preciznostjo, saj LHC lahko 
smatramo za pravo tovarno kvarkov top. Kvark top, ki izstopa s svojo veliko maso mt = 173.2 it 0.9 
GeV [20] in veliko razpadno sirino glavnega razpadnega kanala, 

r(* - w,) ^ W^^ ''-''^^"'' ~ 1.5 GeV . (.4, 

ki nam omogoca, da kvark top v razpadih obravnavamo kot prost delec, postaja vse bolj zanimiv in 
dostopen za iskanje nove fizike (NF) onkraj SM. Osrednje vprasanje, ki ga bomo poskusali raziskati 

V tern delu, je, na kaksen nacin se NF lahko manifestira in opazi v razpadih kvarka top. Glavno 
vodilo so redki razpadi, za katere SM napoveduje zelo majhno verjetnost. Posledicno so odstopanja 
od napovedi SM lahko jasno opazljiva. 

Za iskanje NF so zanimivi razpadi kvarka top, ki potekajo preko nevtralnih tokov, ki spreminjajo 
okus (FCNC). Ti razpadi v okviru SM niso mogoci na drevesnem redu, zato so njihova razvejitvena 
razmerja neopazljivo majhna [36, 37] 

Br[t cy] ~ 10"^^ , Br[t cZ] ~ lO"^"* , Br[t eg] ~ 10"^^ . (5.5) 

Potencialna detekcija taksnih razpadov bi nedvomno pomenila prisotnost NF. 

Po drugi strani lahko odstopanje od napovedi SM iscemo tudi v glavnem razpadnem kanalu 
kvarka top. Ta poteka preko nabite sibke interakcij e na nivoju drevesnega reda (5.2) in velja za 
eksperimentalno signaturo kvarka top. V kolikor bi struktura nabitih kvarkovskih tokov odstopala 
od strukture, ki jo poznamo v SM, bi se to poznalo na sucnostni delezih (eng. helicity fractions) 
bozona W, ki nastane pri razpadu. Definiramo jih tako, da razpadno sirino glavnega razpadnega 
kanala razdelimo na tri dele 

T{t ^ wb) = Fl + r+ + r_ , (5.6) 
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kjer L oznacuje longitudinalno, + in — pa pozitivno in negativno transverzalno stanje sucnosti 
bozona W. Sucnostne deleze nato vpeljemo kot = Tl^+^^/T. Ker SM napoveduje zelo 

majhno vrednost J"+ [28, 29, 30, 31] 

j|M = 0.687(5) , Tf^ = 0.0017(1) , (5.7) 

so te opazljivke zanimive za iskanje NF. Zgoraj navedene vrednosti vkljucujejo kvantne popravke 
visjega reda, najpomembnejsi od katerih so popravki kvantne kromodinamike (QCD). Intuitivna 

Slika 5.1: Ilustracija razpada kvarka top v mirovnem sistemu po 
glavnem razpadnem kanalu v limiti mb = 0. Siroke puscice pred- 
stavljajo tretjo komponento spina, tanke puscice pa smer gibanja. V 
brezmasni limiti sucnost in rocnost polj sovpadata. Ker je sibka inter- 
akcija v SM izkljucno levo-rocna, je sucnost kvarka b v omenjeni limiti 
vedno negativna. Tretja slika prikazuje situacijo, ki je zaradi ohran- 
itve spina prepovedana, saj ima kvark top spin 1/2. To nam orse, 
zakaj je napovedana vrednost za J-4. v SM majhna. Ta enostavna 
slike se podre, ko opustimo brezmasno limito ali v proces vkljucimo 
kvantne popravke visjega reda. 




razlaga za majhnost je podana na sliki Fig. 5.1. V primeru, da se izmerjene vrednosti sucnostnih 
delezev ujemajo z napovedmi SM, meritve sluzijo omejevanju NF, v primeru, da bi se izmerjene 
vrednosti bistveno razlikovale od napovedanih (zlasti v primeru signifikantno nenicelnega J^+), pa 
bi to lahko pomenilo odkritje NF v nabitih tokovih s kvarkom top. 

Pri obravnavi NF v nabitih in nevtralnih tokovih, ki vsebujejo kvark top, ne moremo mimo 
dejstva, da kvark top igra zelo pomembno vlogo tudi v fiziki nizjih energij, kjer se pojavlja kot 
virtualen delec. Zlasti v teoreticnih napovedih za redke procese mezonov B in K lahko v primeru 
NF V fiziki kvarka top pricakujemo spremembe. V primeru, da se meritve ujemajo z napovedmi 
SM, lahko sluzijo za postavitev indirektnih omejitev na prispevke NF. Po drugi strani, v kolikor 
bodo natancno merjene v prihodnje, lahko NF postavlja zanimive napovedi. 

Z ozirom na to v primeru obravnave NF v nabitih tokovih podrobno analiziramo implikacije NF 
na opazljivke v procesih mesanja nevtralnih mezonov B {\AB\ = 2 procesi) in njihovih razpadih 
(|AS| = 1 procesi), ki jih opisemo z naslednjimi efektivnimi Lagrangeovimi funkcijami, v katerih 
ni polj, ki bi imela mase vecje od mase kvarka b 



£|AB|=2 



r2 Trj2 



<? 

1 ' 



4r^^ = ^[i:C.(A.0(") + A.0f)) 



10 



V 



j=3 



(5.8a) 
(5.8b) 



Izraze za Wilsonove koeficiente Cj na elektro-sibki skali izracunamo s postopkom ujemanja polne 
teorije, ki vsebuje vse prostostne stopnje, na efektivne teorije, ki jih opisujeta zgoraj navedeni 
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Lagrangeovi funkciji. S pomocjo enacb renormalizacijske grupe, ki izvirajo iz anomalnih dimenzij 
efektivnih operatorjev, nato spustimo skalo do reda mase kvarka b, kjer je z razlicnimi neperturba- 
tivnimi mogoce izracunati matricne elemente efektivnih operatorjev in s tern napovedati amplitude 
za razlicne procese. 

V kolikor NF posega v proces ujemanja na visokih skalah, se ves njen vpliv na procese v fiziki 
B manifestira kot sprememba Wilsonovih koeficientov na visoki skali ujemanja. 

Tudi pri obravnavi NF v nevtralnih in nabitih kvarkovskih tokovih se posluzimo metod efek- 
tivnih teorij. S pomocjo efektivnih operatorjev visjih dimenzij, kljub nepoznavanju fizike onkraj 
SM, lahko sistematicno parametriziramo ucinke NF na omenjene tokove 

Ces = Csu + ^Y.'^i^i + ^-^^ + 0(1/A^) , (5.9) 

i 

kjer CsM predstavlja SM del, Qi pa so operatorji dimenzije 6, invariantni na operacije umeritvene 
grupe SM ter so sestavljeni le iz polj, ki jih vsebuje SM. Pri tern smo uporabili zelo mocen in pomem- 
ben koncept efektivnih teorij, ki nam zagotavlja, da visje dimenzionalne operatorje spremljajo visje 
negativne potence skale nove fizike in so zato njihovi prispevki vedno manjsi. 



Integrating out t, W, ... Integrating out NP 




uij fif A 



Slika 5.2: Shematicen prikaz nasega pristopa k obravnavi NF. Prvi korak predstavlja dolocitev baze 
operatorjev v En. (5.9) s cimer parametriziramo vpliv NF, katere skala A je dalec nad elektro-sibko 
skalo nt, na fiziko kvarka top. Drugi korak pa predstavlja nadaljnje ujemanje efektivne teorije (5.9) 
z efektivnima teorijama (5.S), ki nam omogoca obravnavo vplivov NF na opazljivke fizike mezonov 
B. 



SI. 5.2 shematicno prikazuje naso strategijo analize NF. Prvi korak ponazarja nase nepozna- 
vanje NF manifestirane na skali A, ki je precej visja od elektro-sibke skale, katere vplive lahko 
parametriziramo z efektivno Lagrangeovo funkcijo oblike (5.9) z izbiro ustrezne baze operatorjev, 
glede na to kaksne spremembe v fiziki kvarka top zelimo obravnavati. Ce zelimo nadalje obravna- 
vati tudi vplive NF na fiziko mezonov B, moramo izvrsiti se drugi korak - ujemanje Lagrangeove 
funkcije (5.9) z Lagrangeovimi funkcijami (5.8). 
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5.2 NF V razpadih kvarka top: Nevtralni tokovi 
5.2.1 Uvod 

Razvejitvena razmerja za procese 



t^qV, V = Z,j,g, q 



c, u . 



ki so V okviru SM neopazljivo majhna, lahko v okviru rezsirjenih teorij postanejo obcutno vecja [37, 
38, 39, 40] in potencialno opazljiva na LHC, saj ATLAS ocenjuje sposobnost odkritja omenjenih 
razpadov v kolikor bi bila razvejitvena razmerja vsaj reda ~ 10~^ [JO]. 

V tern poglavju s pomocjo efektivnih operatorjev analiziramo dvo-delcne razpade t — )• qZ, 7 in 
tro-delcne razpade t — t- qi'^i" . Pri dvodelcnih razpadih upostevamo popravke prvega reda v QCD 
in analiziramo tako posledice mesanja efektivnih operatorojev pod renormaUzacijo, kot tudi efekte 
koncnih popravkov, vkljucno s tako imenovanimi "bremsstrahlung" procesi. Pri obravnavi tro- 
delcnih razpadov pa se osredotocimo na iskanje opazljivk, ki bi lahko pomagale pri diskriminaciji 
med razlicnimi oblikami NF, ki poveljuje razpadom FCNC. To nam omogoca bogatesi fazni prostor 
tro-delcnega koncnega stanja. 

Pri parametrizaciji NF, ki generira tZq, t^q in tgq vozlisca sledimo Ref. [37, 81] 



eff 



A2 



4Qf 



+ 



A2 



+ {L ^ R) +h.c. 



kjer so operatorji definirani kot 



L,R 



qL,Rl^tLfi 
qL,RCT^''tR,L 



Q 



LR,RL 



gzZ^u 



qL,RCT'"^tR,L 



QlR,RL = 9sG% qL,Ra''''TatR,L 



(5.10) 



(5.11) 



Pri racunanju popravkov QCD se posluzimo dimezijske regularizacije, s pomocjo katere regular- 
iziramo tako UV kot IR divergence. 

Analiza ucinkov teh operatorjev v fiziki mezonov je bila opravljena v Ref. [41] in je v tern delu 
ne ponavljamo ali nadgrajujemo in s tern preskocimo oba koraka shematicno prikazana na SI. 5.2 ter 
se osredotocimo le na obravnavo pojavov s kvarki top na masni lupini. Iz analize [41] povzamemo, 
da obstajajo operatorji, ki generirajo efektivna FCNC vozlisca s kvarkom top, za katere indirektne 
omejitve na njihove prispevke ne izkljucujejo potencialne opazljivosti FCNC razpadov kvarkov top. 



5.2.2 Dvo-delcni razpadi 

Najprej se osredotocimo na virtualne popravke QCD. Feynmanovi diagrami za obravnavo popravkov 
na nivoju ene zanke so prikazani na SI. 5.3. V prvi vrsti ti popravki privedejo do mesanja operatorjev 
pod renormalizacijo. Z uporabo metod efektivnih teorij lahko izpeljemo enacbe renormalizacijske 
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Of 



-'lrml 




Slika 5.3: Feynmanovi diagrami virtualnih popravkov QCD za razpade t — t- qZ,^. Kvadratki 
oznacujejo delovanje efektivnih opreratorjev podanih v En. (5.11), krizci pa dodatne tocke iz katerih 
se lahko izsevajo bozoni Z in 7. 



grupe, ki nam povezuje vrednosti Wilsonovih koeficientov ovrednotenih na razlicnih skalah. Ker 
operatorji nimajo anomalnih dimenzij, lahko preostalih 6 operatorjev zberemo v dva vektorja, 



= (07, Of, Of )^ 



RL, LR , 



ki se med seboj ne mesata. Matrika anomalnih dimenzij na nivoju ene zanke je 



2^ 



Cf 






Cf 








8C7^/3 (:7f(3-8sL)/3 5Cf - 2Ca 



(5.12) 



(5.13) 



Pogosto zaradi strukture NF, ki generira operator LR, ta lahko vsebuje ekspliciten faktor mase 
kvarka top. Da analiziramo, ali to privede do opaznih sprememb v anahzi renormalizacijske grupe, 
definiramo nov operator 

dLR = {mt/v)OLR. 

Mesanje gluonskega dipolnega operatorja s fotonsikm in operatojem z dipolno Z sklopitvijo 
povzemajo spodnje enacbe 



hl{^it) = r?'^i67(A) + ^(r?-i-r/'^^)6f(A) 



r/«i6f(A) + 



2-^^2 



(5.14a) 
(5.14b) 



kjer /it predstavlja skalo mase kvarka top, rj = as{A.)/as{nt), 1^1 = 4/3/3o, ^2 = 2/3(3q in /3o je del 
beta funkcije QCD na nivoju ene zanke [52]. 

Ob predpostavki, da pod UV skalo nimamo dodatnih barvnih prostostnih stopenj, ki bi mod- 
ificirale beta funkcijo, imamo /3o = 7 za skale nad fif Ce za operatorje LR vzamemo redefinirano 
obliko Olr, se ki^2 spremenita v ki = 16/3/3o, K2 = 14/3/3o- 

Posledice spreminjanja renormalizacijske skale ponazarja SI. 5.4, kjer prikazujemo 



ko67'^(A) = 0, 
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0.200 




A [GeV] 



Slika 5.4: Razmerje \b]'^ {fit)/bf{A)\ kot funkcija A ob predpostavki hj'^ {K) = in ^ut f« 200 GeV. 
Polne crte se nanasajo na operatorje brez eksplicitne mase kvarka top, prekinjena crta pa se nanasa 
na Wilsonov coefficient operatorja 0\i^. bf^ ni prikazan, saj na grafu odstopanje od bf^ ni opazno. 

ki nam pove koliksna b]'^ (fit) lahko generirano na skali mase kvarka top fit — 200 GeV, izkljucno 
z mesanjem operatorjev pod renomalizacijo QCD in prisotnostjo gluonskega operatorja na visoki 
skali A. 

Opazimo lahko, da inducirani prispevki bj v primeru energijske skale NF okoli A ~ 2 TeV, 
znasajo 10% vredosti bf koeficienta generiranega na skali A. Po drugi strani so, zaradi odstevanja 
podobnih prispevkov v enacbah En. (5.14), inducirani prispevki k bf mnogo manjsi (pod 1% na 
prikazanem razponu skale A). Vkljucitev eksplicitnega faktorja mase kvarka top v operatorje teh 
zakljuckov ne spremeni. 

Po analizi renormalizacijskih lastnosti operatorjev NF nam preostanejo se popravki matricnih 
elementov {qj\ Oi \t) in {qZ\ Oi ki jih ovrednotimo na skali mase kvarka top, in bremsstrahlung 
popravki, katerih Feynmanovi diagrami so prikazani na SI. 5.5. Razpadne sirine teh diagramov so 




Slika 5.5: Feynmanovi diagrami bremsstrahlung procesov t — )• qgZ, 7. Kvadratki oznacujejo delo- 
vanje operatorja NF, krizci pa oznacujejo dodatne tocke iz katerih se lahko izseva gluon (v prvih 
dveh diagramih) ali Z, 7 (v zadnjem diagramu). 

istega reda v as kot razpadne sirine dvo-delcnih koncnih stanj. Sestevek obeh prispevkov poskrbi. 
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da je rezultat IR koncen. Razpadno sirino s popravki reda Os parametriziramo kot 

4 2 2 3 



(5.16) 



+ 



+ 



v^mt 

2 2 

V 



2Re{o^*6»}r^^ - 2Im{a^*6^}r, 



kjer F = Z, 7 in a'^ = 0. F^^ bg g ^ drugi in tretji vrstici En. (5.16) povzemajo prispevke gluonskega 
operatorja in so zato odsotni v prvem redu (a^) in se pojavijo sele na redu Ug. Analiticni izrazi 
vseh razpadnih sirin so podani v dodatku C. 

Kromodinamski popravki procesa s fotonom, ki nima mase, v koncnem stanju so nekoliko bolj 
kompleksni kot v primeru bozona Z. Rezultat, ki ga dobimo iz obravnave predstavljenih diagramov 
virtualnih in bremsstrahulg korekcij, je IR divergenten. To divergenco lahko odstranimo, ce v 
obravanavo vkljucimo dodaten diagram za proces t ^ qg s fotonskim popravkom ene zanke. Ker 
pa eksperimentalno iskanje procesa t ^ qj vselej vkljucuje detekcijo izoliranega fotona, k cemer 
omenjeni diagram ne prispeva, razpadno sirino raje regulariziramo z vpeljavo reza, ki zagotovi, da 
sta smeri fotona in lahkega kvarka ali gluona dovolj narazen 6rj = 1 — p.^ • pj/E^Ej, kjer j = g, q. 
Izkaze se, da je odvisnost rezultata od reza 5q znatna in lahko privede do povecanja prispevka 
gluonskega operatorja. Poleg tega vpeljemo se eksperimentalno motiviran rez na energijo fotona 
E"^"*. Razpadne sirine za FCNC razpade kvarka top v foton predstavimo kot funkcije definiranih 
rezov. 
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Slika 5.6: Relativna velikost as popravkov k . Prikazana sta reprezentativna intervala za 6rc = 5r 
and S™*. Konture z doloceno velikostjo popravkov so narisane za 6^ = (sivo, pike), = 
(rdece) and W = —h^ (modro, crte). 

Numericno analizo prikazuje SI. 5.6, kjer so narisane konture konstantne relativne velikosti ag 
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popravkov k Br(t — )• 57) v ravnini rezov. Opazimo, da so prispevki gluonskega operatorja lahko 
reda 10 — 15% celotne razpadne sirine, odvisno od relativne faze in velikosti Wilsonovih koeficientov 
operatorjev Oj^^ j^^. To pomeni, da eksperimentalna meja na Br(t — )• 57) lahko dejansko omejuje 
tako b"' kot . To opazanje lahko dodatno podkrepimo z analizo razmerja T{t — t- q^)/T{t — t- 
qg), kjer obe razpadni sirini izracunamo do reda a^^, v odvisnosti od razmerja relevantnih FCNC 
Wilsonovih koeficientov \b'^ /b^\, kar prikazuje SI. 5.7 za dve reprezentativni izbiri kinematicnih 
rezov. Vertikalna dimenzija se oblikuje ob spreminjanju relativne faze med koeficientoma b"' in b^. 




0.0 0.2 OA 0.6 O.S 



\b'y/hf!\ 

Slika 5.7: Razmerje razpadnih sirin T{t — )• q'y)/T{t — )■ qg) v odvisnosti od absolutne vrednosti 
razmerja relevantnih FCNC Wilsonovih koeficientov \b"' /b^\. Prikazani so rezultati za dve reprezen- 
tativni izbiri kinematicnih rezov. Obmocja so koncnih razseznosti v vertikalni smeri, zaradi neznane 
relativne faze med b'^ in b^. Prikazane crte so za maksimalno pozitivno (polni) in maksimalno neg- 
ativno (crtkano) interference b"'b^. 

V primeru razpada t — )• qZ so popravki QCD precej manj dramaticni. Signifikanco popravkov 
povzamemo v Tab. 5.1, kjer navajamo relativno spremembo v razpadnih sirinah in razvejitvenih 
razmerjih ob prehodu iz reda na red a<j. 

Opazimo, da relativna sprememba razpadnih sirin lahko doseze 10%, sprememba v razvejitvenih 
razmerjih pa je mnogo manjsa. Razlog za to je skorajsnje popolno iznicenje as prispevkov k razpadni 
sirini t — )• qZ in razpadni sirini glavnega razpadnega kanala t — >• Wb, na katerega so razvejitvena 
razmerja normirana. V kolikor bi tak rezultat pricakovali za prispevke a^, je podoben rezultat za 
b^ netrivialen. Vidimo tudi, da pri dolocenih faznih odnosih med b^ in b^ sprememba lahko naraste 
na par procentov. 



^T{t — 7> qg) s kromodinamskimi popravki je vzeta iz Ref. [70]. 
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=b9 = 


=b9 = 


= b^,bs = 


= 0, = b9 


=0,6^ = ¥ 




0.92 


0.91 


0.92 


0.95 


0.94 


Br^LO/BrLO 


1.001 


0.999 


1.003 


1.032 


1.022 



Tabela 5.1: Numericne vrednosti razmerji 

pNLO/pLO BrNLO(i ^ qZ)/BT^'^{t qZ) za dolocene 

vrednosti FCNC Wilsonovih koeficientov. 



5.2.3 Tro-delcni razpadi 

V tern poglavju se osredotocimo na razpad t — )• ql'^ll.~ , kjer FCNC tranzicija poteka preko istih 
fotonskih in Z operatorjev kot v prejsnjem poglavju, nadalje pa se bozon sklaplja s parom nabitih 
leptonov, ki jih detektiramo v koncnem stanju. Osrednji cilj te analize je potencialna diskriminacija 
med razlicnimi FCNC operatorji na podlagi razlicnih kinematicnih opazljivk, ki jih lahko definiramo 
na racun vecjega tro-delcnega faznega prostora. Za nas bodo zlasti zanimive asimetrije definirane 
na podlagi smeri, pod katerimi se gibajo delci koncnega stanja, ki so lahko senzitivne na obliko 
FCNC vozlisca. 

V mirovnem sistemu leptonskega para definiramo Zj = cos 6j , ki se nanasa na smeri negativno 
nabitega leptona in lahkega kvarka. V mirovnem sistemu pozitivno nabitega leptona in lahkega 
kvarka pa definiramo zi = cosOg, ki se nanasa na smeri nabitih leptonov. Obe definiciji ponazarja 
SI. 5.8 




Slika 5.8: Definicija dveh kotov na podlagi katerih vpeljemo dve razlicni kinematicni asimetriji. 
Puscice oznacujejo smeri gibalnih kolicin delcev. 

S pomocjo tako definiranih kotov vpeljemo dve asimetriji 

A. = . (5.17) 

i Zi>0 + 1 Zi<0 

Aj = ApB imenujemo asimetrija naprej-nazaj (FBA), Ai = Alr pa proglasimo za asimetrijo levo- 
desno (LRA). 

Analize asimetrij se lotimo v treh korakih. Najprej si ogledamo razpade, ki potekajo preko 
izmenjave fotona, nato razpade, ki potekajo preko bozona Z, na koncu pa analiziramo se razpade, 
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ki potekajo preko obeh kanalov hkrati, kar se zdi vredno raziskati, saj stevilni modeli NF lahko 
generirajo znatne FCNC razpade tako v fotone kot v bozone Z. Skupno leptonsko koncno stanje 
nam omogoca raziskovanje interference obeh pojavov. Analiticne formule so podane v poglavju 2 
in dodatku D. 

Zaradi izkljucno vektorske sklopitve fotona z nabitimi leptoni, v primeru fotonske mediacije 
razpada FBA ne more zavzeti nenicelne vrednosti. Po drugi strani se LRA izkaze za nenicelno in 
povsem neodvisno od paramterov nove fizike. Vsebuje pa mocno odvisnost od reza na invariantno 
maso leptonskega para ^/e, ki jo prikazuje levi graf na SI. 5.9. Vidimo, da je glede na razlicne 
kinematicne reze, LRA lahko celo razhcno predznacena. V kohkor FCNC proces poteka preko 




-0.15 -0.10 -O05 0.00 0.05 0.10 0.15 -0.15 -0.10 -0.05 0.00 0.05 0.10 0.15 



[GeV] ApB Are 

Shka 5.9: Levo: Odvisnost LRA od reza na spodnjo mejo invariantne mase leptonskega para e 
za primer razpadov preko fotona. Sredina: Korelacija med FBA in LRA za razpade, ki potekajo 
preko bozona Z. Sivo obmocje je dobljeno, ko so dovoljene vse sklopitve, ostala obmocja pa, ko 
so dolocene sklopitve postavljene na nic. Desno: Korelacija med FBA in LRA za razpade preko 
fotona in bozona Z, vkljucno z interferencnimi prispevki. Sivo obmocje je dobljeno, ko so nenicelne 
lahko vse FCNC sklopitve s fotonom in bozonom Z, ostala obmocja pa, ko so dolocene sklopitve 
postavljene na nic. 

bozona Z, tudi FBA zavzame netrivialne vrednosti. Obe asimetriji sedaj postaneta odvisni od 
parametrov NF. S pomocjo zreba nakljucnih vrednosti parametrov NF lahko preiscemo razpon 
vrednosti asimetrij in korelacijo med FBA in LRA. To prikazujeta srednji in desni graf SI. 5.9. 
Srednji graf se nanasa na razpade, ki potekajo izkljucno preko bozona Z, medtem ko desni graf 
vsebuje tako razpade preko bozona Z kot tudi razpade preko fotona in interferencne prispevke obeh 
procesov. 

Vidimo lahko, da so velike vrednosti FBA (|^fb| ^ 0.1) nedosegljive v razpadih, ki potekajo 
preko bozona Z, kjer je le ta omejena na interval ^fb S [—0.12,0.12]. Eksperimentalno izmerjena 
tocka V ravnini (ApBj ^lr) bi lahko sluzila za izkljucevanje modelov, ki generirajo le dolocene vrste 
efektivnih FCNC sklopitev. Ce obravnavamo razpade preko bozona Z in fotona kot nelocljiva se 
razpon dosegljivih vrednosti LRA bistveno poveca v smeri pozitivnih vrednosti. 

Predstavljeno analizo razpadov t — t- qi~^i~ lahko v prihodnosti soocimo z eksperimentalnim 
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iskanjem t — )■ qZ, kjer identifikacija bozona Z poteka preko detekcije leptonskega para, v katerega 
bozon Z razpade. Po drugi strani je, kot smo ze omenili, eksprimentalno iskanje t ^ q'j vezano na 
detekcijo izoliranega fotona. V kolikor bi se zeleli oddaljiti od teh omejitev in se osredotociti na 
leptonsko koncno stanje, bi bila potrebna nova podrobna analiza ozadij takega koncnega stanja. 

5.3 NF V razpadih kvarka top: Nabiti tokovi 
5.3.1 Uvod 

V tern poglavju se posvetimo nabitim tokovom, ki vsebujejo kvarke top, in analiziramo posledice 
odstopanja od SM v le teh. V analizi sledimo konceptu orisanem v uvodnem poglavju in najprej 
analiziramo implikacije v fiziki mezonov B. Sele nato se posvetimo glavnemu razpadnemu kanalu 
kvarka top in sucnostnim delezem, ki so obcutljivi na omenjene spremembe v nabitih tokovih. 

Pri oblikovanju baze operatorjev dimenzije 6, s katerimi razsirimo SM, najprej poiscemo vse 
strukture, ki so invariantne na umeritveno grupo SM in vsebujejo nabite tokove s kvarkom top. 
Nato dolocimo se okusno strukturo teh operatorjev, pri cemer se omejimo na okvir minimalne 
krsitve okusa (MFV) [103, 104, 105], v katerem edina krsitev okusa izvira iz Yukawinih sklopitev, 
kakor v SM. To nas privede do sledece baze sedmih efektivnih operatorjev 



Qrr = Vtb[tBYbR]{4>iiD^(l)d) , (5.18a) 

Qll = [Q'srVQsll'/'d^'i^/.'Ad) -[QW^Qsll'/'di^M'/'-i)' (5-18b) 

Q'll = [Q3r''j^Q3]i<t>y^D^(^d)-[Q3l^Q3]i4iD^(t>d), (5.18c) 

QIl = [Q'3rVQ3](4^"iI?M'/'d)-[Q'37^Q3](0:^iI?M</'<i), (5.18d) 

QLRt = [Q3r''a^''tR]cPuW^,, (5.18e) 

Q'LRt = [Q3rV^tn]KW^,, (5.18f) 

QLRb = [QsrV^'^bfil^rfH^;,, (5.18g) 

kjer smo definirali SU{2)l dublete 

Q3 = {V^bULk, bLf , Qs = QiVu = ih, VudiiV , (5.19) 
kovariantne odvode ter tenzor elektrosibkih polj 

D, = d^ + dw^T'^ + i^^B^Y, (5.20) 



in koncno skalarna polja (l)u,d (y SM (pu = = ^'^'^^d)- 



118 



RAZSIRJEN POVZETEK 



5.3. NF V razpadih kvarka top: Nabiti tokovi 



5.3.2 Indirektne posledice v fiziki mezonov B 

Prisotnost operatorjev (5.18) spremeni izraze Wilsonovih koeficientov za \AB\ = 2 in \AB\ = 1 
procese, saj kot prikazuje SI. 5.10 operatorji vstopijo v mesalne diagrame in diagrame za redke 
razpade mezonov B. Vselej obravnavamo vnos le enega operatorja v diagram, s cimer se konsis- 
tentno omejimo na prispevke NF utezene z 1/A^, visje potence pa zanemarimo. Ko izracunamo 




Slika 5.10: Primer dveh diagramov, kjer prispevki efektivnih operatorjev (5.18), oznaceni z oraznim 
kvadratkom, vplivajo na proces mecanja mezonv B preko skatlastega diagrama in razpade mezona 
B preko pingvinskega diagrama. V oznacuje foton ali gluon, kvarki v zankah pa so u, c, t. 



vse diagrame in s tern opravimo ujemanje med naso razsirjeno teorijo na efektivne Lagrangeove 
funkcije (5.10), lahko parametriziramo efekte NF s spremembo Wilsonovih koeficientov 

dip) = Cf'^{fi) + 5a{^i), (5.21) 
6Qif,) = Y.Kjifi)F^^'\xt,fi) + ^*{fi)FpHxt,fi), (5.22) 

3 

kjer j = 1, 6 in tece po operatorjih (5.18), funkcije F, ki so odvisne od xt = m^/m^r in skale na 
kateri smo opravili ujemanje n, pa so podane v dodatkih B.2 za |Ai?| = 2 procese (oznacnene s S^) 
in B.3 za \AB\ = 1 procese (oznacene s f^). Definiramo tudi renormirane Wilsonove koeficiente 
NF, ki jih homo uporabljali v nadaljnji analizi 

A',") _ ^LL ^ _ Crr _ CiRb (/) _ gjm /t- 

" A^V2Gf' ^'^'^ " A22V2GP ' ''^^■'-A^G^' ''^^^ " A^G^ ' ^^"'^^ 

S pomocjo podrobnih studij, kako spremembe v Wilsonovih koeficientih (5.21) vplivajo na 
opazljivke v fiziki mezonov B, ki so bile opravljene v Ref. [64, 110, 120, 123, 117, 118], lahko 
izpeljemo omejitve na parametre kj. V kolikor se omejimo na realne nj, lahko izpeljemo intervale, 
V okviru katerih se nahajajo parametri kj s 95% stopnjo zaupnja (C.L.), ki jih prikazuje Tab. 5.2. 
Omeniti velja, da so meje za parametra krr in K^Rh za red velikosti ostrejse kot za ostale parametre. 
Meji izvirata iz analize 6 — )• 57 razpada, kjer so prispevki operatorjev Qrr in QiRb-, ki vsebujeta 
desno-rocne kvarke 6, efektivno povecani za faktor mwt/mb glede na prispevke ostalih operatorjev. 

V primeru, da sprostimo zahtevo po realnosti kj in dovolimo, da imajo tudi imaginarno kompo- 
nento, lahko s pomocjo analize efektov v mesanju in krsitve simetrije parnosti in konjugacije nahoja 
(CP) V razpadih b ^ izpeljemo dovoljena obmocja v kompleksnih ravninah parametrov Hj. To 
prikazuje SI. 5.11 za vse koeficiente razen ki^r^ in kll, katerih imaginarnih komponent ne moremo 
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B-B 


B^Xs-i 


B ^ Xs^Ji+^M- 


skupno 


C^{2mw) ~ 1 




U.Uo 
-0.09 


U.Uo 
-0.12 


U.4o 
-0.49 


0.04 / 
-0.09 V 


0.03 \ 
-0.10 / 


A > 0.82 TeV 




nil 
U. ii 

-0.11 


U.l / 
-0.04 


U.oi 
-0.30 


0.11 / 
n nfi I 

— U.UD \ 


0.10 \ 

— U.UD / 


A > 0.74 TeV 


^"ll 


0.18 
-0.18 


0.06 
-0.22 


1.02 
-1.04 


0.08 / 
-0.17 I 


0.05 \ 
-0.15 J 


A > 0.60 TeV 


K-RR 




0.003 
-0.0006 


0.68 
-0.66 


0.003 / 
-0.0006 \ 


0.002 \ 
-0.0006 J 


A > 3.18 TeV 


KLRb 




0.0003 
-0.001 


0.34 
-0.35 


0.0003 
-0.001 


f 0.003 \ 

I -0.01 ; 


A > 9.26 TeV 


ULRt 


0.13 
-0.14 


0.51 
-0.13 


0.38 
-0.37 


0.13 / 
-0.07 I 


0.12 \ 
-0.14 ) 


A > 0.81 TeV 


'^LRt 


0.29 
-0.29 


0.41 
-0.11 


0.75 
-0.73 


0.27 / 
-0.07 \ 


0.25 \ 
-0.06 ) 


A > 0.56 TeV 



Tabela 5.2: Zgornje in spodnje 95% C.L. meje za realne dele Kj in = 2mw ter /i = mw (v 
oklepaju). Zadnji stolpec prikazuje oceno spodnje meje na skalo NF A, v primeru da je Cj ~ 1. 



omejiti preko efektov v mesanju, kjer prispevata samo z realnimi deli, niti preko CP krsitve v raz- 
padih 6—7- S7, kjer so njuni prispevki tako majhni, da omejitve niso mogoce. V primeru izpeljave 
mej iz CP krsitve v razpadih 6—7-57 prikazemo tudi projekcijo, koliko lahko pricakujemo, da se bo 
dovoljeno obmocje se skrcilo v prihodnosti, ko bo meritve izboljsal Super-Bell [124, 125]. 

Nadalje lahko rezultate za 6Ci uporabimo za analizo nekaterih opazljivk, ki se niso merjene s 
tako natancnostjo, da bi bistveno prispevale k omejevanju velikosti parametrov Kj. Za te opazljivke 
lahko preverimo napovedi odstopanja od SM, ki so se kompatibilne z izpeljanimi mejami. To 
prikazuje SI. 5.12 za razvejitvena razmerja Br[5^, — )■ /.i^/z^], Br[i? — )• K^*^^!/] ter za asimetrijo 
^fb(q^) V razpadih B^ K*l+£^. 

V prvi vrsti lahko vidimo, da sta koeficienta KLRb in i-^rr tako mocno omejena, da odstopanja 
V omenjenih razpadih ni moc pricakovati. Po drugi strani lahko opazimo, da najnovejse meritve 
razvejitvenega razmerja za Bg — )• postajajo uporabne za omejevanje ostalih parametrov. V 

primeru nadaljnjega spuscanja zgornje eksperimentalne meje bi lahko ta opazljivka postala pomem- 
ben del omejevanja parametrov kj. Analiza Ap-Qio^) razkriva, da bistvenih odstopanj od SM v tej 
opazljivki ni pricakovati. Predstavljen graf je za k^^, katerega efekti so med vecjimi. Za ostale 
parametre velja podobna ugotovitev. 
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Current and projected 95% C.L. regions 




ReM Re[yf] 



Slika 5.11: Levo: 95% C.L. (polna crta) in 68% C.L. dovoljena obmocja za parametre i^ii ■, f^'iRf 
dobljena analize mesanja mezonov B. Desno: 95% C.L. dovoljena obmocja za parametra niRb 
in krr, dobljena iz analize CP krsitve v 6 — t- 57 razapadih. Nepolne crte oznacujejo potencialno 
zozitev obmocji ob izbolsavi eksperimentalnih napak na Super-Bell. 



5.3.3 Sucnostni delezi v razpadih kvarka top 

Kot zadnje predstavimo naso analizo vpliva operatorjev (5.18) na sucnostne deleze v glavnem 
razpadnem kanalu kvarka top. V ta namen vpeljemo najsplosnejso obliko efektivnega tWh vozlisca 
preko sledece Lagrangeove funkcije 



Aff = ^0 



7^ {aLPi + anPR) - {hRLPi + biRPR) tW^, , (5.24) 



V2 

kjer anomalne sklopitve lahko povezemo s parametri Kj 

5aL = Vtbl^LL* , an = VtlK*nn , biR = -£^^t>LRt ' ^RL = ~^^^tb^LRb ■ (5-25) 

V analizo vkljucimo popravke QCD prvega reda, ki jih prikazuje SI. 5.13. Podobno kot pri analizi 
indirektnih efektov v fiziki mezonov B, bomo tudi tukaj (sprva) analizirali prispevke le enega 
operatorja naenkrat. Poleg tega se bomo omejili na realne vrednosti anomalnih sklopitev. Pod 
takimi pogoji Sa^ nima vpliva na sucnostne deleze, saj spremlja vozlisce z enako strukturo kot SM 
in se njegovi prispevki iznicijo. Za ostale tri anomalne sklopitve analiziramo najprej spremembe v 
kar prikazuje SI. 5.14. Rezultati so prikazani kot pasovi, ki ponazarjajo povecanje prispevkov 
k ob vkljucitvi popravkov QCD prvega reda. Indirektne omejitve na parametre kj, ki smo jih 
izpeljali v prejsnjem poglavju se prevedejo v 

- 0.0006 <aR< 0.003 , -0.0004 < bjiL < 0.0016 , -0.14(-0.29) < < 0.08 . (5.26) 
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Slika 5.12: Levo, Sredina: Razpon razvejitvenih razmerij, ko anomalne parametre Kj spremin- 
jamo znotraj 95% C.L. dovoljenih intervalov prikazanih v Tab. 5.2. Pikcaste crte predstavljajo la 
teoreticne negotovosti SM napovedi. Za razpad v mione prikazujemo tudi spodnjo eksperimentalno 
mejo 90% C.L. intervala [127], 95% C.L. zgornjo mejo LHCb [12?^] in najnovejso zdruzeno LHC 
meritev [129]. Desno: Afb{q'^) pas, ki ga dobimo, ko k'^^ spreminjamo znotraj 95% dovoljenih 
intervalov prikazanih v Tab. 5.2. Prikazana je tudi srednja vrednost napovedi SM (crna) in pas la 
teoreticne negotovosti (crtkano) ter eksperimentalno izmerjeni toccki z pripadajocimi napakami iz 
Ref. [1.30]. 




SUka 5.13: Feynmanova diagrama za popravke QCD prvega red k glavnem razpadnem kanalu 
kvarka top t — t- Wb. Kvadratki oznacujejo anomalno vozlisce, ki ga generira NF. Krizec oznacuje 
dodatno tocko iz katere se lahko izseva gluon. 



Iz SI. 5.14 lahko zakljucimo, da tudi ob upostevanju QCD popravkov prvega reda, vrednosti 
reda procent ali vec nikakor ne morem pricakovati. Ker je vrednost J-l izmerjena z bistveno 
boljsno natancnostjo, je vredno pogledati, ali morda izmerjena vrednost lahko sluzi za postavitev 
meje na velikost anomalnih sklopitev. Kako se J-'l spreminja v odvisnosti od anomalnih sklopitev, 
je prikazano na SI. 5.15. Vpliv korekcij QCD je v primeru J-l zanemarljiv. Anomalni sklopitvi an 
in bjiL zaradi mocnih indirektnih omejitev zopet me moreta znatno vplivati na vrednost deleza. Po 
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Slika 5.14: Odvisnost od anomalnih sklopitev, za katere privzamemo, da so realne in da je le ena 
sklopitev razlicna od nic. Levo: Odvisnost J^^ od an (modro, pikcasto), (oranzno, crtkano) 
in biR (crno, polno). Zgornja in spodnja crta pasu pripadata analizi brez in z prvim redom QCD 
popravkov. Krizec oznacuje napoved SM. Desno: Odvisnost J-+ od Blr. Crtkana crta predstavlja 
rezultat brez popravkov QCD, polna crta pa vkljucuje nace popravke prvega reda. Prikazana je 
tudi SM vrednost podana v En. (5.7) in 95% C.L. dovoljena intervala za b^j^ podana v En. (5.26). 



drugi strani manj omejeni biR lahko znatno spremeni vrednost J-l, kar je podrobneje prikazano 
na desnem grafu SI. 5.15. Vidimo, da je meja, ki jo izpeljemo iz razpadov kvarka top v tern 
primeru primer Ijiva z indirektni mejami in v prihodnosti, v kolikor bo predvidena senzitivnost 
Atlasa realizirana, lahko postane tudi dominantna. 

Ce bi analizirali poleg sucnostnih delezev se kaksno drugo opazljivko iz fizike kvarka top, bi 
lahko sprostili omejitev, s katero predpostavimo le eno nenicelno anomalno sklopitev. Ravno to so 
naredili v Ref. [139], kjer so poleg sucnostnih delezev (brez popravkov QCD) v analizo vkljucili tudi 
produkcijo enega kvarka top, ki tudi poteka preko sibke interakcije in je obcutljiva na anomalne 
tWb sklopitve. Avtorji [139] so analizirali dovoljena obmocja v ravninah parov anomalnih sklopitev. 
Zanimivo je njihove ugotovitve, ki temeljijo izkljucno na direktnih procesih, postaviti ob bok nasim 
ugotovitvam indirektnih omejitev. Ker smo z izpeljavo dovoljenih intervalov uporabili vec kot 
eno spremenljivko, lahko prikazemo analogno analizo v ravninah parov anomalnih sklopitev, kar 
prikazuje SI. 5.16. 

Sivi obmocji, dobljeni iz direktnih omejitev, imata obliko pasu, ker sta v smereh krr in kijh, 
mnogo sirsa kot obmocja, ki jih dobimo iz obravnave indirektnih omejitev, kar upravicuje nase 
sklepanje, da je v primeru NF v fiziki kvarka top potrebno vedno analize direktnih opazljivk 
postaviti ob rob indirektnim analizam v fiziki mezonov. 
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Slika 5.15: Odvisnost J-'l od anomalnih sklopitev, za katere privzamemo, da so realne in da je le 
ena sklopitev razlicna od nic. Prikazana je tudi srednja vrednost meritve iz Tevatrona (crtkano) 
in pripadajoci 95% C.L. interval. Poleg tega je prikazana tudi predvidena bodoca velikost 95% 
C.L. intervala iz Atlasa na podlago centralne vrednosti iz Tevatrona. Levo: Odvisnost J^l od an 
(modro, pikcasto), Irl (oranzno, crtkano) and 6lr (crno, polno). Desno: Odvisnost od biji 
in 95% C.L. dovoljeni intervali za Blr podani v En. (5.26). 




Slika 5.16: 95% C.L. dovoljena obmocja v razlicnih (k^, kj) ravninah. Siv pas predstavlja dovoljeno 
obmocje kot dobljeno iz direktnih omejitev v Ref. [139]. Ki so privzeto realni. Levo: k/j^ - k^l 
(polno), k'^^ (crtkano), k'^j^ (pikcasto) ravnina. Skala ujemanja je /i = 2m]y. Desno: klri, - 
i^LRt (polno), K'^jijf (crtkasto) ravnina. Skalo ujemanja spremenimo iz ^ = 2mw (ozji obmocji) na 
/i = niw (sirsi obmocji). 
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5.4 Zakljucki 

Ob koncu dobe Tevatrona smo ze globoko zakorakali v dobo LHC in lov na fiziko onkraj stan- 

dardnega modela je v polnem razmahu. Iskanjc novih dclccv nikakor ni cdini nacin, s katcrim 
LHC teoreticni fiziki osnovnih delcev prinasa nova vprasanja in odgovore. Nadejamo se razjasnitve 
problema fizike okusa, v kateri kvark top s svojo veliko maso igra vodilno vlogo. Ker LHC lahko 
smatramo kot pravo tovarno kvarka top, nam je prvic na voljo raziskava fizike kvarka top z veliko 
natancnostjo. Dolocitve parametrov in inter akcijskih struktur kvarka top nam lahko sluzi kot okno 

V svet nove fizike. V tem delu smo preucili razlicne aspekte razpadov kvarka top in raziskovali, 
kako se NF, ki smo jo parametrizirali s pomocjo efektivnih teorij, lahko v njih manifestira. 

Na eni strani smo raziskovali razpadne sirine in upostevali tudi popravke prvega reda v QCD, kar 
je smiselno, ko imamo opravka s kvarki in smo sooccni z vedno vecjo natancnostjo ckspcrimcntalnih 
meritev. Preucili smo razvcjitvcna razmcrja razpadov t — )■ 97, Z in razlicne kincmaticne opazljivke 

V tro-delcnem razpadu t — t- qi'^i^ ter glavni razpadni kanal kvarka top t — )■ Wb, kjer smo pozornost 
usmerili v sucnostne deleze bozona W, ki so obcutljivi na strukturo tWb vozlisca. 

Po drugi strani smo izpostavljali vlogo, ki jo igra top kvark v fiziki mezonov kot virtualen 
delec. Posledicno dolocene modifikacije fizike kvarka top lahko vplivajo na teoreticne napovedi 
opazljivk V fiziki mezonov. Ker za modifikacije nabitih tokov, ki vsebujejo kvarke top celostne 
analize indirektnih posledic ni moc najti v literaturi, smo podrobno preucili posledice v |AS| = 2 
in \AB\ = 1 procesih. Po pricakovanjih smo lahko na NF postavili ne-trivialnc indirektnc omejitve. 

Ne glcdc na to, ali nam LHC v razpadih kvarka top razkrijc novo fiziko ali nc, bodo bodoce mer- 
itve V fiziki kvarka top igrale pomcmbno vlogo v raziskovanju fizike okusa in grajcnja ali omejevanja 
modelov fizike onkraj standardnega modela. 
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